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1.  Einleitimg 

Konventionelle  Werkstoffe  fur  die  Luft-  und  Raumfahrttechnik  aus  Aluminium  und  Titan 
zeichnen  sich  durch  ihre  Duktilitat  aus.  Wegen  dieser  Duktilitat  ist  es  generell  moglich, 
schadenstolerant  zu  konstruieren.  Bei  der  Fertigung  von  Luftfahrtstrukturen  aus  Faserver- 
bunden  ist  es  notwendig,  ein  schadenstolerantes  („damage  tolerant")  Werkstoffverhalten  zu 
realisieren. 

Diese  Werkstoffe  versagen  haufig  „katastrophal“  ohne  die  fur  konventionelle  Werkstoffe 
bekannte  Duktilitat.  Dies  wird  durch  die  steifen  und  hochfesten  Fasem  verursacht.  Bei  der 
Entwicklung  von  Faserverbundwerkstoffen  mit  keramischer  und  damit  sproder  Matrix  wird 
deshalb  versucht,  von  Schadensmechanismen  Gebrauch  zu  machen,  die  zu  einer  Quasidukti- 
litat  fuhren. 

Eine  Moglichkeit,  dies  zu  realisieren,  ist  das  Zulassen  von  Matrixrissen,  die  einen  gewissen 
Verlust  an  Steifigkeit  verursachen,  wodurch  die  sogenannte  Quasiduktilitat  erreicht  wird.  Fur 
dieses  Verhalten  ist  es  notwendig,  dass  Matrixrisse  an  der  Faseroberflache  abzweigen.  Dies 
wird  generell  realisiert,  indem  man  die  Haftung  zwischen  Faser  und  Matrix  gering  halt. 

Kohlenstofffaserverstarkter  Kohlenstoff  (C/C)  ist  prinzipiell  keine  faserverstarkte  Kera- 
mik.  Er  wird  jedoch  generell  aufgrund  des  sproden  Verhaltens  dazu  gerechnet.  Fur  Faserver- 
bundwerkstoffe  ist  es  daher  notwendig,  Aussagen  zu  treffen,  ob  ein  sich  entwickelnder  Riss 
oder  Schaden  entweder 

•  zu  Sprodbruch  (=Faserbruch)  oder 

•  zu  quasiduktilem  Bruch  fuhrt. 

Im  Gegensatz  zum  Sprodbruch,  der  zu  einem  sofortigen  Versagen  fuhrt,  fuhrt  ein  quasidukti- 
ler  Bruch  lediglich  zur  Eigenschaftsanderungen  im  Faserverbundwerkstoff.  Hierbei  fuhren 
Schadigungen,  wie  beispielsweise  Debonding  oder  Delamination,  zu  einer  Verringerung  der 
Steifigkeit  oder  zu  einer  Verschiebung  der  Eigenfrequenzen. 

Die  Untersuchung  der  Rissausbreitung  in  kohlenstofffaserverstarktem  Kohlenstoff  ist 
Gegenstand  dieser  Arbeit 
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2.  Stand  der  Forschung 


2.1  Technische  Anwendung  von  kohlenstofffaserverstarktem  Kohlenstoff 

Kohlenstofffaserverstarkter  Kohlenstoff  (C/C)  zeichnet  sich  durch  einen  niedrigen  Warme- 
ausdehnungskoeffizienten  und  eine  hohe  Temperaturbestandigkeit  bis  3000°C  aus.  Er  ge- 
wahrleistet  eine  hohe  Steifigkeit  und  hohe  Festigkeit. 

Die  Anwendung  im  Hochtemperaturbereich  erfordert  jedoch  das  Anbringen  von  Oxidations- 
schutzschichten,  da  sich  bei  T>  600°C  eine  hohe  Oxidationsempfindlichkeit  des  Werkstoffes 
einstellt.  Als  Beschichtungsmaterial  kommt  hierbei  folgendes  Material  zum  Einsatz: 
T<1500°C  SiC  oder  Si3N4 

1500°C<T<1800°C  Borsilikate 

T>1800°C  Iridium  und  Platin-Legierungen 

Sein  niedriges  spezifisches  Gewicht  erweitert  sein  Einsatzspektrum  auf  Raumfahrtstrukturen. 
Die  Abschleifbestandigkeit  (Abrasiwerhalten)  fuhrt  zu  einer  Verwendung  als  Bremsschei- 
benwerkstoff.  Fur  das  Kampfflugzeug  F-16  wurden  beispielsweise  2D-C/C-Bremsscheiben 
entwickelt.  Damit  lassen  sich  bis  zu  2000  Landungen  durchfuhren.  Eine  Einfuhrung  von  3D- 
C/C  wird  den  Einsatz  auf  3000  Landungen  erweitem.  Im  Vergleich  dazu  werden  Wechsel 
der  Bremsscheiben  bereits  nach  180  Landungen  notwendig,  wenn  diese  aus  konventionellen 
Werkstoffen,  wie  beispielsweise  beim  Kampfflugzeug  F-4F  Phantom,  gefertigt  sind. 


Abb.  2-la:  Bremsscheiben  fur  F-16 
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2.2  Eigenschaften  von  kohlenstofffaserverstarktem  Kohlenstoff 

Kohlenstofffaserverstarkter  Kohlenstoff  (C/C)  wird  aufgrund  des  sproden  Verhaltens  zu 
den  faserverstarkten  Keramiken  gerechnet.  Um  „damage  tolerantes“  Werkstoffverhalten  zu 
realisieren,  werden  Matrixrisse  zugelassen.  Dieser  Schadensmechanismus  soli  trotz  der  spro¬ 
den  Matrix  zu  einer  Quasiduktilitat  fuhren. 

Der  untersuchte  C/C-Werkstoff  wurde  uber  die  sog.  Polymer-Route  mit  den  unterschiedli- 
chen  Prozessschritten  gemaB  Abb.  2-lb  hergestellt  [1], 


n 

Kohlenstofffasem 

I'-lV;: 

J 

V 


Phenolharz 


Impragnieren 


Harten,  180°C 


Karbonisieren  <1000°C 
(starke  Schrumpfung) 


Graphitieren  ca.  2000°C 
(Strukturanderung) 


Abb.2-lb:  Herstellung  von  C/C 


Als  Ausgangsmaterial  werden  Kohlenstofffasem  auf  Polyacrylnitril-Basis  (PAN)  verwendet. 
Die  Herstellung  der  Kohlenstofffasem,  die  mit  Hilfe  der  Polymer-Route  gemaB  Abb.  2-lb  zu 
C/C  weiterverarbeitet  werden,  erfolgt  durch  das  Verspinnen  des  PAN  aus  dem  schmelzfliis- 
sigen  Zustand.  Nach  zwei  Zipperreaktionen  bei  1 80°C  und  250°C  folgt  eine  Oxidationsreak- 
tion  bei  350°C.  Durch  die  anschlieBende  Karbonisierung  wird  ein  Schichtpolymer  unter 
Freiwerden  von  gasformigen  Reaktionsprodukten  erzeugt.  AnschlieBend  erfolgt  die  Graphi- 
tierung  der  Kohlenstofffasem.  Um  die  Faserfestigkeit  zu  steigem,  miissen  die  hexagonalen 
Schichtebenen  des  Kohlenstoffs  (C-  6er-Ringe)  moglichst  parallel  liegen.  Ein  Strecken  der 
Faser  unter  Anbringen  von  Zugspannung  gewahrleistet  die  erforderlichen  Reorientie- 
mngsprozesse  wahrend  des  „Processings“. 
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Die  so  gefertigten  HT-Fasem  werden  mit  dem  Polymer  Phenolharz  impragniert.  Das  Phe- 
nolharz  wird  bei  180°C  ausgehartet  und  wird  durch  eine  anschlieBende  Warmebehandlung 
bei  Temperaturen  bis  1000°C  karbonisiert.  Bei  diesem  Prozess,  der  unter  Ausschluss  von 
Sauerstoff  stattfmdet,  wird  die  Polymermatrix  in  eine  Glaskohlenstoffmatrix  umgewandelt. 
Dies  geht  mit  einer  erheblichen  Schrumpfung  der  Matrix  und  unter  Entwicklung  von  Pyroly- 
segasen  einher.  Findet  diese  Schrumpfung  ungehindert  statt  (die  hochmodulige  C-Faser  ver- 
hindert  diese  stark),  dann  betragt  die  Schrumpfung  linear  20%.  Dies  bedeutet,  dass  als  Folge 
der  Karbonisierung  das  Polymer  auf  50%  seines  Ausgangsvolumens  schrumpft. 

Als  Folge  der  Behinderung  des  Schrumpfens  im  Faserverbund  und  der  Pyrolysegasentwick- 
lung  treten  bei  der  Karbonisierung  Risse  und  Hohlraume  auf  (Abb.  2-2).  Diese  kdnnen  in 
nachfolgenden  Reimpragnierungsschritten  (teilweise)  geschlossen  werden  (Abb.  2-3),  sodass 
die  Porositat  im  Endprodukt  zwischen  5  und  15%  betragt.  In  einem  letzten  Schritt  wird  der 
karbonisierte  C/C-Werkstoff  bei  Temperaturen  oberhalb  von  ca.  2000°C  behandelt,  haupt- 
sachlich,  um  die  Glaskohlenstoffmatrix  (teilweise)  in  eine  graphitierte  Matrix  mit  verbesser- 
ten  mechanischen  Eigenschaften  umzuwandeln. 


Matrixrisse  in  Querbundeln  Poren  in  Querbiindeln 
Abb.2-2  Poren  und  Risse  in  2D  C/C 
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Abb. 2-3:  Reimpragnierungsschritte 


4000 


2000  - 


1000  H 


• 

4 

Zusatzliche 

\ 

Oliihbehandlung 

Karbonisierung 

j^X  1000  2000  3000 

X*  Behandlungstemperatur  [°C] 


Faserbsvch  nach  der  Karbonisierung 


Abb. 2-4:  Festigkeit  in  Abhangigkeit  von  der  Behandlungstemperatur 
Phenolharz-impragnierte  3K-HTA-Faserbiindel  [1] 
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Charakteristisch  fur  die  Festigkeitsentwicklung  von  C/C  als  Funktion  der  Behandlungstem- 
peratur  (Abb.2-4)  ist  eine  Abnahme  der  Festigkeit  auf  ein  Minimum  nach  der  Karbonisie- 
rung  und  eine  teilweise  Erholung  der  Festigkeit  nach  einer  weiteren  thermischen  Behandlung 
oberhalb  von  1000°C  [1], 

Das  Aussehen  der  Bruchflache  des  karbonisierten  Materials  ist  glatt  (Sprodbruch  als  Folge 
von  Penetration  der  Risse  in  die  Fasem)  wahrend  nach  einer  zusatzlichen  Warmebehandlung 
oberhalb  von  1000°C  ein  Ubergang  zu  einem  Bruchbild  mit  Rissabzweigungsmerkmalen 
auftritt. 

Die  Untersuchung  des  Einflusses  der  Schrumpfspannungen  in  der  Matrix,  der  thermischen 
Spannungen  infolge  unterschiedlicher  Warmeausdehnungskoeffizienten  und  der  Eigen- 
schaftsanderung  in  Faser  und  Matrix  nach  der  zusatzlichen  Warmebehandlung  auf  das  Riss- 
fortschrittsverhalten  sind  Gegenstand  dieser  Arbeit 
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2.3  Die  Energiefreisetzungsrate  fur  die  Rissausbreitung 

Das  Rissausbreitungsverhalten  von  isotropen  homogenen  Werkstoffen  kann  mit  Hilfe  der 
Energiefreisetzungsrate  G  beschrieben  werden.  Dazu  wird  eine  unendlich  groBe  Platte  mit 
einem  Riss  betrachtet. 


JjC 

a 


Dicke“l 


Abb. 2-5:  Unendliche  Platte  mit  Riss 


GemaB  Griffith  [2]  andert  sich  die  Verformungsenergie  einer  angerissenen  zugbeanspruchten 
Scheibe  von  unendlicher  Ausdehnung  gemaB  Abb.2-5  um  den  Betrag  AUe  mit: 


A  U 


e 


2  2 

mo 


E 


(2-1) 


worin  a  die  Risslange,  a  die  angelegte  Spannung  und  E  den  Elastizitatsmodul  des  Werkstof- 
fes  darstellt.  Zur  Bildung  der  erweiterten  Rissflachen  muss  die  Oberflachenenergie  [3]  von 

^0=4  a7s  (2-2) 


aufgebracht  werden.  Hierbei  ist  ys  die  spezifische  Oberflachenenergie. 
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Der  kritische  Fall  der  instabilen  Rissausbreitung  wird  dann  erfolgen,  wenn 


d(AUe)  =  d(AU0) 

d(2a)  2  (da) 


Hierbei  ist  der  Term 


d(AUe)_c 

d(2a) 


(2-4) 


die  Energiefreisetzungsrate.  Der  Term 

d(AU0)_R 

d(2a) 


(2-5) 


stellt  die  Bruchenergie  dar  und  ist  der  Risswiderstand,  den  der  Werkstoff  der  Ausbreitung 
eines  Risses  entgegensetzt.  Beim  ideal  sproden  Bruch  fuhrt  eine  infinitesimale  kleine  Riss- 
verlangerung  zur  Instability.  Fur  diesen  Fall  ist  die  Energiefreisetzungsrate  G  gleich  der  kri- 
tischen  Energiefreisetzungsrate  Gc.  Bei  einem  sich  plastisch  verformenden  Werkstoff 
muss  die  angelegte  Belastung  erhoht  werden,  um  eine  einmal  eingeleitete  Rissausbreitung 
aufrechtzuerhalten,  d.h.  der  Risswiderstand  steigt  mit  der  Rissverlangerung  Aa  an. 


dG  dR 
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Tragt  man  nun  R  als  Funktion  von  Aa  gemaB  Abb.  2-6  auf,  so  erhalt  man  die  entsprechende 
R-Kurve  oder  auch  Risswiderstandskurve  genannt.  Solange  die  Energiefreisetzungsrate 
schwacher  mit  a  zunimmt  als  der  Risswiderstand  bleibt  der  Riss  stabil: 


dR  dG 
da  da 

Der  Ubergang  zur  instabilen  Ausbreitung  findet  an  folgendem  Punkt  statt: 

dR  _  dG 
da  da 


(2.6) 


(2.7) 


Das  Rissausbreitungsverhalten  in  Faserverbundwerkstoffen  ist  auBerst  komplex.  In  faserver- 
starkten  Werkstoffen  tritt  generell  erste  Rissbildung  in  der  Matrix  auf. 

Die  Ermittlung  der  Energiefreisetzungsrate  in  kohlenstofffaserverstarktem  Kohlenstoff  muss 
folglich  naher  untersucht  werden,  um  Aussagen  iiber  das  in  der  Einleitung  geschilderte  Ver- 
halten,  ob  der  Riss  in  die  Faser  eindringt  oder  an  der  Faser-Matrix-Grenzflache  abzweigt,  zu 
treffen. 

2.4  Energiekriterium  fur  das  Bifurkationsproblem  der  Rissabzweigung  bzw.  Riss- 
penetration 

2.4.1  Arten  des  Rissfortschrittes 


Ein  Matrixriss,  der  auf  ein  Interface  trifft,  das  sich  zwischen  Faser  und  Matrixwerkstoff  be- 
findet,  kann  auf  zwei  Arten  weiterlaufen: 

a)  Rissabzeigung  am  Interface  um  90°  („Single  Deflection41)  nach  einer  Seite 

b)  Rissabzweigung  am  Interface  um  90°  nach  zwei  Seiten  (..Double  Deflection11) 

c)  Risspenetration  am  Interface  (..Penetration11) 

Dieser  Wettbewerb  zwischen  den  zwei  Rissfortschrittsarten  wird  von  unterschiedlichen  Fak- 
toren  beeinflusst.  Eigenspannungen  infolge  von  herstellungsbedingten  Schrumpfprozessen, 
unterschiedliche  Ausdehnungskoeffizienten,  Orientierungseffekte  in  der  Faser  und  Matrix 
(Graphitieren)  konnen  hier  wesentliche  Einflusse  ausiiben.  Abb.  2-7a  bis  Abb.  2-7c  zeigt  die 
unterschiedlichen  Rissfortschrittsmoglichkeiten: 
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Abb.  2-7c:  Penetration  in  die  Faser  („Penetration“) 
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2.4.2  Energiekriterium 

Um  Aussagen  uber  das  o.a.  Rissfortschrittsverhalten  zu  erlangen,  ist  es  notwendig,  ein  Krite- 
rium  fiir  die  Risserweiterung  zu  entwickeln.  He  und  Hutchinson  [5]  haben  hierzu  ein  Ener¬ 
giekriterium  aufgestellt,  anhand  dessen  Aussagen  iiber  den  Rissfortschritt  gemacht  werden 
konnen.  GemaB  He  und  Hutchinson  [5,6]  wird  der  Matrixriss  in  das  Interface  abzweigen 
(„Deflection“),  wenn  folgende  Bedingung  erfullt  ist: 


wobei  Gc'  die  Bruchenergie  des  Interfaces  und  Gcf  die  Mode  I  Bruchenergie  der  Faser 

ist.  Gd  stellt  die  Energiefreisetzungsrate  fur  einen  abzweigenden  Riss  bzw.  Gp  die  Energie- 
freisetzungsrate  fur  einen  Matrixriss,  der  in  die  Faser  eindringt,  dar. 

Der  nach  einer  Seite  hin  abzweigende  Riss  („Single  Deflection11)  ist  im  Regelfalle  energe- 
tisch  giinstiger  als  der  nach  beiden  Seiten  hin  abzweigende  Riss  („Double  Deflection"),  weil 
er  zu  den  groBten  Werten  von  Gd  fiihrt. 

Der  Riss  wird  in  die  Faser  eindringen  (Penetration),  wenn  die  Ungleichung  nicht  erfullt  wird. 
Die  einzelnen  Komponenten  der  Ungleichung  sind  fur  die  entsprechende  Faser-/Matrix- 
Zusammenstellungen  unter  den  gegebenen  Belastungen  zu  ermitteln. 

Im  Verlauf  dieser  Arbeit  wird  anhand  verschiedener  Verfahren  fur  den  kohlenstofffaserver- 
starkten  Kohlenstoff  (C/C)  in  Abhangigkeit  von  der  Behandlungstemperatur  die  Bruchener¬ 
gie  Gc'  des  Interfaces  bzw.  die  Bruchenergie  der  Faser  Gcf  ermittelt  (siehe  dazu  auch  [50]). 
Um  die  Materialeigenschaften  des  Interfaces  zu  bestimmen,  wird  ein  Druckscherversuch 
ausgewertet.  Zur  Ermittlung  der  Faserbruchenergie  Gcf  werden  Ergebnisse  aus  der  Literatur 
nach  Ochiai  und  Osamura  [7]  bzw.  Honjo  und  Shindo  [8]  herangezogen. 
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Die  Energiefreisetzungsraten  fur  abzweigende  Risse  Gd  und  Penetrationsrisse  Gp  werden  auf 
unterschiedliche  Arten  bestiramt.  Zunachst  werden  einige  Ergebnisse  aus  der  Literatur  be- 
schrieben.  Diese  Ergebnisse  werden  fur  unterschiedliche  Voraussetzungen  analytisch  be- 
stimmt  und  sind  somit  nur  fur  bestimmte  Falle  giiltig.  1 

Nach  ausfuhrlicher  Darstellung  der  diversen  Theorien  und  analytischer  Ansatze  aus  der 
Literatur  werden  diese  miteinander  verglichen.  Die  Ergebnisse  dienen  als  Ausgangswerte  fur 
die  Ermittlung  der  Werte  fur  den  kohlenstofffaserverstarkten  Kohlenstoff.  Um  die  exakten 
Werte  fur  Gd  und  Gp  zu  erhalten,  wird  schlieBlich  eine  Finite  Element  Analyse  durchgefuhrt. 
Damit  wird  das  in  der  Einleitung  beschriebene  Problem  charakterisiert  bzw.  gelost. 


1  So  ergibt  sich  gemaB  einer  analytischen  Losung  nach  He  und  Hutchinson  [5,6]  fur  das  Verhaltnis 
Gd/Gp  ungefahr  14,  wenn  Faser  und  Matrix  aus  dem  gleichen  homogen  isotropen  Material  bestehen. 
Wenn  der  Elastizitatsmodul  der  Faser  den  dreifachen  Wert  der  des  Matrixwerkstoffes  ist,  wachst  die- 
ser  Wert  auf  14  an. 
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2.5  Gegeniiberstelhmg  von  verschiedenen  Theorien  zur  analytischen  Bestimmung 
von  Gn/Gn 


Im  folgenden  werden  nun  Theorien  gegeniibergestellt,  die  auf  analytischem  Wege  das  Ver- 

haltnis  der  Energiefreisetzungsrate  Gd/Gp  bestimmen.  In  den  einzelnen  Abschnitten  werden 

diese  Theorien  vorgestellt  und  zusammengefasst. 

Dabei  wird  nun  auf  folgende  theoretische  Ansatze  eingegangen: 

1)  He  und  Hutchinson  [5]  beschreiben  die  Rissabzweigung  an  einem  Interface  zwischen 
unterschiedlichen  elastischen  Werkstoffen.  Hierbei  wird  sowohl  von  einer  isotropen  Fa- 
ser  als  auch  von  einer  isotropen  Matrix  ausgegangen. 

2)  Tullock,  Reimanis,  Graham  und  Petrovic  [9]  diskutieren  Rissabzweigung  und  Risspe- 
netration,  wenn  Faser  und  Matrix  aus  unterschiedlichen  Materialien  bestehen.  Dabei 
wird  eine  Boundary  Element  Methode  genutzt.  Die  Ergebnisse  werden  anschlieBend  mit 
denen  von  He  und  Hutchinson  verglichen. 

3)  He,  Evans  und  Hutchinson  [10]  ermitteln  den  Einfluss  von  Restspannungen  auf  den 
Rissfortschritt  in  isotropen  Werkstoffkombinationen. 

4)  Martinez  und  Gupta  [1 1]  erarbeiten  eine  Theorie  zur  Ermittlung  eines  Energiekriteriums 
fur  Rissabzweigung  am  Interface  fur  Faserverbundwerkstoffe,  die  aus  orthotropen  Fa- 
sem  und  einer  orthotropen  Matrix  bestehen. 

5)  Charalambides  [12,13]  untersucht  den  Einfluss  von  thermischen  und  mechanischen  Las- 
ten  auf  die  Energiefreisetzungsrate  von  Faserverbunden  mit  sproder  Matrix 
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2.6  Energiekriterium  fur  Faserverbundwerkstoffe  mit  isotropen  Materialkompo- 
nenten  nach  He  und  Hutchinson  [51 


2.6.1  Rissabzweigung/Risspenetration 

Im  folgenden  wird  nun  analysiert,  welche  Faktoren  einen  Einfluss  auf  die  Risserweiterung 
(Penetration  oder  Rissabzweigung)  haben.  Hierbei  ist  besonders  der  Verlauf  der  Energiefrei- 
setzungsrate  mit  zunehmender  Risslange  fur  die  beiden  Risserweiterungsarten  interessant. 
Die  Entscheidung,  ob  ein  Riss  beim  Auftreffen  auf  das  Interface  in  die  Faser  eindringt  oder 
ob  er  am  Interface  abzweigt,  wird  diskutiert.  Im  folgenden  werden  die  unterschiedlichen  Fal- 
le  von  Abb.  2-7a  bis  2-7c  analysiert  und  diskutiert.  Grundlage  der  ersten  Problemstellung  ist 
ein  halbunendlicher  Riss.  Dieser  verlauft  in  der  Matrix  und  trifft  in  einem  rechten  Winkel  auf 
das  Interface.  Durchdringung  tritt  durch  eine  Erweiterung  des  Risses  in  der  gleichen  Ebene 
auf  (Abb.  2-7c),  wahrend  Abzweigung  sowohl  einseitig  (Abb.  2-7a)  als  auch  zweiseitig 
(Abb.  2-7b)  stattfmden  kann. 

Die  isotrope  Faser  bzw.  die  isotrope  Matrix  wird  durch  den  Schubmodul  g,  und  die  Querkon- 
traktionszahl  V,  charakterisiert,  wobei  der  Index  i  fur  1  und  damit  fur  die  Faser  bzw.  2  fur  die 
Matrix  steht.  Fur  den  ebenen  Dehnungszustand  hangt  die  Losung  von  den  Materialparame- 
tem,  den  sogenannten  Dundursparameter  a  und  P  [14],  ab: 


a_  v2)-/x2(l- v,) 

/t,(l-v2)  +  /t2(l-v]) 

(2-9) 

;  Ju1(l-v2)  +  ^2(l-v1) 

(2-10) 

Der  Parameter  a  wird  in  der  englischsprachigen  Literatur  als  „elastic  mismatch14  bezeichnet. 

2.6.2  Herleitung  der  Energiefreisetzungsrate  ernes  Matrixrisses  fur  Abzweigen  am 
Interface  oder  Eindringen  in  die  Faser 

Als  erster  und  haufigster  Fall  soil  die  Energiefreisetzungsrate  fur  das  Eindringen  eines  Mat¬ 
rixrisses  in  die  Faser  (Abb.2-7c)  bzw.  das  Abzweigen  um  90°  (einfach  gemaB  Abb.  2-7a  oder 
doppelt  gemaB  Abb.2-7b)  theoretisch  hergeleitet  werden. 

Hierbei  wird  von  einem  halbunendlichen  Riss,  dessen  Spitze  am  Interface  anstoBt,  ausge- 
gangen.  Die  Risserweiterung  a  ist  somit  gleich  0. 
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Eine  symmetrische  Last  wird  in  Bezug  auf  die  Rissebene  angenommen.  Die  Spannungsver- 
teilung  an  der  Rissspitze  in  der  Faser  (Material  1)  wird  charakterisiert  durch 

°I,(0,y)  =  k,(2xyrl  (2-11) 

Hierbei  ist  der  Amplitudenfaktor  ^  proportional  zu  der  angenommenen  Last. 


Abb.:  2-8  Spannungen  an  der  Rissspitze 

Der  Exponent  der  Spannungssingularitat  X  ist  reell  und  hangt  von  den  elastischen  Konstan- 
ten  a  und  [3  gemaB  Zak  und  Williams  [15]  wie  folgt  ab: 


cosA/r  = 


2(0  -  a)  n  ^i.a  +  P2 

1 TT(  )  TT 


(2-12) 


Anhand  dieser  Gleichung  ergibt  sich  unter  Vemachlassigung  der  Orthotropie  der  Kohlen- 
stofffaser  mit  a=0,23  ein  X  von  0,49.  Fur  Faserverbundwerkstoffe,  bei  denen  Faser  und  Mat¬ 
rix  aus  identischem  Material  bestehen,  ergibt  sich  mit  a=0  ein  X  von  0,5. 

Es  wird  angenommen,  dass  der  Bezugsriss,  wie  in  Abb.  2-7a  bis  2-7c  gezeigt,  auf  drei  ver- 
schiedene  Arten  fortschreitet: 

1)  Abzweigen  nach  Auftreffen  auf  das  Interface  fur  einfaches  Abzweigen  gemaB  Abb.  2- 
7a  und  fur  doppeltes  Abzweigen  gemaB  Abb.  2-7b 

2)  DurchstoBen  des  Interfaces  (Penetration)  gemaB  Abb.  2-7c 
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Aus  der  Gleichung.  2-12  ist  ersichtlich,  dass  die  GroBe  der  Spannungssingularitat  von  der 
Steifigkeit  der  beteiligten  Materialien  abhangt.  Wenn  ein  Riss  von  einem  Medium  in  ein 
zweites  steiferes  eindringt  (a>0),  sinkt  die  Spannungssingularitat.  Im  umgekehrten  Fall 
steigt  die  GroBe  der  Spannungssingularitat  an.  Die  nachsten  Abschnitte  behandeln  die  Er- 
mittlung  der  Energiefreisetzungsraten  fiir  die  o.a.  drei  Falle. 

2.6.3  Eindringen  in  die  Faser  nach  Durchstofien  des  Interfaces 

Im  Fall  des  DurchstoBens  des  Interfaces  gemaB  Abb.  2-7c  wird  der  Spannungszustand  an 
der  fortschreitenden  Rissspitze  mit  einfacher  Mode  I  Bedingung  beschrieben.  Der  Span- 
nungsintensitatsfaktor  hangt  daher  wie  folgt  von  ki  und  a  ab: 

K,  =c(a,p)k,am-1  (2-13) 


wobei  c  dimensionslos  ist.  Die  Energiefreisetzungsrate  ergibt  sich  zu: 


G 


p 


1-v, 

2/i, 


K 


LlLc2kia'~2X 

111, 


(2-14) 


2.6.4  Energiefreisetzungsrate  bei  Abzweigen  des  Matrixrisses  am  Interface  (Deflecti¬ 
on) 

Im  Fall  des  Abzweigen  am  Interface  ergibt  sich  fiir  Abb.  2-7a  und  2-7b  die  Spannungsver- 
teilung  am  Interface  unmittelbar  vor  der  Spitze  nach  Rice  [16]  wie  folgt: 

Oyy  (x,0)  +  i(J  xy  (x,0)  =  (K,+iK2)(27tryU2rl£  (2-15) 


mit  r=x-a,  i=V-l  und 


1  1-/3 

£  - — In — — 

2k  1  +  p 


(2-16) 


In  diesem  Fall  fiihren  weitere  Uberlegungen  zu: 


Kx  +iK2  =  k,aV2  [d(a,p)aiE  +e(a,P)a  '*] 


(2-17) 
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Hierbei  sind  d  und  e  dimensionslose  komplexe  Funktionen  der  elastischen  Konstanten  a  und 
P,  die  sich  aus  der  Losung  der  Integralgleichung  ergeben.  Auf  diese  soil  hierbei  jedoch  nicht 
naher  eingegangen  werden.  Die  Energiefreisetzungsrate  ftir  Rissabzweigung  (deflection)  ist: 

Gd  =  [(1  —  v,) / /i!  +(1- v2)/ (i2](Kf  +  AT22)/(4cosh2  kb)  (2-18) 

K\  +  K\  =  kfal~n[\d\2  +  |ef  +  2  Re(de)\  (2-19) 

2.6.5  Grenzwertbetrachtung  in  der  Nahe  des  Interfaces 


In  jedem  der  drei  Falle  (Abb.  2-7a  bis  2-7c)  geht  die  Energiefreisetzungsrate2  mit  a->0  in 
Abhangigkeit  von  dem  Exponent  der  Spannungssingularitat  (A,  kleiner  oder  grofier  als  A  ) 
gegen  0  bzw.  gegen  unendlich.  Die  Abhangigkeit  der  Energiefreisetzungsrate  G  von  a  ist 
sehr  gering,  da  A  nur  leicht  von  A  abweicht  (siehe  Gleichung  2-12).  Lediglich  fur  kleine 
Werte  von  a  (  a<-0,7  )  weicht  A  bedeutend  von  A  ab.  Im  Bereich  des  Interfaces  konvergiert 
die  Energiefreisetzungsrate  nach  dieser  Betrachtung  folglich  nicht  in  einem  bestimmten 
Wert,  wenn  A  stark  von  A  abweicht. 

2.6.6  Ermitteln  des  Verhaltnisses  Ga/Gn  aus  den  beiden  Energiefreisetzungsraten 

Um  im  Laufe  der  Arbeit  Aussagen  zu  machen,  welchem  Verlauf  der  Matrixriss  nach  Errei- 
chen  des  Interfaces  folgt,  ist  es  notwendig,  das  Verhaltnis  aus  den  entsprechenden  Energie¬ 
freisetzungsraten  zu  bestimmen.  Ermittelt  man  dieses  nun  unter  o.a.  Vorraussetzungen,  so 
ergibt  sich,  dass  das  Verhaltnis  Gd/Gp  unabhangig  von  a  (  und  k,)  ist: 

TT = [Tr^rlMf  +  M!  +  ^  (2-2°) 

Gp  (1-a) 


In  der  spateren  FEM-Berechnung  ist  dies  zu  berucksichtigen.  Die  Modellierung  eines  Rissfortschrittes  in  der 
Nahe  des  Interfaces  wird  somit  ebenso  keinen  absoluten  bzw.  mathematischen  Grenzwert  fur  die  Energiefreiset¬ 
zungsrate  Gd  oder  Gp  bringen.  Da  gcmaG  Gleichung  2-12  der  Exponent  X  nur  unwesentlich  von  'A  fur  das  betrach- 
tete  Material  abweicht,  kann  aus  physikalischen  Uberlegungen  heraus  [17]  trotzdem  bei  einer  diskreten  Annahe- 
rung  an  das  Interface  ein  Wert  fur  Gd  und  Gp  festgelegt  werden.  Aus  diesen  beiden  Werten  ist  es  dann  zulassig 
das  Verhaltnis  Gd/Gp  zu  bestimmen.  Im  folgenden  wird  zusatzlich  auch  mathematisch  gezeigt,  dass  fur  das  Ver¬ 
haltnis  der  Energiefreisetzungsraten  Gd/Gp  ein  Grenzwert  existiert,  da  hierbei  keine  Abhangigkeit  von  a  bzw.  X 
vorhanden  ist. 
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Mit  Hilfe  dieses  Verhaltnisses  (Gleichung  2-20)  kann  nun  bestimmt  werden,  ob  ein  Matrix- 
riss  sich  durch  Abzweigen  an  der  Faser-Matrix-Grenzflache  oder  durch  Eindringen  in  die 
Faser  erweitert.  Die  Funktion  c(a,P)  fur  den  Fall  des  Eindringens  und  d(a,(3)  und  e(a,p)  fur 
die  zwei  Falle  des  abzweigenden  Risses  (deflection)  ermitteln  sich  aus  der  Losung  von  Integ- 
ralgleichungen. 

Die  Abb.  2-9  zeigt  das  Verhaltnis  Gd/Gp  in  Abhangigkeit  von  a  fur  P=  0.  Positive  a-Werte 
bedeuten  eine  hohere  Steifigkeit  der  Faser  gegeniiber  der  Steifigkeit  der  Matrix. 


Abb. 2-9:  Verhaltnis  der  Energiefreisetzungsraten  in  Abhangigkeit  von  a 

Da  der  Einfluss  von  p  auf  das  Verhaltnis  der  Energiefreisetzungsraten  nicht  besonders  groB 
ist,  wurde  P=0  gesetzt,  was  fur  die  meisten  Materialkombinationen  angenommen  werden 
kann.  Die  roten  Markierungen  geben  an,  welcher  Wert  des  Energiefreisetzungsratenverhalt- 
nisses  Gd/Gp  fur  den  kohlenstofffaserverstarkten  Kohlenstoff  unter  Vemachlassigung  der 
Faserorthotropie  ermittelt  werden  kann. 
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2.6.7  Zusammenfassung  der  Theorie  von  He  and  Hutchinson 

He  und  Hutchinson  [5]  ermitteln  mit  Hilfe  eines  theoretischen  Ansatzes  heraus  Werte  fur  das 
Energiefreisetzungsratenverhaltnis  Gd/Gp  .  Dieses  Verhaltnis  kann  jedoch  nur  fur  Material- 
kombinationen  aus  isotropen  Werkstoffen,  wie  z.B.  Glasfaser/Epoxy,  verwendet  werden. 
Lediglich  unter  Vemachlassigung  der  Orthotropie  der  Kohlenstofffaser  konnen  Anhaltswerte 
fur  das  Energiefreisetzungsratenverhaltnis  von  kohlenstofffaserverstarktem  Kohlenstoff  ge- 
wonnen  werden. 

Anhand  der  Grenzwertbetrachtung  gemaB  Abschnitt  2.6.5  muss  man  fur  die  spatere  FEM- 
Analyse  berucksichtigen,  dass  es  bei  der  Annaherung  des  Risses  an  das  Interface  keinen 
Grenzwert  fur  die  Energiefreisetzungsrate  (auBer  null  und  unendlich)  gibt.  Lediglich  fur  den 
Quotienten  besteht  ein  Grenzwert.  Die  Annaherung  des  Risses  mit  beliebigem  inkrementel- 
lem  Rissfortschritt  Aa  ist  somit  unnotig.  Es  muss  daher  bei  der  Wahl  des  inkrementellen 
Rissfortschrittes  im  Bereich  des  Interfaces  ein  Kompromiss  beziiglich  physikalisch  sinnvol- 
ler  GroBen  und  numerischer  Genauigkeit  bestimmt  werden. 


19 


Rissausbreitung  in  kohlenstofffaserverstarktem  Kohlenstoff  (C/C) 
Thorsten  Weber 


DLR  -  Institut  fur  Werkstoff-Forschung 


2.7  Ermitteln  der  Energiefreisetzungsrate  mit  Hilfe  einer  Boundary  Element  Me- 
thode  nach  Tullock.  Reimanis.  Graham  und  Petrovic  [91 

2.7.1  Allgemeines 

Die  bisherigen  analytischen  Techniken  zielten  auf  Studien  ab,  die  die  Interaktion  zwischen 
Riss  und  Interface  untersuchten,  wobei  der  Riss  als  eine  kontinuierliche  Verteilung  von  Ver- 
schiebungen  entlang  einer  halbunendlichen  Strecke  gesehen  wurde.  Die  analytischen  Metho- 
den,  die  in  den  bisher  vorgestellten  Verfahren  benutzt  wurden,  waren  bezuglich  der  gewahl- 
ten  Modellgeometrie  und  der  Modellierung  von  infinitesimalen  Rissen  stark  eingeschrankt. 
Dadurch  war  es  nicht  moglich,  geometrische  Einfliisse  beispielsweise  die  Verstarkungsgeo- 
metrie,  die  Orientierung  von  Fasem  oder  die  Variation  des  Faservolumengehaltes  einfliefien 
zu  lassen.  Dies  fuhrt  dazu,  dass  es  einer  besonderen  Analyse  bedarf,  mit  der  schwierige  Mo- 
dellgeometrien  und  andere  Randbedingungen  beriicksichtigt  werden  konnen.  Die  o.a.  Auto- 
ren  wahlten  hierzu  im  Gegensatz  zu  meiner  noch  ausstehenden  Finiten  Element  Berechnung 
(FEM)  als  Ansatz  die  Boundary  Element  Methode  (BEM),  die  sich  fur  elastisches  Material 
eignet. 

Die  BEM  formuliert  die  Ausgangsdifferentialgleichungen  in  ein  System  um,  das  aus  Inter- 
gralgleichungen  besteht.  Dieses  System  ist  auf  dem  Rand  definiert.  Daraus  ergibt  sich  ein 
System  von  niedriger  Ordnung,  wobei  die  numerische  Losung  und  die  Vemetzung  mit  ge- 
ringem  Aufwand  verbunden  ist. 

Im  folgenden  wird  die  Boundary  Element  Methode  (BEM)  genutzt,  um  Aussagen  fiber  das 
Rissfortschrittsverhalten  bezuglich  Risserweiterung  in  die  Faser  oder  Risserweiterung  ins  In¬ 
terface  gegeben. 

2.7.2  Aufstellen  der  Grundgleichungen 

Ausgehend  von  einer  statischen  und  isothermen  Deformation  eines  elastischen  Kontinuums 
D  mit  dem  Rand  S  und  den  Lame-Konstanten  A  und  (X  werden  die  Naviergleichungen  als  re- 
levante  Feldgleichungen  aufgestellt.  Sie  werden  wie  folgt  bestimmt: 

Vcr  =  0;<t  =  AVwl +  /t(Vw  +  D  (2-21) 

wobei  o  der  Cauchyspannungstensor  und  u  der  Verschiebungsvektor  ist.  1  ist  der  Identi- 
tatstensor  und  der  Index  T  bestimmt  die  Transponierte. 
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2.7.3  Die  Boundary  Element  Methode 

Im  folgenden  soil  nun  nur  kurz  auf  die  BEM  eingegangen  werden,  da  diese  nur  als  Vergleich 
bzw.  Bestatigung  der  von  He  und  Hutchinson  [5]  erzeugten  Ergebnisse  herangezogen  wird. 
Eine  komplette  Ableitung  ist  daher  hier  nicht  erforderlich.  Es  wird  nur  eine  kurze  Zusam- 
menfassung  gegeben.  Ausgehend  von  Betti’s  reziprokem  Arbeitstheorem  kann  mit  x  — »  S 
das  Verschiebungsfeld  wie  folgt  ermittelt  werden: 

cv  ( x)Uj  (x)  =  £  u*j  ( x ,  y)tj  ( y)dS  (y)  f*  (x,  y)uj  ( y)dS  ( y )  (2-22) 

wobei  ye  S,  Uj(y)  und  tj(y)  die  j-Komponenten  der  Verschiebungen  und  Haftungen  bei  y  sind. 
Uij*(x,y)  ist  die  i-Komponente  der  Verschiebungen  bei  x,  was  durch  eine  Einheitskraft  in  j- 
Richtung  bei  y  verursacht  wird  und  ty*(x,y)  ist  die  damit  zugeordnete  Haftung.  cy(x)  ist 
abhangig  von  der  Lage  von  x  und  der  lokalen  Geometrie  und  wird  wie  folgt  bestimmt: 


JCG  S 
xe  D-S 


(2-23) 


Gleichung  2-22  und  2-23  sind  als  Basis  fur  die  sogenannten  direkten  Boundary  Element  Me¬ 
thode  zu  sehen,  die  in  diesem  Abschnitt  angesprochen  wird.  Ausdriicke  fur  die  fundamentale 
Losung  Uy*(x,y)  und  ty’(x,y)  werden  in  der  Literatur  gegeben.  Die  BEM  kann  stiickweise  zu 
einem  kontinuierlichen  Gebiet  (Domane)  aufgeweitet  werden,  wobei  jedes  Untergebiet 
(Subdomane)  unterschiedliche  aber  homogene  Materialeigenschaften  haben  kann.  Durch 
Anwendung  von  Gleichung  2-22  auf  jedes  Untergebiet  und  Ausformulierung  der  Interface- 
bedingungen,  um  die  zusatzlichen  Unbekannten  zu  eliminieren,  wird  ein  komplettes  System 
von  Gleichungen  aufgestellt,  um  die  Losungen  des  Randwertproblems  zu  erhalten.  Fur  ein 
perfektes  Interface  ergibt  sich  aufgrund  der  Kontinuitat  der  Verschiebungen  und  Gleichge- 
wichtsbedingungen  am  Interface: 

u‘  =uu  (2-24) 

t‘  =t"  (2-25) 

Hierbei  beziehen  sich  die  Indizes  auf  das  individuelle  Untergebiet  (Subdomane). 
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Reibungseffekte  am  Interface  konnen  ebenso  eingeschlossen  werden,  wie  spezielle  Abhlin- 
gigkeiten  zwischen  den  Interfaceverschiebungen  und  der  Haftung.  Nun  wird  auf  die  numeri- 
sche  Implementierung  eingegangen. 

2.7.4  Numerische  Implementierung 

Analytische  Losungen  der  Boundary  Integralgleichurg  2-22  sind  nicht  vollstandig  durch- 
flihrbar,  ausgenommen  fur  einfache  Gebiete  und  Randbedingungen.  Lediglich  erne  Na  he  rung 
fur  die  Gleichung  2-22  kann  hier  erfolgen.  In  der  ersten  Stufe  der  BEM  wird  der  Rand  in  ei- 
ne  Scrie  von  M-Elementen  eingeteilt,  in  denen  Geometrie,  Gesehwindigkeit  und  Haftung 
durch  stiickweise  Polynome  angenahert  werden.  In  der  hier  vorliegenden  Literaturstelle  wer¬ 
den  lineare  quadratische  Element  genutzt,  um  die  Geometrie  zu  erzeugen.  Uber  diesen  Ele- 
menten  wird  die  Haftung  und  Verschiebung  als  jeweils  konsiant  angenomrnen.  Die  Bounda¬ 
ry  Integralgleichung  ist  dann  wie  folgt  diskretisiert: 

cu  =  X  jjj5?  ul(x,y)dS(y)y<ll  -ft  tl(x,y)dS(y)  jrj  |  (2-26) 

Die  Anwendung  der  Gleichung  auf  eine  Serie  von  bestimmten  Punkten  auf  S  iuhrt  zu  einer 
Serie  von  linearen  algebraischen  Gleichungen,  die  fur  die  iibrigen  Randwerte  gelost  werden 
konnen. 


2.7.5  Ermittlung  der  Dehnungsenergiefreisetzungsrate 

Fur  die  Anwendung  der  BEM  stehen  grundsatzlich  zwei  Methoden  zur  Ermittlung  der  Ener- 
giefreisetzungsrate  zur  Verfugung: 


1)  die  sogenannte  differentielle  Energiefreisetzungsmethode 

2)  die  Ermittlung  der  Energiefreisetzungsrate  nach  Kaninen  [30] 


Die  differentielle  Energiefreisetzungsmethode  beruht  auf  der  direkten  Berechnung  der 
Dehnungsenergiefreisetzungsrate: 


(2-27) 
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G  =  Ul  —  +  O(Ad)  (2-28) 

A A 

Hierbei  ist  U  die  Dehnungsenergie  des  Systems  und  A  ist  die  Oberflache  des  Risses.  Die 
zweite  Gleichung  ist  eine  numerische  Implementierung  einer  Definition,  wie  sie  aus  einem 
System  vor  und  nach  dem  Rissfortschritt  extrahiert  wird.  Der  Rissfortschritt  betragt  dabei  Aa 
und  erzeugt  dabei  eine  Anderung  der  Oberflache  mit  AA.  In  der  Praxis  wird  die  inkrementel- 
le  Risslange  solange  reduziert,  bis  G  konstant  ist. 

Bei  der  Methode  nach  Kaninen  [30]  wird  die  Energiefreisetzungsrate  als  Arbeitsterm  for- 
muliert,  der  aufgewendet  werden  muss,  um  den  Riss  wieder  zu  schlieBen: 

1  Aa 

G  =  lim - f  t(l)Au(2)da  (2-29) 

2A a  l 

wobei  t(1)  die  Kraft  entlang  der  Risserweiterung  der  Lange  Aa  ist,  nachdem  der  Riss  ge- 
schlossen  wurde,  und  Au(2)  die  Rissoffnungsverschiebung  des  Risses  ist,  bevor  dieser  ge- 
schlossen  wurde.  Fiir  die  Anwendung  der  Boundary  Element  Methode  wird  die  Methode  des 
RissschlieBens  angewendet,  da  diese  monoton  konvergiert  und  numerisch  als  recht  stabil  er- 
wiesen  hat. 

2.7.6  Durchfiihrung  der  Boundary  Element  Methode 

Bei  der  Diskretisierung  wurden  die  acht  nachsten  Elemente  die  radial  von  der  Ursprungsriss- 
spitze  positioniert  sind,  mit  einer  konstanten  Lange  angenommen. 


Abb.2-10:  Diskretisierung  des  Rissfortschrittes  mit  Hilfe  der  BEM 
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Diese  Elemente  werden  genutzt,  um  die  Dehnungsenergiefreisetzungsrate  mit  Hilfe  der  Glei- 
chung  2-29  zu  bestimmen.  Die  sich  anschlieBenden  Elemente  haben  gemaB  Abb.  2-10  unter- 
schiedliche  GroBen.  Die  Abb.  2-11  vergleicht  die  Werte  von  He  und  Hutchinson  [5]  mit  de- 
nen  aus  der  BEM  gemaB  Tullock  und  Petrovic  [9]  errechneten. 


Abb. 2-1 1 :  Vergleich  des  Energiefreisetzungsratenverhaltnisses  fur  einfache  und  doppelte 
Rissabzweigung  am  Interface  gemaB  He  und  Hutchinson  [5]  und  Tullock  et  al.[9] 

2.7.7  Zusammenfassung  der  Rechnungen  von  Tullock.  Reimanis,  Graham  und  Petro¬ 
vic 

Anhand  der  Analyse  zeigt  sich,  dass  die  BEM-Methode  die  analytischen  Berechnungen  von 
He  und  Hutchinson  bestatigt.  GemaB  Tullock  und  Petrovic  [9]  ist  die  BEM  eine  gute  Metho- 
de,  um  Energiefreisetzungsraten  zu  ermitteln.  Im  Verlauf  dieser  Arbeit  wird  als  numerisches 
Losungsverfahren  die  Finite  Element  Methode  verwendet.  Hierzu  muss  allerdings  noch  an¬ 
hand  von  bekannten  Losungen  in  Analogic  zu  den  o.a.  Rechnungen  gepriift  werden,  ob  diese 
sich  ebenfalls  fur  die  Ermittlung  der  Energiefreisetzungsrate  eignet. 
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2.8  Rissabzweigung  an  einem  Interface  unter  Beriicksichtigung  des  Einflusses  von 
Restspannungen  nach  He,  Evans  und  Hutchinson  [101 

2.8.1  Einleitung 


In  den  bisher  vorgestellten  analytischen  Ansatzen  wurden  Verbundwerkstoffe  betrachtet,  de- 
ren  Faser  und  Matrix  isotrope  Eigenschaften  haben.  Hierbei  wurde  ausschliefilich  die  Ab- 
hangigkeit  von  Gd/Gp  vom  Dundursparameter  [14]  a  betrachtet,  ohne  jedoch  jegliche  Rest¬ 
spannungen  zu  beriicksichtigen.  Im  folgenden  wird  nun  ein  analytischer  Ansatz  diskutiert, 
der  auch  Restspannungen,  die  sich  aus  der  axialen  Komponente  Ot  und  der  Querkomponente 
o„  gemaB  Abb.  2-12  zusammensetzen,  beriicksichtigt. 


Abb. 2-1 2:  Spannungen  a,  und  an  in  der  Rissebene 


Hierzu  werden  zwei  zusatzliche  dimensionslose  Parameter  eingefuhrt: 

nn 


k, 


n,  = 


G.Vp 


kI 


(2-30) 

(2-31) 
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Hierbei  ist  a  die  Lange  des  Rissfortschrittes,  die  entweder  entlang  am  Interface  mit  ad  oder  in 
der  Faser  mit  ap  bezeichnet  werden.  X  ist  der  Exponent  der  Spannungssingularitat  fur  den 
Hauptriss.  Die  Variable  ki  ist  ein  Faktor,  der  proportional  zu  dem  entsprechenden  Span- 
nungsfeld  im  ungestorten  Bereich  ist.  Die  Auswirkungen  dieser  Parameter  auf  das  weitere 
Verhalten  von  Matrixrissen  wird  nun  diskutiert.  Auch  hier  wird  wieder  betrachtet,  ob  ein 
Riss  zum  Faserbruch  fuhrt  oder  lediglich  ins  Interface  abzweigt.  Neben  den  Dundursparame- 
tem  a  und  (3  wird  auch  hier  wieder  8  in  Abhangigkeit  von  (3  gemaB  Gleichung  2-16  einge- 
ftihrt. 


2.8.2  Spannungsintensitatsfaktoren  und  Energiefreisetzungsraten 

Der  Hauptriss  wird  als  halbunendlicher  Riss  angenommen,  der  senkrecht  zum  Interface  ver- 
lauft  und  dessen  Rissspitze  am  Interface  anstoCt.  Die  Spannungen  oxx  im  Material  unmittel- 
bar  vor  dem  Riss  (in  der  Faser)  werden  nach  Gleichung  2-11  charakterisiert.  Wie  von  He  und 
Hutchinson  [5]  bereits  diskutiert,  wird  angenommen,  dass  der  Rissfortschritt  des  Rissastes 
(ap  oder  ad)  klein  im  Vergleich  zur  Lange  des  Hauptrisses  (Matrixrisses)  ist.  Im  Fall  der 
Risspenetration  ist  der  Spannungszustand  an  der  fortschreitenden  Spitze  die  Folge  einer 
Mode  I-Belastung.  Der  Spannungsintensitatsfaktor  K  hangt  dann  von  k,,  ap  und  den  Rest- 
spannungen  parallel  am  Interface  ot  wie  folgt  ab: 

1  I 

Kx  =  c(a)k,ap2~x  +h{a)atap2  (2-32) 


Hierbei  sind  c  und  h  dimensionslose  Funktionen  von  den  Dundursparametem  a  und  (3.  Um 
die  Anzahl  der  Parameter  in  den  derzeitigen  Ergebnissen  zu  reduzieren,  wird  ausschlieBlich 
der  Einfluss  von  a  fur  (3=0  untersucht. 

Die  Energiefreisetzungsrate  des  Penetrationsrisses  ergibt  sich  zu: 


<?,=■ 


1-U. 

2jtii 


9  l-V. 
K ,2  =  1 


2  fii 


c2k)a\  11  +2 chk,G,a{p  x  +h2ofap] 


(2-33) 
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Fur  abzweigende  Risse  sind  die  Spannungen  am  Interface  direkt  vor  der  rechten  Rissspitze 
nach  Rice  [16]  wie  folgt  charakterisiert: 

a„(xf>)+i<yv(xfi)  =  (K,  +iK1'pnr)iri’  (2-34) 

wobei  r=x-a.  Mit  p=0  konnen  Ki  und  K2  als  konventionelle  Mode  I  und  Mode  II  Spannungs- 
intensitatsfaktoren  behandelt  werden,  sodass: 


ayy=Kl(27tf)-2  (2-35) 

axy  =  K2(2nr)~2 

Durch  Betrachtungen  analog  zu  He  und  Hutchinson  [5]  miissen  Ki  und  K2  von  kj,  a  und  den 
in  der  Ebene  senkrecht  zum  Interface  liegenden  Restspannungen  Gn  gemaB 


i-A  i 

K]  +iK2  =k,aj  d(cc)+Gnajg(a) 


(2-36) 


abhangen.  Die  Parameter  d,  g  und  h  sind  dimensionslose  komplexe  Funktionen  von  a,  auf 
deren  Herleitung  nicht  weiter  eingegangen  wird.  Die  Spannungskomponente  Gt  besitzt  kei- 
nen  Einfluss  auf  Ki  und  K2,  solange  sie  parallel  zum  Rissfortschritt  wirkt.  Die  Energiefrei- 
setzungsrate  betragt  dann  an  der  Spitze  des  abzweigenden  Risses: 

Gd  +kIena"-l)(dg  +  dg)  +  \g\2<J2nad}  (2-37) 


E,  ist  gegeben  durch: 


Et 


2 


M2) 


l-/32_ 

(\-a)E{ 


(2-38) 
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Das  Verhaltnis  Gd/Gp  ergibt  sich  zu: 


d\2  +  Tin 

dg  +  dg)  +  rj2  g2 

(  V 

ac 

[c2  +  2ritch  +  rfh2  ](l  -  a) 

KUP  > 

(2-39) 


Das  Verhaltnis  der  Energiefreisetzungsraten  hangt  damit  nur  von  den  dimensionslosen  Pa- 
rametem  a,  T|t,  T|n,  und  ad/ap  ab. 

Die  Losungen  der  analytischen  Gleichungen  bzw.  die  Bestimmung  der  Unbekannten  soli 
hier  nicht  beschrieben  werden.  Sie  erfolgte  anhand  einer  Finiten  Element  Methode  unter 
Aufbringen  der  entsprechenden  Spannungen.  Das  Verhaltnis  Gd/Gp  fur  ap=ad  wird  im  folgen- 
den  als  Funktion  von  a,  r|t  und  r|n,  wobei  Gd  ftir  den  einfach  abzweigenden  Riss  ermittelt 
wurde,  dargestellt.  Die  Abb.  2-13  bis  2-15  zeigen  die  Abhangigkeiten  fur  die  Sonderfalle 
T|t=0,  0,1  und  -0,1  als  Funktion  von  r|n  und  a. 


Abb. 2-1 3:  Energiefreisetzungsratenverhaltnis  als  Funktion  von  a  in  Abhangigkeit  von  r|n  fur 

Tit  =“0,1 
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Auffallend  ist,  das  bei  negativem  G„,  d.h.  das  Interface  unterliegt  einer  Druckbeanspruchung, 
die  Energiefreisetzungsverhaltnisse  kleiner  werden,  was  auf  eine  Tendenz  zur  Penetration 
des  Matrixrisses  in  die  Faser  hindeutet.  Positive  Werte  von  on  und  damit  positives  T|n  (Zug- 
beanspruchung  im  Interface  )  erzeugen  ein  groBes  Gd/Gp-Verhaltnis,  was  zu  Rissabzweigen 
fiihrt. 


2.8.3  Zusammenfassung  der  analytischen  Losung  fur  Restspannungsbehaftete  Faser- 
Matrix-Kombinationen 


Anhand  der  Restspannungen,  die  sich  entlang  des  Defektes  am  Interface  ad  und  in  der  Faser 
ap  ermitteln,  konnen  die  Beiwerte  rit  und  r|n  bestimmt  werden.  Aufgrund  dieser  Beiwerte 
wird  von  He,  Evans  und  Hutchinson  [10]  eine  Analyse  durchge fiihrt,  die  einen  signifikanten 
Einfluss  der  Restspannungen  zeigt,  wenn  diese  in  dem  unmittelbaren  Bereich  des  fortschrei- 
tenden  Rissastes  ap  oder  ad  auftreten.  Fiir  die  verwendete  Materialkombination  (Dundurspa- 
rameter  a)  wird  das  Energiefreisetzungsratenverhaltnis  in  Abhangigkeit  von  den  entspre- 
chenden  Beiwert  ermittelt.  Durch  die  Kenntnis  der  Bruchenergie  der  Faser  bzw.  des  Inter¬ 
faces  kann  ermittelt  werden,  ob  der  Riss  im  Penetrationsgebiet  oder  im  Debondinggebiet 
liegt. 

Der  Einfluss  der  Restspannungen  auf  Gd/Gp  des  kohlenstofffaserverstarkten  Kohlenstoffs 
wird  in  einer  FE-Analyse  in  dieser  Arbeit  analysiert. 
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2.9  Energiekriterium  fiir  eine  Rissabzweigung  an  einem  Interface  zwischen  zwei 
orthotropen  Medien  gemafi  Gupta  und  Martinez  [111 

Wie  bereits  erwahnt  sind  die  bisherigen  analytischen  Ansatze  nur  fiir  Faserverbundwerkstof- 
fe  anwendbar,  die  aus  isotropen  Fasem  und  isotroper  Matrix  bestehen.  Ihre  Anwendung  auf 
den  dieser  Arbeit  zugrundeliegenden  kohlenstofffaserverstarkten  Kohlenstoff  ist  somit  nur 
beschrankt  moglich.  Die  Autoren  Gupta  und  Martinez  [11,19]  haben  jedoch  einen  theoreti- 
schen  Ansatz  aufgestellt,  mit  dem  Rissabzweigen  und  Penetration  in  anisotropen  Werkstof- 
fen  beschrieben  werden  kann.  Im  folgenden  wird  dieser  Ansatz  im  wesentlichen  beschrieben. 

2.9.1  Allgemeines 

Auch  fur  Faserverbundwerkstoffe  mit  orthotropen  Materialkomponenten  wird  wieder  von 
zwei  grundsatzlichen  Moglichkeiten  des  Rissfortschritts  eines  Matrixrisses  ausgegangen.  Der 
Riss  kann  entweder  in  die  Faser  fortschreiten  oder  er  kann  entlang  des  Interfaces  abzweigen. 
Die  Art  und  Weise  des  Rissfortschrittes  kann  mit  Hilfe  von  zwei  Kriterien  vorausgesagt 
werden: 

•  Auf  der  Grundlage  der  Spannungen  an  der  lokalen  Rissspitze 

•  Aufgrund  des  Unterschiedes  in  der  Rissenergie  entlang  eines  moglichen  Risspfades 
GemaB  des  Spannungskriteriums  wird  der  Riss  entlang  des  Interfaces  abzweigen,  wenn  fol- 
gende  Beziehung  gilt: 


a;  <  ojsn 

o'f  crxt(90°) 


(2-40) 


wobei  Gj*  und  Of*  die  Interface  bzw.  die  Faserfestigkeit  ist  und  Oyy(0°)  und  oxx(90°)  die 
Spannungen  im  Interface  und  in  der  Faser  darstellen,  welche  durch  die  Losung  von  Integral- 
gleichungen  nach  Gupta  [19]  berechnet  werden  konnen.  Die  Analyse  setzt  Anisotropie  in  Fa¬ 
ser  und  Matrix  voraus  und  ist  deshalb  fur  kohlenstofffaserverstarkten  Kohlenstoff  (C/C)  von 
Interesse.  Das  oben  beschriebene  Kriterium  ist  fiir  die  Materialentwicklung  von  Bedeutung  , 
da  es  festlegt,  welche  Interface-  und  Fasereigenschaften  fur  ein  positives  Schadigungsverhal- 
ten  (Rissabzweigen)  zu  realisieren  sind. 

Neben  diesem  Spannungskriterium  haben  Gupta  und  Martinez  ein  Energiekriterium  fiir  die 
Rissentwicklung  in  Verbundwerkstoffen  mit  orthotropen  Komponenten  entwickelt. 
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Dieses  fordert  wiederum,  dass  die  Energiefreisetzungsrate  fiir  eine  Risserweiterung  entlang 
des  Interfaces  Gd  im  Verhaltnis  zur  Penetration  Gp  grofler  ist  als  die  Interfacefestigkeit  Gc‘ 
im  Verhaltnis  zur  Festigkeit  des  nichtgerissenen  Materials  Gcf  (Gleichung  2-6).  Im  folgenden 
wird  nun  die  Ermittlung  der  Energiefreisetzungsraten  durchgefuhrt,  indem  die  bekannten  Lo- 
sungen  fur  isotrope  Medien  nach  He  und  Hutchinson  [5]  auf  orthotrope  Medien  erweitert 
werden.  Wiederum  ausgehend  von  den  drei  unterschiedlichen  Arten  des  Rissfortschrittes 
(Penetration,  einfaches  Abzweigen  und  doppeltes  Abzweigen  am  Interface  )  erfolgt  die  Er¬ 
mittlung  der  Energiefreisetzungsrate,  wobei  fur  die  beiden  letzten  Falle  derjenige  ausschlag- 
gebend  ist,  der  den  groBten  Wert  von  Gd  erzeugt. 

2.9.2  Riss  senkrecht  zu  einem  Interface 

2.9.2.1  Problemstellung  und  Losungsansatz 

Eine  symmetrische  Last  greift  an  einem  Verbund  aus  zwei  orthotropen  Medien  an,  wobei  ein 
halbunendlicher  Riss  senkrecht  auf  das  Interface  zwischen  beide  stoBt.  Beide  Medien  werden 
als  homogen  und  linear  elastisch  angenommen.  Das  ungeschadigte  Material  (Faser)  erhalt 
den  Index  I,  das  rissbehaftete  Material  (Matrix)  bekommt  den  Index  II.  Fiir  anisotrope  Mate- 
rialien  gilt  das  Hooksche  Gesetz: 

6 

£.  =  ^SyUj  fur  i  =  1  bis  6  (2-41) 

j= i 

Hierbei  wird  die  folgende  Schreibweise  verwendet: 

®  22  ^33  *^23  ^31  &  12}  (2-42) 

{e/}  =  {^ll  £22  £33  ^23  £31  £  12  }  (2-43) 

Die  Elemente  sy  sind  Elemente  der  Nachgiebigkeitsmatrix,  welche  21  Elemente  fur  den  all- 
gemeinen  Fall  besitzt.  Fur  den  orthotropen  Fall  reduziert  sich  die  Spannungsdehnungsbezie- 
hung  auf: 

£,=  YabiJ°j  Sir  i=  1,2,  6  (2-44) 

7=1, 2, 6 
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mit 


bu  = 


Sjj . . .  ebener  Spannungszust. 

S 

s - — . ..ebener  Dehnungszust. 

*33 


(2-45) 


mit  nur  vier  unabhangigen  elastischen  Nachgiebigkeiten:  bn,bi2=b2i>b22  und  b^e-  GemaB 
Stroh  [20]  kann  das  elastische  Spannungsfeld,  im  folgenden  auch  Potential  genannt,  fur  die¬ 
ses  ebene,  elastische  Problem  anhand  von  zwei  komplexen  Funktionen  dargestellt  werden: 

(f>,  (z, )  mit  i=l,2  (2-46) 

Jede  dieser  Funktionen  ist  holomorph  in  ihrem  Argument  z^x+Kjy.  Fiir  ein  orthotropes  Ma¬ 
terial  gibt  es  zwei  komplexe  Zahlen  Kj  mit  positivem  Imaginarteil,  welche  die  Wurzeln  der 
charakteristischen  Gleichung  2-47  darstellen. 

AjfC4  +  2p4XxK2  +1  =  0  (2-47) 

Hierbei  sind  A,]  und  p  dimensionslose  elastische  Konstanten,  die  gemaB  Suo  [21]  wie  folgt 
definiert  sind: 


2bn  +  b66 
2yjbnb22 


(2-48) 


Diese  Konstanten  zeigen  den  Grad  der  Anisotropie  auf.  Fiir  isotrope  Materialien  gilt  A,  =p 
=  1.  Die  vorliegende  Berechnungen  sind  auf  l<p<oo  begrenzt,  da  die  meisten  Materialien  sich 
in  diesem  Bereich  bewegen.  Fiir  diesen  Fall  konnen  die  Wurzeln  nach  Suo  [21]  wie  folgt 
ermittelt  werden: 


(2-49) 

k2  -  iX^ /4(n  —  m)  =  it2 

(2-50) 
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Hierbei  sind  n  und  m  wie  folgt  festgelegt: 


(2-51) 


m  = 


(2-52) 


Mit  diesen  Losungen  fur  die  charakteristische  Gleichung  kann  das  elastische  Feld  nun  durch 
die  zwei  holomorphen  Funktionen 

0;.(z(.)  miti=l,2  (2-53) 


ausgedriickt  werden.  Nach  Stroh  [20]  konnen  die  Verschiebung  u,,  die  Spannungen  Gy  und 
die  resultierenden  Krafte  f;  wie  folgt  angenahert  werden: 


u .  = 


<*2,  = 


2R 


j= 1 


G\i 


(2-54) 


(2-55) 

(2-56) 

(2-57) 


Hierbei  stellt  <J)’  die  Ableitung  nach  dem  entsprechenden  Argument  dar.  Ay  und  Ly  sind  die 
Elemente  der  Matrizen  A  und  L,  die  von  den  elastischen  Konstanten  abhangen.  Unter  der 
Annahme  des  ebenen  Spannungszustandes  gilt: 


A  = 


snKi  +s\ 


Sl\^2  5I2 


■^12  "k 


’22 


K, 


S\2^2  4" 


’22 


(2-58) 
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(2-59) 


Stroh  [20]  zeigt,  dass  A  und  L  nicht  singular  fur  die  Wurzel  Kj  und  K2  sind.  Gupta  [19]  pos- 
tulierte,  dass  fur  zwei  orthotrope  Medien  die  Losung  von  p  und  den  zwei  von 
Dundursparametem  [14]  a  und  (3  abhangen,  die  fur  orthotrope  Werkstoffe  wie  folgt  defmiert 
sind: 


a  = 


{■\JSUS22  )/;  (V5 11522  )y 

[■\lSl  lS 22  )/7  +  (yJsuS22  ^ 

{■\ISUS22  "*■  ^12  )/7  ~  {tylSHS22  +  ‘^12  ), 


*Jhuh22 


(2-60) 


(2-61) 


i/n  —  (2«A,  Js\is22  )7  +  {in A,  (2-62) 

H21={lnXl  yjsus22  );  +  {in A,  ^sns22  \i  (2-63) 


Die  Parameter  a  und  [)  driicken,  wie  bereits  von  He  und  Hutchinson  [5]  fur  isotrope  Materia- 
lien  beschrieben,  den  Steifigkeitsunterschied  („mismatch“)  zwischen  den  zwei  Materialien 
aus.  Sie  nehmen  Werte  zwischen  -1  und  1  an  und  sind  im  homogenen  Fall  0.  Bei  bekanntem 
Potential  0(.(z,.)  mit  i=l,2  konnen  anhand  der  Gleichungen  2-56  und  2-57  die  Spannungen, 
aus  den  Gleichungen  2-55  die  resultierenden  Krafte  und  gemaB  Gleichung  2-54  die  Ver- 
schiebungen  berechnet  werden.  Martinez  und  Gupta  [11]  verwenden  zur  Herleitung  des  Po¬ 
tentials  die  Darstellung  eines  Risses  gemaB  Abb.  2-16.  Die  Risse  werden  durch  Versetzun- 
gen  ersetzt.  Diese  Versetzungen  werden  in  jedem  Medium  nahe  des  Interfaces  eingefugt,  um 
einen  Penetrationsriss  zu  simulieren.  Die  Simulation  des  Interfacerisses  erfolgt  durch  das 
Einbringen  von  zwei  Versetzungen  im  Interface. 
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Abb.2-16:  Einfuhrung  von  Versetzungen  zur  Risssimulation 

a)  einfache  Versetzung  in  der  Matrix 

b)  einfache  Versetzung  in  der  Faser 

c)  Versetzungen  zur  Simulation  von  Interfacerissen 

Anhand  dieser  Versetzungen  ist  es  moglich,  das  Potential  und  damit  das  ebene  Spannungs- 
feld  zu  bestimmen.  Fur  die  ausfuhrliche  Herleitung  sei  hierbei  jedoch  auf  die  Literatur  [19] 
verwiesen,  da  sich  die  notwendige  Losung  der  Integralgleichung  als  recht  umfangreich  dar- 
stellt.  Im  folgenden  wird  daher  nur  auf  die  Darstellung  der  Ergebnisse  eingegangen,  die  hin- 
sichtlich  des  kohlenstofffaserverstarkten  Kohlenstoffs  von  Bedeutung  sind. 

2.9.3  Energiekriterium  fur  Rissabzweigung  Oder  Risseindringen 

Im  folgenden  wird  ein  Ausdruck  Gd/Gp  auf  der  Grundlage  von  dimensionalen  Betrachtungen 
ermittelt.  Fiir  alle  Arten,  des  hier  behandelten  Rissfortschrittes  befindet  sich  der  halbunendli- 
che  Referenzriss  mit  a=0  senkrecht  zum  Interface,  sodass  er  mit  seiner  Spitze  das  Interface 
beriihrt. 

2.9.3. 1  Ermittlung  der  Energiefreisetzungsrate  Gn  fiir  Penetration 

Fur  eine  symmetrische  Mode-I-Belastung  in  Bezug  auf  die  Rissebene  kann  die  an  Spannung 
bei  0=90°  vor  der  Rissspitze  in  der  Faser  (Material  I)  durch  folgenden  Ausdruck  beschrieben 
werden: 


0, 7)^(2/^ 


(2-64) 
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wobei  kj  der  Amplitudenfaktor  und  y  die  Spannungssingularitat  unter  der  Annahme  von  rea- 
len  Werten  zwischen  -l<y<0  ist.  Gemafi  der  Auswertung  der  Integralgleichungen  [1 1]  fur  die 
Losung  des  Potentials  kann  y  wie  folgt  bestimmt  werden: 

t2Dn -(,£>„  +r,D12(V'  -(2)cos(^r)  =0  (2-65) 

t2  *1 

wobei  tj=n+m,t2=n-m  ist.  Die  Gleichung  2-65  kann  fur  die  Matrix  (Material  II)  analog  be- 
schrieben  werden.  Die  Konstanten  Djj  2sind  Elemente  der  Matrix,  die  von  a,  P,  X  Und  P  ab- 
hangen  und  die  die  dazugehorige  Losung  einer  Randverschiebung  nahe  dem  Interface  in  der 
Matrix  (Material  II)  beschreiben.  Sie  werden  ebenfalls  aus  der  Losung  der  Integralgleichun¬ 
gen  gewonnen. 

Fur  den  halbunendlichen  Riss  wird  entweder  angenommen,  dass  er  eindringt  oder  einfach 
bzw.  zweifach  abzweigt. 

Im  Fall  des  Eindringens  wird  der  Spannungszustand  an  der  fortschreitenden  Rissspitze  als 
einfacher  Modus  I  beschrieben. 

Der  Spannungsintensitatsfaktor  K,  wird  in  Abhangigkeit  vom  Amplitudenfaktor  k,  und  der 

Risslange  a  aus  dimensional en  Grunden  wie  folgt  beschrieben: 

i 

Kj  =  ckia2+r  (2-66) 

Hierbei  ist  c  eine  dimensionslose  Konstante,  die  von  den  elastischen  Parametem  a,  (3,  X  und 
p  abhangt.  Die  Energiefreisetzungsrate  fur  den  eindringenden  Riss  wird  gemaB  Sih  [24]  wie 
folgt  angegeben: 

zi  _I 

Gp  =  (sun){  X*  K 1  -  {snn)XAc2k[ax+2y  (2-67) 


2  Die  Gleichungen  fur  Dy  werden  von  Gupta  [22,  23]  mit  einer  anderen  Bezeichnung,  namlich  Cij* 
gefuhrt. 
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2.93.2  Ermittlung  der  Energiefreisetzungsrate  fur  Deflection 

Die  Spannung  am  Interface  an  der  Rissspitze  an  jedem  der  abzweigenden  Risse  (siehe  hierzu 
auch  Abb.2-7b)  wird  nach  Suo  [21,23]  wie  folgt  bestimmt: 


(J°22(0,x)  +  i 


Ht 


\h22 


II  0” 

kJ  n 


(0,x)  = 


Kr,E 

-JlTtr 


(2-68) 


Hierbei  ist  K  der  komplexe  Spannungsintensitatsfaktor  fur  den  Interfaceriss.  Hn  und  H22 
sind  die  Koeffizienten,  welche  durch  die  Gleichung  2-62  und  2-63  mit  r=x-a,  i=V-l  ermittelt 
werden.  Die  Bestimmung  von  E  erfolgt  aus  Gleichung  2-16.  Aus  einer  dimensionalen  Be- 
trachtung  heraus  ergibt  sich  schlieBlich  folgender  Ausdruck: 

i 

K  =  kxa?\dai£  +  ea-ie )  (2-69) 


Hierbei  sind  d  und  e  dimensionslose  komplexe  Konstanten. 

Die  zugehorige  Energiefreisetzungsrate  fur  den  abzweigenden  Riss  ist: 


G.  = 


H. 


22 


W2  =4,fhi +N2  +2Re(<te)} 


4  cosh  K£ 


4  cosh  K£ 


(2-70) 


2.9.3.3  Ermittlung  des  Quotienten  der  Energiefreisetzungsraten 

Der  Quotient  der  Energiefreisetzungsraten  Gj/Gp  ergibt  sich  aus  Gleichung  2-67  und  Glei¬ 
chung  2-70  und  ist  unabhangig  von  a  und  kg 


5i 


h22  V 

4 cosh2  Tt£  (5-, ,/z) , c2 


[\d\2  +|ef  +2R  e(fife)] 


(2-71) 


Es  kann  ein  Wert  der  Energiefreisetzungsrate  bestimmt  werden,  der  wiederum  die  Wahl  des 
Rissverlaufes  bestimmt.  Fur  eine  Auswahl  von  Werkstoffen  wurde  das  entsprechende  Ver- 
haltnis  gemaB  Gleichung  2-71  ermittelt  und  in  Abb.  2-17  in  Abhangigkeit  vom  Dundurspa- 
ramter  a  [14]  eingetragen. 
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C/C  (XFase  =0,067,  pFase  =  1,1 62)  a  =  0,248  Gd/Gp  =  0,118 

SiC  Ni;.\l  a  - 0.385  Gd/Gp  -  0,3859 

E-Glass/Epoxy  a  =  0,9091  Gj/Gp^  1,9315 


C/C  (AFaser=0,067,  pFase  =  1,162)  oc  =  0,248  Gd/Gp  =  0,118 

SiC  Ni;.\l  a  - 0.385  Gd/Gp  -  0.3859 

E-Glass/Epoxy  a  =  0,9091  G^Gp^  1,9315 


Abb.2-17:  Gd/Gp  fur  orthotrope  Materialien  [1 1,  19]  in  Abhangigkeit 
vom  Dundursparameter  a  und  Orthotropiefaktor  Xt 

Abb.  2-17  zeigt,  dass  fur  bestimmte  Materialkombinationen  das  Verhaltnis  Gd/Gp  stark  vari- 
iert.  Liegen  die  Bruchenergieverhaltnisse  Gc‘/Gcf  unterhalb  der  fair  die  entsprechende  Kombi- 
nation  ermittelten  Datenpunkte,  so  erfolgt  Rissabzweigung.  GroBe  Gd/Gp-Werte  deuten  auf 
eine  Tendenz  zur  Rissabzweigung  hin,  wohingegen  kleine  Verhaltnisse  von  Gd/Gp  eher  ein 
Eindringen  des  Matrixrisses  in  die  Faser  begiinstigen. 


39 


Rissausbreitung  in  kohlenstofffaserverstarktem  Kohlenstoff  (C/C) 
Thorsten  Weber 


DLR  -  Institut  fur  Werkstoff-Forschung 


2.9.4  Zusammenfassung  der  Theorie  von  Gupta  und  Martinez 

Anhand  des  Verhaltnisses  Gd/Gp  in  Abhangigkeit  vom  Dundursparameter  fur  Faserverbund- 
werkstoffe  aus  orthotropen  Komponenten  zeigt  sich,  dass  im  Gegensatz  zu  der  Betrach- 
tung  von  Faserverbunden  aus  isotropen  Materialien  nach  He  und  Hutchinson  [5]  die  Kur- 
ven  fur  das  Verhaltnis  Gd/Gp  durch  das  MaB  der  Anisotropie  beeinflusst  werden. 

Aus  dem  o.a.  Diagramm  kann  man  beispielsweise  fur  Faserverbundwerkstoffe  mit  orthotro- 
per  Faser  und  isotroper  Matrix  das  Verhaltnis  der  Energiefreisetzungsraten  Gd/Gp  bestim- 
men,  wie  es  beispielweise  bei  kohlenstofffaserverstarktem  Kohlenstoff  der  Fall  ist.  Das  blaue 
Kreuz  stellt  hier  den  Wert  fur  kohlenstofffaserverstarkten  Kohlenstoff  (C/C)  dar. 

Nach  Bestimmung  des  Bruchenergieverhaltnisses  Gc'/Gcf  konnen  Aussagen  iiber  das  Riss- 
fortschrittsverhalten  getroffen  werden. 

Allerdings  bleiben  hierbei  herstellungsbedingte  Eigenspannungen  im  Material  (z.B.  aufgrund 
von  unterschiedlichen  Warmeausdehnungskoeffizienten)  unberiicksichtigt. 
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2.10  Analvtische  Bestimmung  des  Faserdebondings  in  restspannungshehafteten 
sproden  Faserverbundmaterialien  von  Charalambides  [131 

2.10.1  Allgemeines 

Im  folgenden  wird  nun  ein  theoretischer  Ansatz  vorgestellt,  der  sich  mit  dem  Einfluss  von 
thermischen  Restspannungen  in  Faserverbundwerkstoffen  mit  sproder  Matrix  beschaftigt. 
Die  von  Charalambides  [13]  aufgestellte  Theorie  ist  als  Erweiterung  der  Untersuchung  von 
He  und  Hutchinson  zu  sehen. 

Wie  bereits  von  He  und  Hutchinson  berichtet,  tritt  bei  einer  unter  90°  zur  Matrixrissebene 
angreifenden  Last  Rissabzweigen  auf,  wenn  Gc‘/Gcf<  !4  [25,  26,  27]  ist.  Die  GroBen  Gc‘  und 
Gcf  stellen  die  Bruchenergie  des  Interfaces  und  der  Faser  dar.  Im  folgenden  wird  nun  dieses 
Kriterium  unter  der  Annahme  thermischer  Spannungen  verifiziert.  Faserverstarkte  Werkstof- 
fe  entwickeln  thermische  Restspannungen  wahrend  der  Herstellung.  Dies  ist  auf  die  Unter- 
schiede  in  den  thermischen  Warmeausdehnungskoeffizienten  in  Faser  und  Matrix  zuruckzu- 
fiihren. 

Die  im  folgenden  durchgefuhrte  Analyse  gilt  fur  Faserverbundwerkstoffe  mit  Hauptrissen, 
die  mit  komplexen  Spannungsintensitaten  beaufschlagt  werden,  wobei  thermische  und  me- 
chanischen  Lasten  sich  entweder  in  Phase  Oder  Gegenphase  befinden  konnen.  Dabei  werden 
besondere  Falle  der  Rissabzweigung,  wie  z.B.  Faserdebonding  in  unidirektionalem,  faserver- 
starkten  Verbundwerkstoffen  und  der  Delamination  von  Schichten  gezeigt. 

2.10.2  Mechanisches  Verhalten  beim  Debonding  (Matrixriss  senkrecht  zur  Faser) 

Wie  bereits  ausfuhrlich  diskutiert,  kann  das  konkurrierende  Verhalten  zwischen  Faserdebon¬ 
ding  und  Faserversagen  untersucht  werden,  indem  man  die  Energiefreisetzungsraten  GK  an 
der  Rissspitze  des  abzweigenden  Risses  (Risswinkel  betragt  90°)  mit  der  des  Hauptrisses  GT 
unter  der  aufgebrachten  mechanischen  und  thermischen  Belastung  vergleicht. 

Die  theoretischen  Uberlegungen  hierzu  werden  anhand  eines  zylindrischen 
Einzelfaserverbundes  mit  konzentrischem  Matrixriss  aufgestellt. 

Interessant  sind  hierbei  besonders  Werkstoffe,  bei  denen  mechanische  und  thermische  Lasten 
neben  einem  Modus  I  Spannungsintensitatsfaktor  fur  den  Matrixriss  noch  weitere  Moden  er- 
zeugen. 
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Solche  Bedingungen  treten  bei  der  Delamination  von  Schichten  bzw.  Debonding  von  geneig- 
ten  und  nebeneinander  liegenden  Laminaten  auf.  Fur  diese  Systeme  ist  eine  genaue  Analyse 
notwendig,  um  thermische  Effekte  auf  das  Rissfortschrittsverhalten  wie  Debonding  bzw.  De¬ 
lamination  zu  untersuchen.  Diese  Analyse  wird  im  folgenden  durchgefuhrt. 

2.10.3  Rissabzweigung  unter  Mixed-Mode-Bedingitng 

In  diesem  Abschnitt  wird  eine  Erweiterung  der  bisherigen  Theorien  beschrieben,  um  ther- 
mische  Effekte  auf  das  Debonding  oder  Rissabzweigen  an  der  Spitze  eines  Mixed  Mode 
Risses  zu  beschreiben.  Es  wird  eine  Situation  angenommen,  in  der  eine  aufgebrachte  mecha- 
nische  Spannung  t  und  eine  thermische  Spannung  0R,  mit  0R=EAaAT  als  Restspannung, 
vorhanden  ist.  Es  existieren  gemaB  Abb.  2-18  zwei  Anteile  zu  dem  komplexen  Spannungsin- 
tensitatsfaktor  KT  an  der  Spitze  des  Hauptrisses: 

1) K'  in  Bezug  auf  die  aufgebrachte  mechanische  Spannung  t 

2) Kr  in  Bezug  auf  die  thermische  Spannung  <JR 

Aus  existierenden  Losungen  und  Dimensionsanalysen,  kann  der  o.a.  Spannungsintensitats- 
faktor  in  der  komplexen  Form  geschrieben  werden: 

K‘  =K‘,  +iK‘,  =  ta^F'  (cosy/1  +  ismy/‘)  (2-72) 

\_ 

Kr  =  Kr  +  iKf,  =  oRa2FR(cosy/R  +  isiny/R)  (2-73) 

Hierbei  sind  Fl  und  FR  im  allgemeinen  Fall  eines  Zweikomponentensystems  Funktionen  der 
Geometrie  und  der  elastischen  Komponenten  der  beteiligten  Werkstoffe.  Durch  Superpositi¬ 
on  wird  der  kombinierte  komplexe  Gesamtspannungsintensitatsfaktor  KT  an  der  Spitze  des 
Hauptrisses  wie  folgt  erhalten: 


K1  =  K]  +iKR  =(l  +  —F-)Kt,+i(l  + 

K  r 


krR 

Kn 


(2-74) 


Fiir  den  Matrixriss  wird  die  totale  Energiefreisetzungsrate  in  Abhangigkeit  von  |  KT  |  gemaB 
Irwin  [29]  wie  folgt  bestimmt: 


Gt  = 


1  —  V' 


1  +  tanV) 


(2-75) 
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Hierbei  ist  E  der  Elastizitatsmodul  und  v  die  Querkontraktionszahl  des  Materials. 


Abb.2-18:  Spannungssingularitaten  am  Riss 

GemaB  Cotterrell  und  Rice  [28]  und  Suresh  und  Shih  besitzt  der  Spannungsintensitatsfaktor 
an  der  Spitze  des  abzweigenden  Risses  [51]  zwei  Komponenten,  Kik  und  K2\  und  kann  aus 
Kt,  dem  totalen  Spannungsintensitatsfaktor  an  der  Spitze  des  Hauptrisses,  erhalten  werden. 

K*  —  cnKj  +cnK'n  (2-76) 

K\  =c2]Kj  +  c22K'„  (2-77) 

Hierbei  sind  Cjj(i,j=l,2)  geometrische  Faktoren,  die  fur  ein  homogenes  System  in  Abhangig- 
keit  von  0  (Abb.2-18)  bestimmt  werden  und  mit  denen  von  Cotterell  und  Rice  [28]  fur  kleine 
Werte  von  0  identisch  sind. 

Fur  den  abzweigenden  Riss  ergibt  sich  die  Energiefreisetzungsrate  Gk  nach  Irwin  [29]  und 
unter  der  Einbeziehung  von  o.a.  Gleichungen  zu: 

G‘  =^[(Kl‘)!+(x2*)!]=E2l(A:j'')!c(c,)(e),^r)  (2-78) 

Jb  b 
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wobei 


c(cv  (0),y/T)=  c2  +  4  +  2  (cncl2  +  c2lc22  )tmy/T  +  (c22,  +  4  )tan2  y/T  (2-79) 

Der  Phasenwinkel  des  totalen  Spannungsintensitatsfaktors  an  der  Spitze  des  Matrixrisses 
tan\|/T=KnT/KiT  wird  aus  Gleichung  2-74  wie  folgt  erhalten: 


,  Fr  siny/*  oR 

t  1  - : - : - 

,  Kn  F  smy/  t  , 
tany/  =~t  = - * - —R — j-  tan  yr 

K]  t  |  F  cosy/*  oR 

F‘  cosy/'  t 


(2-80) 


2.10.4  Debondingbedingungen 


Der  Wettbewerb  zwischen  Rissabzweigung  (Faserdebonding)  und  Faserversagen  kann  da- 
durch  beschriebcn  werden,  dass  anhand  von  Gleichung  2-75  und  2-78  Bedingungen  fur  das 
Bruchenergieverhaltnis  Gc7Gcf  bestehend  aus  der  Bruchenergie  fur  das  Interface  und  der 
Bruchenergie  fur  die  Faser  bestimmt  werden 

Fur  den  Sonderfall,  dass  die  unter  90°  von  der  Ebene  des  Hauptrisses  abzweigende  Rissspit- 
ze  mit  Faserversagen  konkurriert,  wird  ein  Debonding  der  Faser  eintreten,  wenn 


Gc  ^  [^n  4~ c2i  2(cnc12  +  c21c22 )tany/  +(c]2+c22}t<m  y/  9=g/2 
GJC  l  +  tan2y rT 


(2-81) 


Fur  den  Fall,  dass  Rissabzweigung  unter  einem  anderen  Winkel  als  90°  zur  Ebene  des 
Hauptrisses  auftritt,  wird  die  giinstigere  Variante  dadurch  erhalten,  indem  man  die  Losungen 
fur  die  Energiefreisetzungsraten  der  zwei  Risse  (Faserdebonding  oder  Faserversagen)  ver- 
gleicht.  Faserdebonding  tritt  im  Gegensatz  zum  Faserversagen  eher  auf,  wenn: 


G‘c  ^  [cn  4" c21  +  2(cnc12  +  c2xc22 )tan y/  +  (cn  +c21)tan“  y/  9=0j 
G{  [ch  +c21  +2(cuc12  +  c21c22  )tany/  +  (cl2  4-  c22  )tan  y/  9=e 


(2-82) 


44 


Rissausbreitung  in  kohlenstofffaserverstarktem  Kohlenstoff  (C/C) 
Thorsten  Weber 


DLR  -  Institut  fur  Werkstoff-Forschung 


Hierbei  sind  0;  und  0f  die  Abzweigungswinkel  fur  Faserdebonding  und  Faserversagen.  Ge¬ 
maB  Gleichung  2-81  und  2-82  beeinflussen  die  thermischen  Lasten  die  Debondingbedingun- 
gen  iiber  den  Phasenwinkel  \|/T,  der  aus  Gleichung  2-80  ermittelt  wird.  Sobald  was 

einer  „in  Phasen“-Wirkung  der  thermisch/mechanischen  Spannungsintensitat  entspricht, 
fuhrt  Gleichung  2-80  dazu,  dass  \j/T=V|/.  Die  thermischen  Spannungen  wiirden  das  Verhaltnis 
Gc‘/Gcf  gemaB  Gleichung  2-81  und  2-82  nicht  beeinflussen. 

Wenn  allerdings  y*  dann  werden  die  o.a.  Bedingungen  gemaB  Gleichung  2-80  abhangig 
von  den  thermischen  Spannungen.  Es  ist  eine  zusatzliche  Analyse  erforderlich,  um  die  ther¬ 
mischen  Einfliisse  zu  spezifizieren.  Im  folgenden  werden  zwei  Falle  unterschieden. 

Im  ersten  Fall  wird  angenommen,  dass  thermische  und  mechanische  Spannungen  unabhan- 
gig  voneinander  sind,  was  mit  einer  „in  Phase“  der  Spannungsintensitaten  an  der  Spitze  des 
Hauptrisses  entspricht:  \j/T=\(/t==\(rR.  Die  thermischen  Effekte  fallen  aus  der  Debondingbedin- 
gung  heraus. 

Im  anderen  Fall  sind  thermische  und  mechanische  Lasten  nicht  „in  Phase“:  V)/*  Debon¬ 
ding  wird  abhangig  von  den  thermischen  Lasten. 

2.10.5  Debonding  unabhangig  von  thermischen  Lasten 

Matrixriss  senkrecht  zur  Faser 


Wenn  mechanische  und  thermische  Spannungen  unabhangig  voneinander  sind,  was  einer  „in 
Phase“  der  Spannungsintensitaten  an  der  Hauptrissspitze  entspricht,  fuhrt  Gleichung 

2-80  zu  \|rT=\pt=V)/R.  Als  Resultat  werden  die  Debondingbedingungen  iiber  Gleichung  2-81 
unabhangig  von  den  thermischen  Spannungen  und  konnen  aus  den  mechanischen  Spannun¬ 
gen  alleine  ermittelt  werden. 

Dieser  Fall  stellt  den  Faserbondingprozess  in  einem  Einzelfaserverbund  dar.  Der  Matrixriss 
wird  von  einem  Mode  I  Spannungsintensitatsfaktor  bestimmt,  der  entweder  von  einer  me¬ 
chanischen  oder  thermischen  Last  hervorgerufen  wird,  sodass  \)/T=\J/t=\|/R=0. 
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4  Punkt  Biegeprobe  aus  zwei  Komponenten 

Unter  der  Annahme  des  ebenen  Dehnungszustand  in  einer  vorgekerbten  Zweikomponen- 
tenprobe  existiert  eine  Mode  I  Bedingungen  fur  die  thermische  und  aufgebrachte  mechani- 
sche  Spannung 


Abb.  2-19:  4-Punktbiegeprobe 

In  diesem  Fall  wird  unter  Beriicksichtigung  von  Gleichung  2-81  Delamination  der  Ober- 
schicht  gemaB  Abb.  2-19  auftreten,  wenn: 

fMc"+4,L,2  P-83) 

2.10.6  Debonding  abhangig  von  thermischen  Lasten 

Im  Allgemeinen  ist  der  Spannungsintensitatsfaktor  an  der  Spitze  des  Hauptrisses  in  Bezug 
auf  die  thermischen  Spannungen  „in  Gegenphase“  zu  den  mechanischen  Spannungen: 

Die  Debondingbedingungen  gemaB  Gleichung  2-81  und  2-82  sind  tatsachlich  abhan¬ 
gig  von  den  thermischen  Effekten  gemaB  Gleichung  2-80. 
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Um  den  Einfluss  thermischer  Effekte  auf  die  obigen  Bedingungen  zu  demonstrieren,  wird 
die  4-Punkt  Biegeprobe  gemaB  Abb. 2-1 9  herangezogen.  Hierbei  handelt  es  sich  um  einen 
homogenen  Faserverbundwerkstoff  mit  Ej=E2  und  Vi=v2,  wobei  die  Indizes  1  und  2  fur  das 
obere  und  untere  Material  stehen.  Die  beiden  thermisch  gebundenen  Schichten  besitzen  al- 
lerdings  verschiedene  thermische  Ausdehnungskoeffizienten  oil  und  a2  ,  wodurch  kein 
ebener  Dehnungszustand  mehr  erreicht  wird.  Unter  diesen  Bedingungen  und  fur  den  Spe- 
zialfall,  dass  beide  Probenschichten  die  gleiche  Dicke  hi=h2  besitzen,  existieren  analytische 
Losungen  fiir  eine  stationare  („steady-state“)  Delamination  in  Abhangigkeit  von  der  thermi- 
schen  und  mechanischen  Belastung. 

Charalambides  [13]  ermittelte  einen  komplexen  mechanischen  Spannungsintensitatsfak- 
tor  fur  die  anliegende  mechanische  Last  P: 

K‘  =i(-L  +  i'W  (2-84) 

4{S  i  2  V 48 

wobei  t=3Pl/bh22  die  maximale  Normal spannung  hervorgerufen  durch  die  Biegung  der  unte- 
ren  Schicht  ist.  Nach  einigen  Berechnungen  ermittelt  man  fiir  die  mechanische  Belastung 
einen  Phasenwinkel  \|/  ~  40,89°. 

Der  Spannungsintensitatsfaktor  in  Bezug  auf  die  thermischen  Spannungen  aR=E(oc2-ai)AT 
wird  von  Charalambides  [13]  als  reiner  Mode  II  wie  folgt  bestimmt: 

Kr  =  i-oRh J  =  -oRh}ei(Kl2)  (2-85) 

4  4 

Der  thermische  Phasenwinkel  ist  i|/R=7u/2.  Die  thermischen  und  mechanischen  Spannungen 
bewirken,  dass  die  Spannungsintensitaten  sich  nicht  in  Phase  (\|/R^V|/t)  an  der  abzweigenden 
Delaminationsrissspitze  befinden. 

Fiir  die  Auswertung  solcher  Systeme  miissen  in  Abhangigkeit  von  Phase  komplizierte 
Bruchdiagramme  erstellt  werden,  um  eine  Aussage  treffen  zu  konnen,  ob  ein  Debonding 
(Dekohasion)  oder  Eindringen  des  Risses  in  die  zweite  Schicht  eintritt 
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2.10.7  Verhalten  von  Matrixrissen  bei  Rissfortschritt  in  geneigten  Fasern 

Matrixriss  in  einem  sproden  Verbund  mit  zur  Belastungsrichtung  geneigten  Fasem  (Faser- 
versagen  vor  Faser debonding) 

Im  Gegensatz  zu  Matrixrissen,  die  sich  unter  90°  zur  Faserrichtung  ausbreiten,  liegt  eine 
komplexere  Situation  in  einem  System  vor,  bei  dem  die  Faserverstarkung  in  Bezug  auf  die 
Matrixrissebene  geneigt  und  nebeneinanaer  angeordnet  sind.  Abb.  2-20  zeigt  hierzu  ein  Bei- 
spiel.  Bei  sproden  Verbundsystemen,  bei  dem  Faserversagen  gemaB  Abb.  2-20  gegeniiber 
Faserdebonding  gemaB  Abb.  2-21  bevorzugt  wird,  muss  der  Spannungsintensitatsfaktor  an 
der  mechanischen  Mode  I  Matrixrissspitze  unter  der  Anwesenheit  von  thermischen  Span- 
nungen  ein  Mode  I  bleiben. 


/  / 

1 


Abb. 2-20:  Faserversagen,  (Rissauftreffwinkel<90°) 

Dies  ist  besonders  fur  lange  Matrixrisse  der  Fall  (Risse  laufen  iiber  eine  groBe  Anzahl  von 
Fasem  hinweg),  wobei  die  durchschnittliche  thermische  Schubhaftung  auf  der  Rissoberfla- 
che  liber  die  gesamte  Risslange  Null  sein  muss. 

Als  Resultat  muss  ein  Faserversagenskriterium  unabhangig  von  den  thermischen  Spannun- 
gen  gefimden  werden. 
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Matrixriss  in  einem  sproden  Verbund  mit  zur  Belastunesrichtung  geneigten  Fasem  (Faserde- 
bonding  vor  Faserversagen) 

Andererseits,  wenn  Faserdebonding  vor  Faserversagen  gemaB  Abb.  2-21  auftritt,  werden  die 
in  Linie  laufenden  Fasem  den  Riss  iiber  eine  bestimmte  Lange  von  der  Matrixrissspitze  aus 
uberbriicken.  Debonding  von  geneigten  Fasem  ist  eher  wahrscheinlich,  und  es  kann  eine 
Mode  II  thermische  Komponente  an  der  ansonsten  M 

ode  I  Matrixrissspitze  existieren.  Unter  diesen  Bedingungen  wird  das  Faserdebonding  von 
den  thermischen  Spannungen  beeinflusst. 


Abb.2-21:  Faserdebonding,  (Rissauftreffwinkel<90°) 

Unter  Beriicksichtigung  der  o.a.  Beobachtungen  konnen  wir  nun  den  Prozess  der  anfangli- 
chen  Faserdebonding  und  weiterer  Debondingerweiterung  (Delamination)  in  faserverstarkten 
Verbundwerkstoffen  mit  sproder  Matrix  gemaB  Abb.  2-21  betrachten. 

Abb.  2-22  zeigt  den  Phasenwinkel  an  der  Debondingrissspitze  in  Abhangigkeit  von  der  auf- 
gebrachten  Last  t  und  den  thermischen  Spannungen  ctr  als  Funktion  der  normierten  Debon- 
dingrisslange  a/R. 
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Die  Kurve  fur  den  Phasenwinkel  \|r*,  der  sich  aufgrund  der  mechanischen  Spannungen  er- 
gibt,  ist  Ergebnis  einer  Finiten  Element  Berechnung  von  Charalambides  [13]. 

Die  Kurve  fur  den  Phasenwinkel  xj/*",  der  sich  aufgrund  des  thermischen  Spannungen  er- 
gibt,  gilt  fur  Faserverbundwerkstoffe,  deren  Interface  mit  einer  Restzugspannung  beauf- 
schlagt  wird.  Diese  Werte  sind  aus  einer  „steady  state“  Losung  aus  der  Literatur  genommen 
worden  [13],  wobei  der  Wert  bei  a/R=0  von  (\|/R)a=o=42,2°  die  Losung  fur  die  initiale  Rissab- 
zweigung  „kink  crack"  darstellt. 


Aus  Abb.  2-22  erkennt  man,  dass  bei  a-*0  der  Phasenwinkel  fur  den  90°  Debondingriss 
(„kink  crack"),  der  aufgrund  von  thermischen  bzw.  mechanischen  Lasten  initiiert  wird, 
gleich  groB  ist:  V)/R=\|/t. 

Dies  ergibt  sich  daraus,  dass  am  Hauptmatrixriss  die  Mode  I  Singularity  den  xiberwiegenden 
Anteil  darstellt. 

Aus  diesen  Ergebnissen  geht  hervor,  dass  im  Anfangsrissstadium  des  Matrixrisses  die  ther¬ 
mischen  Spannungen  eine  geringe  Rolle  spielen  und  nur  der  Mode  I  der  mechanischen  Be- 
lastung  eine  Rolle  spielt.  Hier  gelten  die  Annahmen  von  Fie  und  Hutchinson  [5],  Erst  mit  zu- 
nehmender  Risslange  und  damit  Ubergang  in  den  „steady  state“-Bereich  mit  a>5,0R  wird  der 
thermische  Einfluss  relevant. 
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Die  thermische  und  mechanischen  Phasenwinkel  sind  nicht  mehr  gleich,  sodass  fur  das  Auf- 
stellen  des  Energiekriteriums  auch  der  Einfluss  der  thermischen  Lasten  beriicksichtigt  wer- 
den  muss. 

Fur  die  Aussage,  ob  Faserversagen  oder  Rissabzweigung  am  Interface  auftritt,  muss  ein 
kompliziertes  Bruchdiagramm  erstellt  werden,  das  dann  Aussagen  iiber  die  Art  des  Rissfort- 
schrittes  zulasst. 

2.10.8  Zusammenfassung  der  Untersuchungen  von  Charalambides 

Die  von  Charalambides  [13]  durchgefuhrte  Analyse  zeigt,  dass  sich  thermische  Einfliisse 
durchaus  auf  die  Art  eines  Rissfortschrittes  auswirken  konnen.  Um  eine  Aussage  zu  erhalten, 
welche  Art  des  Rissfortschrittes  eintritt,  miissen  komplizierte  Bruchdiagramme  erstellt  wer- 
de,  die  dann  eine  Analyse  zulassen.  Der  Einfluss  der  thermischen  Faktoren  kann  gemaB 
Charalambides  fur  Matrixrisse,  die  sich  senkrecht  zur  Faserrichtung  ausbreiten, 
vemachlassigt  werden,  da  sich  hierbei  nur  der  Mode  I  auswirkt.  Dies  gilt  jedoch  nur,  wenn 
man  absolute  Flaftung  zwischen  Faser  und  Matrix  voraussetzt. 

Fur  den  hier  untersuchten  kohlenstofffaserverstarkten  Kohlenstoff  spielen  die  thermischen 
Einfliisse  vermutlich  eine  groBe  Rolle,  da  sie  hier  ein  Aufschrumpfen  auf  die  Faser  bewirken 
und  sich  somit  ein  Mixed  Mode  auspragen  kann.  Da  fur  den  verwendeten  kohlenstofffaser¬ 
verstarkten  Kohlenstoff  auch  Reibungseffekte  zwischen  Faser  und  Matrix  untersucht  wer¬ 
den,  ist  eine  Analyse  der  thermischen  Einfliisse  in  dieser  Arbeit  unabdingbar.  Dies  geschieht 
jedoch  nicht  mit  Hilfe  von  komplizierten  Bruchdiagrammen,  die  eine  Auswertung  des  Pha- 
senwinkels  von  thermischer  und  mechanischer  Last  voraussetzen,  sondem  anhand  von  Fini- 
ten  Element  Berechnungen. 
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2.11  Vergleich  der  vorgestellten  Theorien  unter  Beriicksichtigung  des 
kohlenstofffaserverstarkten  Kohlenstoffs 


Grundlegend  fur  die  Experimente  und  Berechnungen,  die  im  Rahmen  dieser  Arbeit  erbracht 
werden,  ist  die  Betrachtung  der  Aussagen  aus  den  theoretischen  Ansatzen.  Aus  diesem 
Grund  wird  nun  jede  Theorie,  die  zum  Anfang  dieser  Arbeit  ausgiebig  dargestellt  wurde, 
hinsichtlich  ihrer  Bedeutung  fur  den  kohlenstofffaserverstarkten  Kohlenstoff  ausgewertet. 

Auswertung  der  analytischen  Betrachtung  nach  He  und  Hutchinson  hinsichtlich  des  kohlen¬ 
stofffaserverstarkten  Kohlenstoff 

He  und  Hutchinson  [5]  haben  die  Ausbreitung  eines  Matrixrisses  in  einem  faserverstarkten 
Werkstoff,  der  eine  isotrope  Faser  und  Matrix  besitzt,  untersucht.  Fur  a=0  (Matrixelastizi- 
tatsmodul=Faserelastitzitatsmodul)  wird  der  Matrixriss  nach  Auftreffen  auf  das  Interface  um 
90°  abzweigen,  sobald  das  Verhaltnis  der  Bruchenergie  des  Interfaces  zu  der  Bruchenergie 
der  Faser  Gc‘/Gcf  groBer  als  1 'A  [25-27]  ist.  Bei  kohlenstofffaserverstarktem  Kohlenstoff  han- 
delt  es  sich  jedoch  um  einen  Werkstoff  mit  einer  orthotropen  Kohlenstofffaser.  Somit  kann 
anhand  der  orthotropen  Materialeigenschaften  der  Dundursparameter  a  nicht  direkt  ermittelt 
werden.  Um  dennoch  mit  Hilfe  der  analytischen  Ergebnisse  nach  He  und  Hutchinson  einen 
Ausgangswert  fur  das  Verhaltnis  der  Energiefreisetzungsraten  zu  bestimmen,  war  daher  eine 
Anpassung  notwendig.  Die  orthotrope  Kohlenstofffaser  wurde  als  isotrop  angenommen 
(Schubmodul  |ii  =  24000  MPa,  Vi=0,2).  Mit  |i2=  15000  MPa  und  der  Querkontraktionszahl 
v2=0,2  fur  die  Matrix  ergibt  sich  gemaB  Gleichung  2-9  ein  Dundursparameter  a  von  0,23. 
Der  Exponent  fur  die  Spannungssingularitat  X  wird  zu  0,49  bestimmt.  In  Abb.2-9  wird  ein 
Wert  von  Gd/Gp  von  0,248  fur  den  doppelt  abzweigenden  Riss  und  0,32  fur  den  einfach  ab- 
zweigenden  Riss  ermittelt. 

Eine  Rissabzweigung  des  Matrixrisses  am  Interface  wird  somit  in  dem  kohlenstofffaserver¬ 
starktem  Kohlenstoff  dann  erfolgen,  wenn  das  Verhaltnis  der  Bruchenergien 


G‘ 

<  0,248 

G/ 

ist.  Die  doppelte  Rissabzweigung  wird  gegeniiber  der  einfachen  bevorzugt  auftreten,  da  der 
Schwellwert  hierfur  eher  erreicht  wird. 
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Auswertung  der  Boundary  Element  Methode  nach  Petrovic 

Tullock,  Reimanis,  Graham  und  Petrovic  [9]  beschreiben  das  Problem  des  Rissfortschrittes 
eines  Matrixrisses  beim  Auftreffen  auf  ein  Interface  mit  Hilfe  einer  Boundary  Element  Me¬ 
thode.  Um  die  Energiefreisetzungsrate  zu  bestimmen,  wurden  innerhalb  dieser  Arbeit  zwei 
grundsatzliche  Ansatze  diskutiert. 

Die  erste  Moglichkeit  ist  die  Ermittlung  der  Energiefreisetzungsrate  anhand  der  Anderung 
der  Gesamtenergie  eines  Einzelfaserverbundes  wahrend  des  Rissfortschrittes  um  einen  in- 
krementellen  Betrag. 

Die  zweite  Moglichkeit  beschreibt  die  Methode  nach  Kaninen  [30].  Hierbei  wird  die  Arbeit 
ermittelt,  die  notwendig  ist,  um  einen  Riss  wieder  zu  schlieBen,  nachdem  er  sich  um  einen 
bestimmte  inkrementelle  Lange  erweitert  hat. 

Tullock,  Reimanis,  Graham  und  Petrovic  [9]  wahlten  zur  Ermittlung  der  Energiefreiset¬ 
zungsrate  die  Methode  nach  Kaninen  [30]  aus.  Die  erfolgreiche  Anwendung  der  numeri- 
schen  Methode,  namlich  der  Boundary  Element  Methode,  hat  gezeigt,  dass  sich  numerische 
Verfahren  sehr  gut  fur  die  Untersuchung  des  Rissfortschrittes  von  faserverstarktem  Werk- 
stoffen  eignen. 

Auswertung  der  analytischen  Betrachtung  nach  He,  Evans  and  Hutchinson  hinsichtlich  des 
kohlenstofffaserverstarkten  Kohlenstoffs 


He,  Evans  und  Hutchinson  [10]  ermittelten  auf  analytischem  Weg  die  Abhangigkeit  des  E- 
nergiefreisetzungsratenverhaltnisses  Gd/Gp  von  Restspannungen,  die  entweder  radial  auf  das 
Interface  wirken  oder  axial  und  damit  senkrecht  zum  Matrixriss  .  Das  geschilderte  Verhalten 
von  kohlenstofffaserverstarktem  Kohlenstoff  wahrend  der  Herstellung  kann  mit  Hilfe  des 
Ansatzes  zumindest  qualitativ  nachvollzogen  werden: 

Unmittelbar  nach  der  Karbonisierung  bis  ca.  1000°C  liegen  im  Interface  aufgrund  von 
Matrixschrumpfen  in  der  Kohlenstoffmatrix  radiale  Druckspannungen  vor.  Der  Wert  T|n  ist 
somit  negativ.  Der  Matrixriss  liegt  unter  Zugspannung,  was  mit  einem  positiven  r|t  in  Ver- 
bindung  steht.  GemaB  Abb.  2-15  nimmt  der  Wert  Gd/Gp  ab,  wenn  der  Wert  T|n  zunehmend 
negativer  wird,  was  aufgrund  einer  Zunahme  radialer  Druckspannungen  im  Interface  infolge 
von  Matrixschrumpfspannungen  verursacht  wird. 

Diese  Matrixschrumpfspannungen  bewirken  schlieBlich,  dass  der  Wert  Gd/Gp  das  Bruch- 
energieverhaltnis  Gc'/Gcf  unterschreitet.  Tendenziell  fuhrt  dies  zu  einer  groBeren  Anfalligkeit 
der  Faser  gegenuber  Risspenetration. 
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Die  zusatzliche  Warmebehandlung  nach  der  Karbonisierung  vermindert  die  Schrumpf- 
spannungen  wieder.  Der  Werkstoff  relaxiert.  Der  Wert  T|n  wird  somit  weniger  negativ.  Gj/Gp 
wachst  an,  was  mit  einer  zunehmenden  Tendenz  zur  Rissabzweigung  im  Faser- 
Matrixbereich  einhergeht.  Anhand  dieses  analytischen  Ansatzes  wird  zumindest  tendenziell 
das  Verhalten  des  kohlenstofffaserverstarkten  Kohlenstoffs  geklart.  Eine  quantitative  Aussa- 
ge  kann  jedoch  nicht  erfolgen,  da  auch  hier  wieder  der  Dundursparameter  aufgrund  von  i- 
sotropen  Werkstoffkomponenten  ermittelt  werden  muss,  was  allerdings  der  orthotropen  Koh- 
lenstofffaser  widerspricht. 

Auswertung  der  analytischen  Betrachtung  nach  Gupta  und  Martinez  Till  hinsichtlich  des 
kohlenstofffaserverstarkten  Kohlenstoffs 


Mit  Hilfe  der  Analyse  von  Gupta  und  Martinez  ist  es  moglich,  den  Einfluss  der  Orthotropie 
auf  das  Rissfortschrittsverhalten  von  C/C  zu  beriicksichtigen.  Hierfur  werden  folgende  Mate- 
rialkennwerte  (auch  fiir  die  quasidreidimensionale  Berechnung  des  Verhaltnisses  Gd/Gp) 
gemafi  Tab.  2-1  verwendet: 


Materialeigenschaften 

Kohlenstofffaser 

Glaskohlenstoffmatrix 

EH 

225000  MPa 

35000  MPa 

EJ_ 

15000  MPa 

35000  MPa 

HI-1 

15000  MPa 

14580  MPa 

v||l 

0,2 

0,2  =V_L|| 

Tab.  2-1:  Materialeigenschaften  von  kohlenstofffaserverstarktem  Kohlenstoff 


Die  Nachgiebigkeitsmatrix  fur  den  zweidimensionalen,  orthotropen  Fall  sieht  wie  folgt  aus: 


1 


0 


0 

0 

1 

/%± 


(2-86) 
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Die  dimensionslosen  Grofien  (a,  (3,  etc.)  ermitteln  sich  unter  Anwendung  der  Gleichungen  2- 
48,  2-60,  2-61,  2-62  und  2-63  auf  die  Materialwerte  gemafi  Tab.  2-1.  Daraus  ergeben  sich 
u.a.  Werte  (Tab.2-2)  fur  die  Bestimmung  des  Energiefreisetzungsratenverhaltnisses  Gd/Gp. 
Aus  der  graphischen  Darstellung  (Abb.  2-17)  kann  man  durch  Interpolation  zwischen  den 
Werten  fur  \\  von  0,1  und  0,02  fur  X, =0,067  und  dem  a  von  0,248  einen  Wert  fur  Gd/Gp 
von  0,1 18  aus  Abb.  2-17  (blaues  Kreuz)  ermitteln. 


Xi  (Faser) 

0,067 

p  (Faser) 

1,162 

a  (C/C) 

0,248 

Gd/Gp  (aus  Abb.  2-17) 

0,118 

Tab.  2-2:  Bestimmung  des  Energiefreisetzungsratenverhaltnisses 

Es  zeigt  sich,  dass  die  theoretischen  Werte  nach  Gupta  und  Martinez  [11,  22],  die  eine  Or- 
thotropie  der  Faser  beriicksichtigt,  deutlich  geringer  sind  als  die,  die  nach  der  Theorie  von 
He  und  Hutchinson  [5]  berechnet  wurden.  Offensichtlich  spielt  die  Orthotropie  der  Faser  ei¬ 
ne  wesentliche  Rolle  beim  Rissfortschritt  und  darf  somit  nicht  vemachlassigt  werden.  Dieser 
theoretische  Ansatz  behandelt  allerdings  keine  Schrumpf-  oder  thermische  Spannungen,  wie 
sie  wahrend  des  „Processing“  auftreten.  Die  Werte  sind  daher  ebenso  wie  die  von  He  und 
Hutchinson  [5]  ermittelten  Werte  als  AusgangsgroBe  zu  sehen  und  konnen  nur  eine  Tendenz 
anzeigen. 

Auswertung  der  analvtischen  Betrachtung  nach  Charalambides  hinsichtlich  des  kohlenstoff- 
faserverstarkten  Kohlenstoffs 


Charalambides  [13]  untersuchte  den  Einfluss  von  thermischen  und  mechanischen  Spannun¬ 
gen.  Hierbei  wurde  besonders  auf  thermische  Restspannungen  eingegangen,  die  beispielwei- 
se  nach  der  Herstellung  im  Faserverbundwerkstoff  auftreten.  Kohlenstofffaserverstarkter 
Kohlenstoff  unterliegt  unmittelbar  nach  der  Karbonisierung  bis  1000°C  dem  Einfluss  hoher 
Matrixschrumpfspannungen.  Im  Gegensatz  zu  den  Betrachtungen  nach  He,  Evans  und  Hut¬ 
chinson  wird  hier  nicht  der  resultierende  Spannungszustand  an  der  Rissspitze  (unmittelba- 
re  Spannungssingularitat)  betrachtet,  sondem  der  Spannungszustand  getrennt  nach  einer  me¬ 
chanischen  und  thermischen  Belastung  untersucht.  Hierbei  wird  unterschieden,  ob  sich  die 
thermische  Belastung  in  Phase  oder  in  Gegenphase  mit  der  aufgebrachten  mechanischen  Be¬ 
lastung  befmden. 
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Charalambides  stellt  in  Abhangigkeit  vom  Phasenwinkel  der  Belastungen  zwei  Gleichungen 
fur  die  Bestimmung  des  Energiefreisetzungsratenverhaltnisses  Gd/Gp  auf.  Die  beiden  Glei¬ 
chungen  2-81  und  2-82  dienen  somit  in  Abhangigkeit  des  Phasenwinkels  zwischen  thermi- 
scher  und  mechanischer  Spannungen  zur  Bestimmung  des  Rissverhaltens.  Um  ein  Problem 
zu  analysieren,  muss  daher  der  Phasenwinkel  fur  thermische  und  mechanische  Spannungen 
fur  einen  bestimmten  Belastungszustand  ermittelt  werden.  Anhand  von  Bruchdiagrammen  ist 
der  weitere  Rissverlauf  zu  bestimmen. 

Bei  Matrixrissen,  die  sich  senkrecht  zur  Faser  befinden,  befmden  sich  die  thermischen  Span¬ 
nungen  mit  den  mechanischen  Spannungen  nur  dann  in  Phase,  wenn  der  Matrixriss  von  einer 
dieser  beiden  Spannungen  alleine  hervorgerufen  wird.  Die  thermische  Spannung  und  die 
mechanische  Spannung  sind  aber  nicht  in  Phase,  wenn  sie  gleichzeitig  auf  den  Rissfortschritt 
des  Matrixrisses  einwirken  und  sich  iiberlagem.  Die  Abhangigkeit  vom  Phasenwinkel  wird 
noch  groBer,  wenn  Schub spannungen  im  Debondingbereich  auftreten.  Fur  diesen  Fall  sind 
dann  komplexe  Bruchdiagramme  notwendig,  anhand  derer  dann  die  Versagensart  Faserbruch 
oder  Debonding  im  Interfacebereich  ermittelt  wird. 

Fur  den  hier  behandelten  kohlenstofffaserverstarkten  Kohlenstoff  ist  die  Erstellung  eines 
solchen  Bruchdiagramms  aufgrund  der  Komplexitat  der  unterschiedlichen  Einfliisse  jedoch 
nur  unter  erheblichen  Aufwand  zu  erstellen.  Die  Finite  Element  Methode  stellt  hier  eine 
schnellere  Alternative  dar. 

Im  folgenden  wird  nun  die  Ermittlung  der  Energiefreisetzungsrate  mit  Hilfe  der  Methode  der 
Finiten  Elemente  dargestellt. 
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2.12  Bestimmung  der  Energiefreisetzungsrate  G  mit  Hilfe  der  Methode  der  finiten 
Elemente 

2.12.1  Allgemeines 

GemaB  Gleichung  2-8  besteht  die  Notwendigkeit,  die  Energiefreisetzungsrate  bzw.  die 
Bruchenergie  der  Faser  und  des  Interfaces  in  Abhangigkeit  von  bestimmten  Einfliissen 
(Temperatur  und  Restspannungen)  zu  bestimmen.  Um  die  Energiefreisetzungsraten  fur  den 
behandelten  C/C-Werkstoff  zu  ermitteln,  bedarf  es  der  Methode  der  Finiten  Elemente. 

Grundsatzlich  konnen  fur  die  Ermittlung  der  Energiefreisetzungsraten  unterschiedliche  Me- 
thoden  verwendet  werden. 

•  Ermitteln  der  Energiefreisetzungsrate  anhand  der  Differenz  der  Gesamtenergie  fur  zwei 
unterschiedliche  Risslangen  mit  Aa— >0 

•  Ermitteln  der  Energiefreisetzungsrate  nach  Kaninen  [30] 

•  Ermitteln  der  Steifigkeit  fur  zwei  unterschiedliche  Risslangen. 

•  Anderung  der  Steifigkeitsmatrix  wahrend  des  Rissfortschritts 

•  Bestimmung  des  J-Integrals 

Fur  die  Problemstellung  dieser  Arbeit  haben  sich  die  beiden  ersten  als  zweckmaBig  heraus- 
gestellt.  Diese  werden  im  Anschluss  ausfuhrlich  beschrieben.  Yorher  wird  allerdings  noch 
auf  die  grundsatzliche  Vorgehensweise  bei  der  Modellierung  des  Faserverbundwerkstoffes 
eingegangen. 

2.12.2  Modellierung  von  Faserverbundwerkstoffen  als  Einzelfaserverbund 

Fur  die  weiteren  Betrachtungen  des  Rissfortschritts  in  C/C -Faserverbundwerkstoffen  war  es 
deshalb  notwendig,  geeignete  Modelle  zu  schaffen.  Mit  Hilfe  dieser  Modelle  ist  es  moglich, 
eine  Beschreibung  der  Risserweiterung  in  Faserverbundwerkstoffen  speziell  in  kohlenstoff- 
faserverstarkten  Kohlenstoff  (C/C)  darzustellen. 
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Grundlage  aller  FEM-Rechnungen  war  die  Modellierung  von  Einzelfaserverbunden  mit  Hil- 
fe  des  PRE-/POST-PROCESSORS  PATRAN  und  anschlieBender  Berechnung  unter  Ver- 
wendung  des  FEM-Programms  MSC  NASTRAN  [31].  Um  Rechenzeiten  ohne  Verringerung 
der  Genauigkeit  zu  reduzieren,  wurde  aufgrund  der  Rotationssymmetrie  nur  eine  zweidimen- 
sionale  Rechnung  durchgefuhrt.  Die  Verwendung  von  CTRIAX  6-Elementen  (Ringelemen- 
ten)  erweitert  die  Rechnung  praktisch  zu  einer  quasidreidimensionalen. 


Rotationssymmetrie 


Z  y 


Vemetzung  mittels 
CTRIAX-6-Elementen 


Abb.  2-23:  CTRIAX-6  Element 

Folgende  Randbedingungen  wurden  verwendet:  Die  Knotenpunktsverschiebungen  in  der  Ro- 
tationsachse  des  Modells  wurden  in  y-  und  x,  r  -  Richtung  verhindert.  Der  untere  Rand  des 
Modells  wurde  nur  in  z-  und  y-Richtung  arretiert.  Der  Rissfortschritt  am  unteren  Modellrand 
wurde  durch  „Offnen“  der  Elemente  Knoten  fur  Knoten  realisiert.  Die  Krafteinleitung  er- 
folgte  durch  eine  konstante  Verschiebung  am  oberen  Modellrand.  Dazu  wurde  NASTRAN- 
seitig  ein  sogenanntes  Multiple  Point  Constraint  (MPC)  angebracht,  das  eine  konstante  z- 
Verschiebung  der  oberen  Knotenreihe  ermoglicht.  Die  Belastung,  die  auf  einen  einzelnen 
Knoten  (den  unabhangigen)  dieser  Reihe  angebracht  wird,  wird  so  uber  die  einzelnen  Knoten 
verteilt,  dass  eine  einheitliche  z-Verschiebung  stattfindet. 
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Die  Vemetzung  wurde  so  gewahlt,  dass  eine  Verdichtung  zum  Interface  sowohl  von  der 
Matrix  als  auch  von  der  Faser  her  stattfand.  Eine  Verdichtung  zum  Rissgrund  hin  wurde  e- 


benfalls  realisiert. 


^993^t 


Glaskohlenstoff- 

Matrix 

Kohlenstofffaser 


Konzentrischer 
Rundriss  in  der 
Glaskohlenstoff- 
Matrix 


Quasidreidimensionales  Modell 


Dreidimensionales  Modell 


Abb.2-24:  Modell 

Im  nachsten  Abschnitt  werden  nun  wie  bereits  erwahnt,  die  Methoden  dargestellt,  die  eine 
Ermittlung  der  Energiefreisetzungsraten  aus  den  Modellen  ermoglichen. 

2.12.3  Ermitteln  der  Energiefreisetzungsrate  bei  Rissfortschritt  anhand  der  Gesamt¬ 
energie  W  des  Modells 

Als  erstes  soil  die  Methode  beruhend  auf  der  Ermittlung  der  Gesamtenergie  dargestellt  wer¬ 
den.  Anhand  der  FEM  wurde  vor  und  nach  Rissfortschritt  um  eine  Lange  Aa  die  Gesamt¬ 
energie  W  des  Modells  ermittelt,  die  sich  aus  der  Summe  der  auBeren  und  der  inneren  Ener- 
gie  ergibt.  Die  Gesamtenergie  wird  bei  Verwendung  des  linear  elastischen  Solvers 
(NASTRAN,  SOL  101)  im  Ergebnis-File  (.f06-File)  ausgegeben. 
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Die  Energiefreisetzungsrate  pro  Einheitsdicke  ermittelt  sich  wie  folgt: 


G(a  + 


W(a  +  A  a)  -  W(a ) 
A  a 


(2-87) 


Fiir  das  axialsymmetrische  Zylindermodell  wird  die  Energiefreisetzungsrate  pro  Flachenan- 
derung  AA  (mit  AA=2ra:Aa)  berechnet. 

aa  (ffo  +  A aj-ma) 

2  AA 


NH 


j  j+1 


theoretischer  Rissfortschritt 


realer  Rissfortschritt 


Abb.  2-25:  theoretischer  und  realer  Rissfortschritt 


Die  FE-Methode  ermittelt  aufgrund  des  quadratischen  Verschiebungsansatzes  im  rotations- 
symmetrischen  6-Knoten  Dreieckselement  einen  linearen  Verlauf  der  Spannungen.  Folglich 
ergibt  die  Energiefreisetzungsrate  beim  „Offnen  eines  Elementes“  nicht  die  Energiefreiset¬ 
zungsrate,  die  dem  Rissfortschritt  von  dem  Knoten  j  zum  nachsten  Knoten  j+1  entspricht. 
Dies  ist  nur  der  Fall,  wenn  G  iiber  die  Knotenkrafte  und  der  Knotenpunktsverschiebungen 
ermittelt  wird. 
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Man  erhalt  somit  nur  die  Energiefreisetzungsrate,  die  dem  Rissfortschritt  von  Elementmitte 
zur  nachsten  Elementmitte  frei  wird.  Aa  stellt  hierbei  den  Rissfortschritt  in  radialer  Richtung 
dar,  wahrend  r  der  effektive,  momentane  Radius  von  der  Symmetrieachse  bis  zur  Rissspitze 
der  Ausgangsrisslange  ist. 

2.12.4  Ermitteln  der  Energiefreisetzungsrate  nach  der  Methode  nach  Kanincn 

Als  zweite  Moglichkeit,  die  Energiefreisetzungsrate  zu  bestimmen,  wird  in  diesem  Abschnitt 
auf  die  Methode  von  Kaninen  [30]  eingegangen.  Bisher  wurde  die  Energiefreisetzungsrate 
direkt  aus  dem  ,f06  File  gewonnen.  Dies  ist  aber  nur  bei  einer  linearen  Rechnung  moglich. 
Der  NASTRAN-Loser  SOL  101  gibt  die  Gesamtenergie  des  Modells  aus.  Im  weiteren  Ver- 
lauf  der  Arbeit  wird  auch  der  Einfluss  der  Reibung  im  Interface  untersucht  werden.  Dazu 
werden  NASTRANseitig  sog.  GAP  Elemente  verwendet  werden,  die  als  Loser  den  SOL  106 
Pur  nicht  lineare  Berechnungen  voraussetzen.  Dieser  errechnet  jedoch  nicht  die  Gesamtener¬ 
gie  des  Modells,  sodass  auf  eine  andere  Methode  zuriickgegriffen  wird,  die  im  folgenden  be- 
schrieben  wird. 


2.12.4.1  Grundgleiehungen 

Im  folgenden  wird  die  Methode  nach  Kaninen  beschrieben,  die  auf  der  Bestimmung  der  E- 
nergie  beruht,  die  notwendig  ist,  um  den  Riss  zu  schlieRen.  Die  Energiefreisetzungsrate  G 
kann  durch  folgende  Gleichung  ermittelt  werden: 


dW  dU 
dc  dc 


(2-89) 


In  Gleichung  2-89  ist  U  die  Dehnungsenergie  des  Korpers  pro  Einheitsdicke,  W  ist  die  Ar¬ 
beit  der  auBeren  Krafte  pro  Einheitsdicke  und  c  ist  die  Risslange.  Die  Beziehung  zwischen  G 
und  dem  Spannungsintensitatsfaktor  fur  ein  isotropes  Material  im  ebenen  Dehnungszustand 
wird  nach  Irwin  [29]  mit  Hilfe  des  RissschlieRintegrals  berechnet: 


G  = 


(1-v2) 


(1-v2) 


K 


2 

II 


E 


E 


(2-90) 
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Die  Methode,  die  hier  beschrieben  wird,  beruht  auf  der  direkten  Berechnung  des  RissschlieB- 
integrals  und  erlaubt  Ki  und  Kn  durch  eine  einzige  Analyse  zu  berechnen.  Sie  setzt  keine  Be¬ 
rechnung  von  Spannungen  voraus. 

2.12.4.2  Beschreibung  der  Methode 


Nach  Irwin  [29]  ist  die  freiwerdende  Energie  beim  Rissfortschritt  um  ein  infmitesimales  Ac 
gleich  der  Arbeit,  die  notwendig  ist,  um  den  Riss  zu  seiner  Ursprungslange  zu  schlieBen.  Bei 
der  Benutzung  eines  Polarkoordinatensystems  mit  dem  Ursprung  an  der  erweiterten  Riss- 
spitze  wird  diese  Annahme  in  folgender  Gleichung  ausgedriickt: 


1  ■  1  rAc 

G  =  lim -  (jJAc-r,0)v(r,7t)dr  +  lim -  ' 

Ac-»0  0  Ar  *0  y  Ac->0  0 \r  **0 


Ac->0  2  Ac 


2  Ac  Jo 


(Ac  -  r,0)u(r,K)dr 


(2-91) 

Hierbei  ist  G  die  Energiefreisetzungsrate,  o  und  T  die  Spannungen  in  der  Nahe  der  Rissspit- 
ze,  u  und  v  sind  die  Rissoffnungsverschiebungen  zwischen  den  Punkten  der  Rissoberflache. 
Ac  ist  die  Risserweiterung  an  der  Rissspitze.  Die  Integrale  sind  Ausdriicke  fur  die  Energie¬ 
freisetzungsrate  fur  den  Mode  I  und  Mode  II.  Folglich  kann  die  o.a.  Gleichung  auch  wie 
folgt  geschrieben  werden: 


G,  =  lim— (Ac-r,0)v(r,7t)dr  (2-92) 

Ac-*0  2  Ac  y 

Gn  =  lim  f  T  (Ac  -  r ,0)u(r ,  k) dr  (2-93) 

ac->o2Acj° 

Hierbei  sind  G|  und  Gn  die  Energiefreisetzungsraten  fur  den  Mode  I  bzw.  Mode  II.  Einen 
Ausdruck,  um  den  Zusammenhang  zwischen  Gi  und  Ki  zu  erhalten,  wird  durch  Einbeziehen 
Westergaard's  Gleichung  [32]  fur  die  Spannungen  und  Verschiebungen  im  Rissgrund  in  die 
o.a.  Gleichung  erreicht.  Nach  einer  Grenzwertbetrachtung  erhalt  man: 


Gj 


G 


II 


K 

—(p 

E 


Ki 


(p 


(2-94) 

(2-95) 
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Hierbei  sind  <p=l  fur  den  ebenen  Spannungs-  und  1-v2  fur  den  ebenen  Verzerrungszustand 
und  E  stellt  den  Elastizitatsmodul  dar. 

Die  Energiefreisetzungsrate  kann  nun  anhand  der  Knotenkrafte  und  Verschiebungen  ermit- 
telt  werden  (siehe  Abb.  2-26). 


Abb.2-26:  Energiefreisetzungsrate  nach  Kaninen 

Physikalisch  kann  das  Integral  als  die  Arbeit  interpretiert  werden,  die  notwendig  ist,  um  ei- 
nen  Riss  zu  schlieBen.  Diese  ermittelt  sich  aus  der  Knotenkraft  F  und  der  Verschiebungsdif- 
ferenz.  In  der  FEM-Berechnung  ist  die  Arbeit  die  Halfte  des  Produktes  aus  der  Kraft  bei  c 
und  d  und  der  Verschiebungen  zwischen  c  und  d,  nachdem  der  Riss  sich  um  einen  Knoten- 
abstand  erweitert  hat.  Dies  erfordert  zwei  Berechnungen: 

•  Ermittlung  von  Fc  und  Tc 

•  Bei  erweitertem  Riss  (Ac)  wird  die  Verschiebung  von  c  und  d  berechnet. 

G'=iS,2VF'(V‘-V")  (2‘96) 

Gn  =  Jim  —-  Tc  (uc  -  ud )  (2-97) 
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Fur  die  Berechnung  wird  ein  Wert  von  Ac  gewahlt,  der  in  der  Abb.  2-26  angegeben  ist.  Die 
Werte  von  F  und  T  entsprechen  den  y-  und  x-Kraften,  die  notwendig  sind,  um  c  und  d  zu- 
sammenzuhalten.  Bei  der  Berechnung  von  F  und  T  ist  es  notwendig,  dass  12  und  Ac  die  glei- 
che  Lange  (Ac  «  Ursprungsrisslange)  besitzen,  da  ansonsten  eine  Ausgleichrechnung  not¬ 
wendig  ist  [30],  Die  Elemente  1,2, 3, 4  sollten  ebenfalls  gleich  groB  sein. 

2.12.5  Zusammenfassung  der  beiden  Methoden 

Im  folgenden  wird  die  Methode  der  Gesamtenergie  fur  die  Falle  angewendet,  die  ei¬ 
ne  absolute  Haftung  zwischen  Faser  und  Matrix  voraussetzen.  NASTRANseitig  wird 
hier  immer  die  Gesamtenergie  des  Modells  berechnet. 

Die  Methode  nach  Kaninen  wird  im  folgenden  ausschlieBlich  fur  die  Falle  verwen- 
det,  bei  der  Nichtlinearitaten  beriicksichtigt  werden  mussen  (Reibung),  die  zwischen 
Faser  und  Matrix  Reibung  zulassen.  Diese  Methode  ist  aufwendiger.  Sie  muss  jedoch 
bei  den  reibungsbehafteten  Modellen  verwendet  werden,  da  hierbei  keine  Gesamt¬ 
energie  des  Modells  seitens  NASTRAN  ausgegeben  wird. 
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3.  Aufgabenstellung 

In  der  vorliegenden  Arbeit  soli  das  Rissausbreitungsverhalten  von  kohlenstofffaserverstark¬ 
tem  Kohlenstoff  beschrieben  werden.  Als  Grundlage  der  Arbeit  dienen  experimentelle  Er- 
gebnisse,  die  im  Rahmen  eines  DFG-Vorhabens  an  unterschiedlichen  C/C-Werkstoffen  er- 
zielt  wurden  [1]. 

Ziel  der  vorliegenden  Arbeit  ist  es  nun  hauptsachlich  mit  Hilfe  der  Methode  der  Finiten  E- 
lemente  zu  klaren,  welche  Werkstoffparameter  verantwortlich  fur  die  dargestellte  Festig- 
keitsentwicklung  von  C/C  sind.  Der  Zusammenhang  zwischen  Festigkeit  von  C/C  (Abb.2-4) 
und  den  Eigenschaften  der  Komponenten  des  Verbundes  (Faser,  Matrix  und  Faser-Matrix- 
Grenzflache)  ist  auBerst  komplex,  da  die  Eigenschaften  der  Komponenten  beim  Processing 
einen  stetigen  Wandel  erfahren. 

Zu  der  geringen  Festigkeit  nach  der  Karbonisierung  konnten  folgende  Phanomene  beitragen: 

•  Schrumpfspannungen  in  der  Matrix  vergroBem  die  Wahrscheinlichkeit  eines  Faserbru- 
ches 

•  Sprodes  mechanisches  Verhalten  der  Glaskohlenstoffmatrix  verursacht  Matrixrisse,  die  in 
die  Faser  penetrieren,  und  somit  zu  einer  Schwachung  des  Gesamtverbundes  fuhren 

•  Degradation  der  Faser  bewirkt  eine  Festigkeitsabnahme  der  Faser 

•  Perfekte  Haftung  zwischen  Matrix  und  Faser  verhindert  Abbau  von  Spannungen  und 
fuhrt  damit  zu  einer  Schwachung  des  Gesamtverbundes  unter  Mode  I-Belastung 

Die  nachfolgende  Erholung  nach  einer  zusatzlichen  Behandlung  oberhalb  von  1000  °C  konn- 
te  durch  folgende  Phanomene  zu  erklaren  sein: 

•  Relaxation  der  Matrixschrumpfspannungen 

•  Graphitierung  der  Glaskohlenstoffmatrix  wodurch  Rissabzweigen  ermoglicht  wird. 

•  Erholung  der  Faserfestigkeit  aufgrund  einer  verstarkten  Graphitierung  der  Faser 

•  Enthaftung,  sodass  Kraftiibertragung  zwischen  Faser  und  Matrix  nur  durch  Coulombsche 
Reibung  oder  andere  Effekte  (ionische,  elektrostatische,  chemische  Bindung  etc.)  mog- 
lich  ist. 
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4.  Berechnungen  und  Versuche 

4.1  Allgemeines 

Die  folgenden  Abschnitte  stellen  nun  die  Berechnungen  und  Versuche  dar,  die  notwendig 
waren,  um  die  herstellungsbedingten  Einfliisse  auf  das  Bruchverhalten  von  kohlenstofffaser 
verstarktem  Kohlenstoff  zu  ermitteln. 

Als  erstes  wird  dazu  ein  Finite  Element  Modell  entwickelt,  mit  dem  die  Eigenschaften  des 
kohlenstofffaserverstarkten  Kohlenstoffs  gut  wiedergegeben  werden.  Die  Giite  des  Modells 
wird  anhand  der  bereits  besprochenen  theoretischen  Ergebnisse  uberpriift. 

AnschlieBend  erfolgt  die  Ermittlung  der  Bruchenergie  der  Faser  Gcf  und  des  Interfaces  Gc‘ 
fur  den  kohlenstofffaserverstarkten  Kohlenstoff.  Dabei  werden  fur  die  Faser  Ergebnisse  aus 
einem  Druckscherversuch  [33-37]  und  fur  das  Interface  die  Ergebnisse  aus  einer  FEM- 
Berechnung  auf  Grundlage  der  theoretischen  Betrachtungen  von  Honjo  und  Shindo  [8]  bzw. 
Ochiai  und  Osamura  [7]  ausgewertet. 

Die  Hauptarbeit  wird  darin  bestehen,  Aussagen  iiber  das  Bruchverhalten  bzw.  den  Rissfort- 
schritt  des  Matrixrisses  nach  Auftreffen  auf  das  Interfaces  zu  erhalten.  Dazu  miissen  in  Ab- 
hangigkeit  von  verschiedenen  Eigenschaften,  die  sich  aufgrund  des  Processings  im  Werk- 
stoff  einstellen,  die  Energiefreisetzungsraten  fur  Penetration  Gp  und  Rissabzweigung  Gd  be- 
stimmt  werden. 

Sowohl  fur  die  Bestimmung  der  Bruchenergien  als  auch  der  Energiefreisetzungsraten  mit 
Hilfe  der  Finiten  Element  Methode  werden  dann  die  zwei  bereits  besprochenen  Methoden 
angewendet. 

Eine  Methode  nach  Kaninen  [30]  wird  im  folgenden  bei  Rechnungen  verwendet,  die  keine 
Haftung  zwischen  Faser  und  Matrix  voraussetzen.  Die  andere  Methode  berechnet  aus  der 
Differenz  der  Gesamtenergie  vor  und  nach  dem  Rissfortschritt  die  Energiefreisetzungsrate. 
Diese  Methode  eignet  sich  gut  fur  die  Berechnung  der  Bruchenergien  und  dient  zur  Bestim¬ 
mung  der  Energiefreisetzungsraten  unter  der  Annahme  einer  perfekten  Haftung  zwischen 
Faser  und  Matrix. 
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4.2  Modellierung  der  Risserweiterung  in  einem  zvlindrischen  Faser-ZMatrixmodell 

und  Verifikation  mit  bekannten  Losungen 

4.2.1  Allgemeines 

Grundlage  fur  die  weitere  Vorgehensweise  ist  die  Erstellung  eines  Modells  fur  die  Methode 
der  Finiten  Elemente.  Mit  Hilfe  dieses  Modells  konnen  im  weiteren  Verlauf  die  Eigenschaf- 
ten  des  kohlenstofffaserverstarkten  Kohlenstoffs  untersucht  werden.  Es  ist  jedoch  unabding- 
bar,  dieses  Modell  anhand  der  bereits  bekannten  Losungen  von  He  und  Hutchinson  [5]  auf 
seine  Richtigkeit  zu  uberpriifen. 

4.2.2  Generieren  eines  Modells  mit  dem  FEM  Pre-/Postprozessor  PDA/PATRAN 

Fur  die  weiteren  Betrachtungen  des  Rissfortschritts  in  kohlenstofffaserverstarktem  Kohlen¬ 
stoff  (C/C)  war  es  deshalb  notwendig,  ein  geeignetes  Modell  zu  schaffen.  Tullock  et  al.  [9] 
nutzte  fur  die  Nachrechnung  der  von  He  und  Hutchinson  [5]  ermittelten  Werte  die  Boundary 
Element  Methode.  Diese  Nachrechnung  bestatigte  die  Anwendbarkeit  der  BEM  fur  dieses 
Problem. 

In  dieser  Arbeit  wird  die  Finite  Element  Analyse  verwendet,  um  die  herstellungsbedingten 
Einfliisse  auf  den  Rissverlauf  zu  untersuchen.  Dazu  wird  ein  entsprechendes  Modell  erzeugt. 
Mit  Hilfe  dieses  Modells  soil  zuerst  nachgewiesen  werden,  dass  die  Finite  Element  Analyse 
ein  geeignetes  Verfahren  zur  Beschreibung  der  Risserweiterung  in  Faserverbundwerkstoffen 
speziell  in  kohlenstofffaserverstarkten  Kohlenstoff  (C/C)  darstellt.  Zunachst  wird  nachge¬ 
wiesen,  dass  mit  der  FEM  einfache  Rissprobleme,  fur  die  es  eine  analytische  Losung  gibt, 
beschrieben  werden  konnen,  wie  z.B.  fur  eine  Rundprobe  mit  konzentrischem  AuBenriss. 
Hierbei  wird  die  „Werkstoffkombination“  als  isotrop  und  homogen  angenommen.  Fur  den 
Rundstab  mit  AuBenriss  soli  die  FE-Losung  mit  der  analytischen  Losung  verglichen  werden: 

•  Abhangigkeit  der  Energiefreisetzungsrate  von  der  Risslange  a  und  Vergleich  mit  der  Lo¬ 
sung  eines  konzentrischen  Rissfortschritts  beim  Rundstab  [38] 

•  Das  Verhaltnis  der  Energiefreisetzungsraten  Gd/Gp  betragt  0,25  fiir  einen 
Dundursparameter  a  =0  (Fasermaterial  =  Matrixmaterial)  gemaB  He  und  Hutchinson  [5, 
25,26,  27], 
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Das  Finite  Element  Modell  wurde  hinsichtlich  seiner  Dimensionen  fur  die  Darstellung  eines 
Einzelfaserverbundes  aus  kohlenstofffaserverstarktem  Kohlenstoff  ausgelegt.  Der  Einzelfa- 
serverbund  hat  einen  AuBendurchmesser  von  8,62  |im  (Matrix).  Die  Faser  hat  einen  Durch- 
messer  von  6,95  pm.  Damit  wird  ein  Faservolumenanteil  von  65  %  realisiert.  Durch  die 
Verwendung  von  CTRIAX6-Elementen  wird  ein  quasidreidimensionales  Modell  geschaffen. 
Der  Rissfortschritt  am  unteren  Modellrand  wurde  durch  die  Aufhebung  der  Randbedingung 
z=0  („Offnen  der  Elemente“)  Knoten  fur  Knoten  realisiert.  Die  Modellhohe  von  20,83  jim 
wurde  gewahlt,  damit  Storspannungen  als  Folge  des  Risses  vor  dem  Modellende  ausklingen. 
Die  Krafteinleitung  von  (willkurlich  gewahlten)  100000N  erfolgt  durch  eine  konstante  "Ver- 
schiebung"  am  oberen  Modellrand,  die  mit  Hilfe  von  NASTRAN  durch  ein  Multi  Point 
Constraint  (MPC)  realisiert  wurde.  Die  Vemetzung  wurde  so  gewahlt,  dass  eine  Verdichtung 
zum  Interface  sowohl  von  der  Matrix  als  auch  von  der  Faser  her  stattfand.  Eine  Verdichtung 
zum  Rissgrund  ist  ebenfalls  vorhanden  (Abb.  4-1).  Um  die  Genauigkeit  der  Finite-Element- 
Berechnung  zu  ermitteln  und  die  erhaltenen  Energiefreisetzungsraten  bei  Rissfortschritt  mit 
den  analytischen  Werten  fur  die  Probe  aus  einem  isotropen  Werkstoff  mit  konzentrischer 
AuBenkerbe  zu  vergleichen,  war  es  notwendig,  dass  beide  Bereiche,  Matrix  und  Faser  die 
gleichen  Materialwerte  (Tab.4-1)  erhalten. 


Interface 


Vernetzung  mittels 
CTRIAX-6-Elementen 


Abb.4-1:  AusschnittvergroBerung  des  Modells  am  Rissgrund 
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Faser 

Matrix 

Elastizitatsmodul 

Ef=  35000  Mpa 

Em  =  35000  Mpa 

Querkontraktionszahl 

vf=  0,2 

Vm  =  0,2 

Wanneausdehnungskoeffozient 

df=  3,5*10-*  1/K 

a„  =  3,5*10^  1/K 

Tab.4-1 :  Materialwerte  fur  den  Rundstab 

Nach  dem  Abringen  der  Zugbelastung  wurde  mit  Hilfe  einer  Finite-Element-Rechnung  vor 
und  nach  Rissfortschritt  um  eine  Lange  Aa  die  Gesamtenergie  W  des  Modells  ermittelt,  die 
sich  aus  der  Summe  der  aufieren  und  der  inneren  Energie  ergibt.  Aus  der  Differenz  der  bei- 
den  Gesamtenergien  konnte  dann,  wie  bereits  erwiihnt,  jeweils  die  Energiefreisetzungsrate 
fur  den  inkrementellen  Rissfortschritt3ermittelt  werden. 

4.2.3  Vergleich  der  Energiefreisetzungsrate  G„  der  FEM-Rechnung  mit  der  Losung 
fiir  den  Rundstab  mit  konzentrischem  Riss 

In  der  Literatur  [38]  wird  fur  den  Rundstab  mit  konzentrischem  Riss  eine  analytische  Losung 
fur  homogenes,  isotropes  Material  angegeben.  In  der  Bruchmechanik  wird  anhand  dieser  Lo¬ 
sung  fur  diese  Zugprobe  der  Spannungsintensitatsfaktor  ermittelt.  Die  Abb.  4-2  zeigt  die 
Abmessungen  der  Rundprobe.  Die  Literatur  gibt  fur  diesen  konzentrischen  Riss  folgende 
Abhangigkeit  zwischen  Spannungsintensitatsfaktor  Ki  und  der  Ligamentlange  a  an: 


3  Unter  der  Berucksichtigung  der  von  He  und  Hutchinson  [5]  gemachten  Aussagen  beziiglich  der 
Grenzwertanalyse  der  Energiefreisetzungsraten  bei  Annaherung  ans  Interface  gemaB  2.6.5  kann  die 
Finite  Element  Methode  fur  die  Anwendung  der  Methode  der  Anderung  der  Gesamtenergie  erfolgen. 
GemaB  He  und  Hutchinson  ist  es  fur  die  Berechnung  von  Gd/Gp  namlich  unerheblich,  welche  inkre¬ 
mentellen  Rissfortschritt  am  Interface  gewahlt  wird.  Die  Wahl  der  inkrementellen  Rissfortschrittswei- 
te  bei  der  Bestimmung  der  Energiefreisetzungsraten  sind  fur  eine  geringe  Abweichung  des  Exponen- 
ten  der  Spannungssingularitat  A  von  'A  namlich  unerheblich.  Eine  Abhangigkeit  von  der  Rissfort- 
schrittsweite  a  verschwindet  namlich  fur  A=l/2  gemaB  Gleichung  2-14  und  2-18.  Der  Wert  A  wurde 
fur  den  kohlenstofffaserverstarkten  Kohlenstoff  mit  0,49  bestimmt.  Dieser  Wert  weicht  um  nur  2% 
von  0,5  ab.  Aus  diesem  Grund  ist  ein  signifikanter  Einfluss  von  der  Rissfortschrittsweite  a  nicht  gege- 
ben.  Die  Finite  Element  Berechnung  lasst  somit  bei  der  Methode  der  Gesamtenergie  gute  Ergebnisse 
erwarten. 
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Aus  dem  Spannungsintensitatsfaktor  K  wird  fur  homogenes  Material  die  Energiefreiset- 
zungsrate  nach  folgender  Gleichung  ermittelt: 


G  = 


K] 


(4-2) 


wodurch  mit  Hilfe  von  Gleichung  4-1  und  4-2  die  Abhangigkeit  zwischen  Energiefreiset- 
zungsrate  und  relativer  Risslange  aR/b  =l-a/b  bekannt  ist. 


Abb.4-2  Geometrie  der  Rundprobe 

Die  Abb.  4-3  stellt  die  ermittelten  Werte  der  analytisch  ermittelten  Energiefreisetzungsrate 
fur  den  Rundstab  und  der  aus  der  FEM-Rechnung  stammenden  Werte  gegeniiber.  Hierbei 
wurde  einmal  eine  Berechnung  mit  Hilfe  von  CTRIAX  6-Elementen  durchgefuhrt,  womit 
ein  quasidreidimensionales  Modell  beschrieben  wird.  Als  Referenz  wurde  zusatzlich  ein  3D- 
Modell  gegeniibergestellt. 

Grundsatzlich  zeigt  sich  eine  gute  Ubereinstimmung  zwischen  den  Werten  aus  der  analyti- 
schen  Losung  fur  den  Rundstab  und  den  FEM-Rechnungen.  Die  Werte  aus  der  analytischen 
Losung  fur  den  Rundstab  sind  jedoch  um  ca.  6%  im  Schnitt  hoher  als  die  aus  den  FEM- 
Rechnungen.  Weiterhin  zeigt  sich,  dass  die  Losungen  der  CTRLAX-6  eine  ausreichende  Ge- 
nauigkeit  im  Vergleich  zu  den  dreidimensionalen  Ergebnissen  aus  der  3D-HEX  Vergleichs- 
rechnung  bieten.  Somit  kann  die  zweidimensionale  (quasidreidimensionale)  Rechnung  fur 
die  weiteren  Berechnungen  genutzt  werden.  Auffallend  ist  hierbei  ein  ZICK-ZACK-formiger 
Verlauf.  Dies  ist  vermutlich  auf  die  Eigenschaften  der  Elemente  zuruckzufuhren. 
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aR 

FE-2D-CTRIAX6 

Rundstab 

FE  3D-HEX,  grab 

FE  3D-HEX,  fein 

0,0247 

23,47 

0,025 

35,68 

0,0747 

89,82 

0,075 

110,14 

0,1125 

110,31 

0,1246 

169,63 

0,125 

181,6 

191,97 

0,1625 

168,05 

0,1746 

260,5 

0,175 

260,2 

285,9 

0,215 

360,8 

0,2246 

368,2 

0,225 

385,86 

398,28 

380,8 

0,2746 

500,56 

0,235 

404,5 

0,245 

429,5 

0,255 

455,8 

0,265 

483,2 

0,275 

493,07 

537,86 

511,6 

0,3245 

668,8 

0,325 

717,3 

Tab.  4-2:  Vergleich  der  Energiefreisetzungsrate  ermittelt  rait  Hilfe  von  Gleichung  6-1  und 
der  FEM-Rechnung  fiir  isotropes  Material 


relative  Risslange  aR/b  [1] 


Abb.4-3:  Vergleich  der  analytischen  Losung  und  FEM-Rechnung 
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4.2.4  Uberpriifen  des  FEM-Modells  mit  den  Ergebnissen  von  He  und  Hutchinson 

GemaB  den  analytischen  Betrachtungen  von  He  und  Hutchinson  [5]  betragt  fur  einen  Dun- 
dursparameter  von  a=0  (Fasermaterial  =  Matrixmaterial)  das  Verhaltnis  der  Energiefreiset- 
zungsraten  Gd/Gp  ca.  14  [5,  25,  26,  27],  Dieses  Ergebnis  muss  ebenfalls  mit  den  vorher  be- 
schriebenen  Modellen  erreicht  werden,  bevor  eine  Erweiterung  auf  die  Werkstoffkombinati- 
on  C/C  moglich  ist.  Hierzu  sind  zwei  Arten  von  Rissverlaufen  mit  Hilfe  der  Methode  der  Fi- 
niten  Elemente  zu  berechnen: 

•  Rissverlauf  vom  Modellrand  (Matrixbereich)  und  DurchstoBen  des  Interfaces  in  den  Fa- 
serbereich  (Penetration) 

•  Rissabzweigen  ins  Interface;  der  Riss  andert  seine  Richtung  um  90°  und  verlauft  entlang 
des  Interfaces  weiter. 


Pro  infinitesimalen  Rissfortschritt  war  jeweils  eine  Rechnung  notwendig,  bei  der  die  ent- 
sprechende  Energiefreisetzungsrate  nach  Gleichung  2-87  und  2-88  ermittelt  wurde.  Fiir  die 
Ermittlung  der  Energiefreisetzungsrate  bei  Rissfortschritt  am  unteren  Modellrand  mussten 
die  Werte  aus  der  FEM-Berechnung  noch  mit  dem  Faktor  zwei  multipliziert  werden,  da  die 
ermittelte  Energie  im  Modell  nur  die  Halfte  der  tatsachlichen  Energie  darstellt.  Fiir  die  Er¬ 
mittlung  der  Energiefreisetzungsrate  beim  Abzweigen  des  Risses  ins  Interfaces  wird  ein  In- 
terfaceriss  durch  die  Verdopplung  der  jeweiligen  Knoten  in  diesem  Bereich  erzeugt.  Eine 
Multiplikation  mit  dem  Faktor  2  der  Energiefreisetzungsraten  fur  die  Rissabzweigung  ist 
nicht  erforderlich,  da  die  Energieanderung  AW  sich  als  Folge  der  beriicksichtigten  Risser- 
weiterung  AaR  nur  im  betrachteten  (halben)  Modell  abspielt. 


relative  Risslange  aR/b  [1] 

Energiefreisetzungsrate  Gp  bei  Penetration 
[J/m2] 

Energiefreisetzungsrate  Gti  bei 
Rissabzweigung  [J/m2] 

0,827 

330,92 

0,850 

362,16 

0,8740 

355,74 

0,899 

389,01 

0,926 

384,93 

0,830 

81,60 

0,854 

70,61 

0,879 

65,52 

0,905 

65,77 

Tab.4-3:  Energiefreisetzungsrate  fur  einen  auBengekerbten  Rundstab  (Werkstoff,  isotrop) 
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Die  folgende  Abb.  4-4  zeigt  die  Energiefreisetzungsrate  in  J/m2  fur  zwei  mogliche  Rissver- 
laufe.  Der  Riss  wachst  von  einer  relativen  Risslange  aR/b=0,72  bis  auf  aR/b=0,82  an.  Am  In¬ 
terface  aR/b=0,82  (Grenze  zwischen  Faser  und  Matrix,  die  aber  aus  gleichem  Material  beste- 
hen)  wurde  das  Problem  in  zwei  Rechnungen  aufgespaltet.  Die  Energiefreisetzungsrate  G 
wurde  fur  den  Fall  des  um  90°  abzweigenden  Risses  (Deflection)  bzw.  fur  den  weiterhin  pe- 
netrierenden  Riss  (Penetration)  ermittelt. 


relative  Risslange  aR/b[l] 


Abb. 4-4:  Energiefreisetzungsrate  bei  Penetration  und  Deflection 

4.2.5  Zusammenfassung  der  Verifikation  des  FEM-Modells  mit  bekannten  Losungen 

Die  Energiefreisetzungsrate  des  penetrierenden  Risses  ist  am  Interface  stetig,  da  sich  die  Ma- 
terialeigenschaften  nicht  andem.  Erst  bei  der  Erweiterung  des  Model  Is  in  den  folgenden  Ka- 
piteln  auf  ein  System  aus  zwei  Komponenten  wird  dort  in  der  Energiefreisetzungsrate  eine 
Unstetigkeit  auftreten.  Die  Abb.  4-4  vergleicht  die  Energiefreisetzungsraten  fur  die  folgen¬ 
den  Falle: 

•  Rissabzweigung  entlang  des  Interfaces  (Deflection) 

•  DurchstoBen  des  Interfaces  und  Eindringen  in  die  Faser  (Penetration) 
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Charakteristisch  ist,  dass  die  Energiefreisetzungsrate  bei  Penetration  um  den  Faktor  4  hoher 
gegeniiber  der  bei  Rissabzweigung  (Deflection)  liegt.  Weiterhin  ist  zu  bemerken,  dass  die 
Energiefreisetzungsrate  bei  Rissabzweigung  im  betrachteten  Bereich  bei  weiterem  Rissfort- 
schritt  abnimmt,  wohingegen  sie  bei  Penetration  leicht  zunimmt.  Vergleicht  man  die  Werte 
der  Energiefreisetzungsrate  von  Rissabzweigung  Gd  mit  denen  von  Penetration  Gp  so  erge- 
ben  sich  folgende  Werte: 

Gd  am  Interface  bei  Rissabzweigung:  8 1 ,59  J/m2 

Gp  am  Interface  bei  Penetration:  330,91  J/m2 

Der  Quotient  der  beiden  Energiefreisetzungsraten  aus  Deflection  und  Penetration,  Gd/Gp, 
wird  somit  mit  dem  Wert  0,246  ermittelt. 

Der  Wert  von  0,25,  der  aus  der  analytischen  Betrachtungen  von  He  und  Hutchinson  [5,  25, 
26,  27]  ermittelt  wurde,  konnte  anhand  dieser  FEM-Rechnung  fur  das  homogene  isotrope 
Material  mit  einem  Fehler  von  1,6%  bestatigt  werden.  Die  Genauigkeit  des  Ergebnisses 
zeigt,  dass  die  gewahlte  Risslange  AaR  in  unmittelbarer  Umgebung  um  das  Interface  fur  die 
Berechnung  Gd/Gp  ausreichend  ist. 

Somit  ist  nachgewiesen,  dass  das  Finite  Elemente  Modell  fur  einen  Dundursparameter  a=0 
(Fasereigenschaft  =  Matrixeigenschaft)  zu  zuverlassigen  Ergebnissen  fuhrt. 

In  den  anschliefienden  Kapiteln  werden  nun  folgende  Werte  bestimmt,  mit  denen  es  moglich 
ist,  Aussagen  uber  den  Rissfortschritt  in  kohlenstofffaserverstarkten  Kohlenstoff  zu  erhalten: 

if 

-  Ermitteln  der  Bruchenergie  der  Faser  Gc 

-  Ermitteln  der  Bruchenergie  des  Interfaces  Gc' 

-  Ermitteln  der  Energiefreisetzungsrate  fur  Rissabzweigen  Gd 

-  Ermitteln  der  Energiefreisetzungsrate  fur  Risspenetration  Gp 
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4.3  Ermitteln  der  Bruchenergie  G,.f  der  Kohlenstofffaser 

f  i  • 

Wie  bereits  erwahnt,  miissen  die  Bruchenergien  der  Faser  Gc  und  des  Interfaces  Gc  ermit- 
telt  werden,  um  eine  Aussage  iiber  die  Art  des  Rissfortschrittes  zu  erhalten.  Diese  Bruch¬ 
energien  werden  zu  diesem  Zweck  mit  den  Energiefreisetzungsraten  verglichen,  um  eine 

Aussage  zu  treffen,  ob  eine  Rissabzweigung  am  Interface  stattfmdet  oder  ein  Faserbruch 

f 

stattfmdet.  Zunachst  wird  hierzu  die  Bruchenergie  der  Faser  Gc  ermittelt. 

4.3.1  Allgemeines 

Ausgehend  von  Literaturdaten  wird  nun  die  Bruchenergie  fur  die  Kohlenstofffaser  bestimmt. 
Ochiai  und  Osamura  [7]  und  Honjo  und  Shindo  [8]  haben  eine  Berechnungsmethode  verof- 
fentlicht,  die  auf  der  Grundlage  der  „Shear  lag  analysis44  beruht. 

Anhand  von  Titancarbid-beschichteten  Siliziumcarbidfasem  und  Titancarbid-beschichteten 
Kohlenstofffasem  wurde  die  Spannung  ermittelt,  die  zu  Faserbruch  fuhrt. 

Die  TiC-Schichtdicken  der  SiC-Faser  lagen  im  Bereich  von  0,003  pm  und  waren  deutlich 
geringer  als  die  der  C-Fasem  (iiber  0,010  pm).  Hierbei  zeigte  sich  bei  den  Kohlenstofffasem 
eine  Abhangigkeit  der  Bmchspannung  von  der  Schichtdicke  des  TiC.  Die  Abnahme  der  Fa- 
serfestigkeit  zeigte  folgende  Proportionalitat  bei  der  Kohlenstofffaser: 

®uts  ~  j  . .  (4'3) 

yanc 

Dagegen  wurde  bei  den  Siliziumcarbidfasem  keine  Abhangigkeit  der  Bmchspannung  von 
der  TiC-Schichtdicke  bestimmt.  Nach  Entfernen  der  TiC-Schicht  stellte  sich  bei  der  C-Faser 
die  urspriingliche  Faserfestigkeit  der  unbehandelten  Ausgangsfaser  ein. 

Die  Erklarung  fur  dieses  Verhalten  liegt  darin,  dass  bei  einer  guten  Haftung  zwischen  Faser 
und  Matrix  ein  Bmch  der  Beschichtung  zu  einem  vorzeitigen  Faserversagen  fuhrt.  Dieser  Ef- 
fekt  tritt  jedoch  nur  auf,  wenn  eine  kritische  DefektgroBe,  die  hierbei  durch  die  Schichtdicke 
aTic  bestimmt  wird,  iiberschritten  wird. 
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4.3.2  Uberlegiingen  nach  Honjo  und  Shindo  [81 

Versagen  der  Faserbeschichtung  („Fibercoating“)  infolge  seiner  geringen  Bruchdehnung 
kann  fruhzeitig  Faserbruch  auslosen.  Die  Festigkeit  der  Faser  Ofu  ist  in  diesem  Fall  eine 
Funktion  der  Dicke  a  der  Faserbeschichtung: 


<7 


/« 


1 

1.12 


Efil_ 

Jta 


fur  a>ac 


(4-4) 


Nach  Substitution  der  kritischen  GroBe  ac  ergibt  sich  folgender  Ausdruck  fur  ac: 


a „  = 


1.12(7® 

V  J 


e,g! 

K 


(4-5) 


o°ft,  stellt  hierbei  die  Faserfestigkeit  ohne  Faserbeschichtung,  Ef  den  Elastizitatsmodul  und 
Gcf  die  kritische  Energiefreisetzungsrate,  die  sog.  Bruchenergie  der  Faser  dar. 

Wenn  die  Dicke  der  Faserbeschichtung  kleiner  als  die  kritische  GroBe  ac  ist,  bleibt  die  Faser¬ 
festigkeit  unbee  intrachtigt.  Es  gilt  folgende  triviale  Beziehung: 


&  fit  ®  fiu 


fur  a<ac 


(4-6) 


Die  kritische  Dicke  des  Faseriiberzuges  ac  reprasentiert  den  Grenzwert,  unterhalb  dem  die 
Faser  aufgrund  von  inneren  Defekten  (z.B  Fehlstellen,  Porenbildung  etc.)  versagt  und  ihre 
Festigkeit  gleich  der  Festigkeit  der  unbehandelten  Faser  ist. 

Die  Gleichung  4-4  und  4-5  sind  frir  den  ebenen  Spannungszustand  dargestellt.  Die  Gleichun- 
gen  4-4  bis  4-6  gelten  fur  den  Spezialfall,  dass  die  elastischen  Eigenschaften  der  Faserbe¬ 
schichtung  identisch  zu  denen  der  Faser  sind.  Wenn  die  Eigenschaften  des  Faseriiberzuges 
sich  von  denen  der  Faser  unterscheiden,  ist  die  Bestimmung  der  Faserfestigkeit  mit  Hilfe  der 
o.a.  Gleichung  prinzipiell  nicht  moglich. 
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Nach  Osamura  und  Ochiai  [7]  fuhren  ein  grofier  E-Modul  und  eine  groBe  Schichtdicke  des 
Faseriiberzuges  zu  einer  erheblichen  Reduktion  der  Faserfestigkeit.  Wie  bereits  erwahnt 
wurde  anhand  von  TiC-beschichteten  Kohlenstoff-  bzw.  Siliziumcarbidfasem  diese  Gesetz- 
maBigkeit  bestimmt  [7].  Dabei  wurde  allerdings  die  Faser  als  isotrop  angenommen.  Auch 
thermische  Spannungen  wurden  nicht  beriicksichtigt. 

Im  folgenden  wird  die  Finite  Element  Methode  fiir  die  Berechnung  der  Faserfestigkeit  von 
beschichteten  orthotropen  Kohlenstofffasem  unter  Beriicksichtigung  thermischer  Span¬ 
nungen  verwendet.  Anhand  der  Ergebnisse  von  Ochiai  und  Osamura  [7]  erfolgt  zuvor  eine 
Verifizierung.  AnschlieBend  wird  die  Methode  auf  die  von  mir  untersuchten  Kohlenstofffa¬ 
sem  umgesetzt. 

4.3.3  Verifizierung  des  finiten  Element  Modells  anhand  der  Literatur 

Im  folgenden  wird  ein  Finite  Element  Modell  entwickelt,  um  die  von  Osamura  und  Ochiai 
[7]  durchgefuhrten  Berechnungen  nachzuvollziehen  und  anschlieBend  auf  die  orthotrope  Fa¬ 
ser  unter  Beriicksichtigung  der  thermischen  Spannungen  zu  erweitem.  Mit  Hilfe  dieser  Er- 
weiterung  lasst  sich  dann  iterativ  sowohl  die  Bruchenergie  der  Faser  als  auch  die  kritische 
Risslange  ermitteln. 

4.3.4  Verifizierung  der  Berechnungen  von  Ochiai  und  Osamura 

FEM-Modelle  von  TiC-beschichteten  Kohlenstofffasem  (Faserradius  3,475  pm)  mit  mehre- 
ren  Schichtdicken  wurden  erzeugt.  Wie  bereits  beschrieben  wurde  hierbei  ein  Einzelfaser- 
verbund  als  quasizweidimensionales  Modell  entwickelt.  Die  Ergebnisse  von  Ochiai  [7]  ge- 
maB  Abb.  4-5  und  die  Materiaidaten  nach  Tabelle  4-4  dienten  als  Ausgangsdaten. 


TiC-Beschichtung 

C-Faser 

E  =  440000  MPa 

E  =  225000  MPa 

< 

11 

© 

v  =  0,2 

a  =  7,4*1  O'6 1/K 

a=  1‘IQ-6  1/K 

Tab.  4-4:  Materiaidaten  der  TiC-Beschichtung  und  der  C-Faser  gemaB  Ochiai  [7] 
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3080  MPa 

3500  (unbeschichtctc  C’-Faser) 


0  0,02  0,04  .  0,06  0,68 

0.0?  uni  0.0S  um 


Abb.  4-5:  Bruchspannung  in  Abhangigkeit  der  TiC-Schichtdicke  nach  Ochiai  [7] 

Zunachst  wurde  eine  Beschichtung  aus  TiC  von  0,08  jam  Dicke  analysiert.  Anschliefiend 
wurde  ein  Zugversuch  mit  einer  Spannung  von  800  MPa  (entspricht  der  Zugfestigkeit  dieser 
Schichtdicke  [7],  angenommene  Daten  sind  Mittelwerte)  fiir  diese  beschichtete  Faser  be- 
rechnet. 


Abb.4-6:  Energiefreisetzungsrate  fair  C-Faser  mit  TiC-Schicht  ohne  thermische  Spannungen 
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Mit  der  bereits  bekannten  Methode,  namlich  Knoten  fur  Knoten  zu  offnen,  wurde  der  Riss 
durch  die  TiC-Beschichtung  in  die  Faser  vorangetrieben.  Anhand  der  Unterschiede  in  der 
Gesamtenergie  wahrend  des  Rissfortschritts  wurde  die  Energiefreisetzungsrate  bestimmt  und 
liber  der  relativen  Risslange  a/a,  (in  Bezug  auf  die  Coatingdicke  von  0,08  pm)  aufgetragen 
(Abb.  4-6). 

Mit  zunehmender  Risslange  steigt  zunachst  die  Energiefreisetzungsrate  an.  Beim  Eindringen 
in  die  C-Faser  besteht  eine  Singularity.  Der  Wert  der  Energiefreisetzungsrate  geht  gegen  un- 
endlich.  Nach  Eindringen  in  die  TiC-Schicht  erfolgt  eine  Abnahme  der  Energiefreisetzungs¬ 
rate.  Erst  bei  groBeren  Risslangen  kann  wieder  ein  Anstieg  der  Energiefreisetzungsrate  er- 
wartet  werden. 

Durch  Interpolation  in  der  Nahe  des  Interfaces  (a/a,=  1 )  erhalt  man  aus  Abb. 4-6  einen 
Wert  von  2,6  J/m2  als  erste  Iteration  fur  die  Bruchenergie. 

Ochiai  [7]  bestimmte  die  Bruchenergie  der  Kohlenstofffaser  anhand  der  Festigkeit  der  unbe- 
schichteten  Faser  mit  Hilfe  des  folgenden  Diagramms: 


a/r  [1] 

Abb.4-7:  Energiefreisetzungsrate  fur  eine  unbeschichtete,  isotrope  Kohlenstofffaser 

Fur  die  Bruchenergie  von  2,6  J/m2  ergibt  sich  aus  Abb.4-7  eine  relative  kritische  Risslange 
von  a/r=0,05  bzw.  eine  kritische  Risslange  von  0,174  pm  bei  r=3,475  pm.  Nach  Ochiai  und 
Osamura  [7]  muss  diese  kritische  Risslange  im  bcschichteten  Fall  (Abb.  4-6)  ebenfalls  zu 
einer  Schadigung  fuhren. 

Daher  ergibt  sich  fur  die  TiC-behaftete  C-Faser  eine  relative  Gesamtrisslange  von 

0,174;Uw-l-0,08jUw  _  2 
0,08  /im 
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Somit  wurde  sich  fur  die  relative  Gesamtrisslangen  von  3,2  eine  Bruchenergie  der  C-Faser 
(ohne  Berucksichtigung  von  Temperatureinfliissen  und  Orthotropie)  von  2,3  J/m2  ergeben 
(Abschatzung  durch  Extrapolation).  Der  nachsten  Iterationsschrittsschritt  ergabe  mit  dieser 
Bruchenergie  aus  Abb.  4-7  eine  relative  kritische  Risslange  von  0,044  bzw.  eine  kritische 
Risslange  von  0,153  pm. 

4.3.5  Ermitteln  der  Bruchenergie  fur  die  orthotrope  Kohlenstofffaser  unter  Beruck¬ 
sichtigung  thermischer  Spannungen 

Im  folgenden  soli  nun  die  o.a.  Methode  fur  die  Berechnung  der  Bruchenergie  der  orthotropen 
Kohlenstofffaser  genutzt  werden.  Hierbei  soil  im  Gegensatz  zu  der  Methode  nach  Ochiai  der 
herstellungsbedingte  Einfluss  thermischer  Spannungen  beriicksichtigt  werden. 

Fur  eine  TiC-Beschichtung  von  0,05  pm  bzw.  0,08  pm  und  die  entsprechende  unbeschichte- 
te  Faser  wurden  FEM-Berechnungen  gemaG  o.a.  Iterationsschema  durchgefuhrt,  um  Aussa- 
gen  iiber  die  kritische  Risslange  und  Bruchenergie  der  orthotropen  Kohlenstofffaser  zu  erhal- 
ten.  Die  Zugspannungen  entsprechen  den  Bruchlasten,  die  von  Ochiai  [7]  ermittelt  wurden 
und  aus  Abb.  4-5  ersichtlich  sind.  In  den  folgenden  Abbildungen  wird  der  Verlauf  der  Ener- 
giefreisetzungsraten  bei  Rissfortschritt  fur  die  drei  Falle  dargestellt. 


Abb.4-8:  Energiefreisetzungsrate  fur  die  TiC-beschichtete  orthotrope 
Kohlenstofffaser  (Schichtdicke  0,05  pm) 
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Fiir  die  beiden  Beschichtungsdicken  wird  nun  das  bereits  beschriebene  Iterationsverfahren 
durchgefuhrt,  sodass  die  kritische  Risslange  und  damit  die  DefektgroBe,  die  zum  Faserbruch 
fuhrt,  ermittelt  werden  kann.  Aus  der  Methode  wird  auch  die  Bruchenergie  Gcf  der  orthotro- 
pen  Faser  erhalten. 


TiC-Beschichtung  0,05  |im 

1.  Iteration 

2.  Iteration 

3.  Iteration 

Bruchenergie  Gc*  (Abb.4-6) 

21  J/m2 

15  J/mJ 

15  J/m2 

relative  kritische  Risslange  (Abb.4-8) 

0,0162 

0,011 

0,011 

kritische  Risslange  (Abb.4-8) 

0,056  |im 

0,038  gm 

0,038  gm 

relative  kritische  Risslange  fur  beschichtete 

Faser  (Abb.4-6) 

2,12 

1,475 

1,475 

Tab.4-5:  Iterationen  fur  eine  TiC-Beschichtung  von  0,05  pm 


TiC-Beschichtung  0,08  gm 

1.  Iteration 

2.  Iteration 

3.  Iteration 

Bruchenergie  Gc‘  (Abb.4-7) 

26,5  J/m2 

15  J/m2  ”  ~ 

15  J/m5  “ 

relative  kritische  Risslange  (Abb.4-8) 

0,0195 

0,011 

0,011 

kritische  Risslange  (Abb.4-8) 

0,068  gm 

0,038  gm 

0,038  gm 

relative  kritische  Risslange  fur  beschichtete 

Faser  (Abb.4-7) 

1,85 

1,475 

1,475 

Tab.4-6:  Iterationen  fur  eine  TiC-Beschichtung  von  0,08  pm 


4.3.6  Zusammenfassung  der  Rechnungen  zur  Bestimmung  der  Bruchenergie  der  Fa¬ 
ser 

Es  zeigt  sich,  dass  die  kritische  Risslange  und  die  Bruchenergie  fur  beide  Schichtdicken 
gleich  ist.  Aus  dem  iterativen  Verfahren,  das  aus  Berechnungen  nach  Honjo  und  Shindo  [8] 
bzw.  Ochiai  und  Osamura  [7]  entwickelt  wurde,  wird  die  Bruchenergie  Gcf  fur  die  Kohlen- 
stofffaser  mit  15  J/m2  ermittelt.  Die  orthotrope  Kohlenstofffaser  wird  dann  unter  Mode  I  (Fa¬ 
serbruch  unter  90°)  versagen,  wenn  bei  einer  kritischen  Risslange  von  0,038  pm  die  Energie- 
freisetzungsrate  Gp  die  Bruchenergie  Gcf  von  15  J/m2  iiberschreitet: 

Gn  >  Gl  mit  Gcf  =  15  J/m2  (4-7) 

P  c 
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Der  Unterschied  zu  dem  Wert  2,3  J/m2,  der  mit  Hilfe  von  Ochiai  ermittelt  wurde,  lasst  sich 
dadurch  erklaren,  dass  Ochiai  die  Orthotropie  der  Faser  vemachlassigt  hat.  Weiterhin  wur- 
den  bei  der  Analyse  von  Ochiai  thermische  Spannungen  aufgrund  der  unterschiedlichen 
Warmeausdehnungskoeffizienten,  die  beim  Abkiihlen  auftreten,  nicht  beriicksichtigt. 

Die  Literatur  [1]  gibt  als  Obergrenze  fur  die  Risszahigkeit  Kc  einer  „gewohnlichen“  Kohlen- 
stofffaser  einen  Wert  von  1  MPaVm  an.  Damit  ergibt  sich  fur  eine  orthotrope  Kohlenstofffa- 
ser  nach  [39, 40]  folgende  Rechnung  unter  Annahme  eines  ideal  elastischen  Verhaltens: 

G{  =  ~  (4-8) 

E.l 


mit  Eeff  zu 


mit  den  by  aus 


1  _  1^11^22  j 

Eeff  V  2  V 


bn  =qn -  — 

(4-10) 

#33 

u  _  „  #13#23 

bn  -  Qn 

(4-11) 

#33 

2 

h  —  n  ^23 

°22  #  22 

(4-12) 

q  33 

C  -  '/  66 

(4-13) 

W-l eh] 

(4-14) 

J22  2^12  +  b66 


2b, . 


(4-9) 


Q  stellt  hierbei  die  Nachgiebigkeitsmatrix  fur  den  anisotropen  Fall  dar. 
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If  xy 


-v. 


-V 


_1^ 


0 

0 

0 

J_ 


0  0 

0  0 
0  0 


0 


0 

0 

1 

Vxy 


yz 


(4-15) 


Unter  Verwendung  der  Materialkonstanten  aus  Tab.  2-1  ergeben  sich  die  Werte  by  zu: 


bn  = 

1,77*  10"6  1/MPa 

b,2  = 

-1 6*  1 0“6  1/MPa 

b22  = 

64*  10'6  1/MPa 

b66  = 

66*1 0‘6  1/MPa 

Tab.  4-7  Materialkonstanten  der  Nachgiebigkeitsmatrix  fiir  die  orthotrope  Kohlenstofffaser 

Der  Wert  1/Eeff  ermittelt  sich  daraus  gemaB  o.a.  Gleichung  zu  29,74*  10"6  1/MPa.  Als  Ober- 
grenze  fiir  die  Bruchenergie  einer  Kohlenstofffaser  gemaB  Literatur  kann  somit  30  J/m2  an- 
genommen  werden.  Die  von  mir  ermittelten  Werte  fiir  die  Kohlenstofffaser  Tenax  HTA  lie- 
gen  mit  15  J/m2  in  dem  Bereich,  der  von  der  Literatur  grob  abgeschatzt  wird. 
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4.4  Ermitteln  der  Bruchenergie  Gf‘  anhand  eines  Druckscherversuches 

Um  Aussagen  iiber  den  Rissfortschritt  zu  erhalten,  muss  das  Verhaltnis  der  Bruchenergien 
Gc‘/Gcf  bestimmt  werden.  Die  Bruchenergie  der  Faser  Gcf  wurde  bereits  ermittelt. 

Um  nun  den  Wert  fur  die  Bruchenergie  des  Interfaces  zu  bestimmen,  musste  eine  Methode 
gewahlt  werden,  die  eine  Schubspannungsverteilung  im  Interface  erzeugt. 

Die  im  folgenden  verwendeten  Proben  bestehen  aus  2D  C/C  (8-H-Satin)  Gewebe  in  einer 
Kohlenstoffmatrix.  Ein  erheblicher  Vorteil  der  unsymmetrisch  gekerbten  Druckscherprobe 
besteht  darin,  dass  keine  Zugnormalspannungen  senkrecht  zu  den  Fasem  auftreten.  Lediglich 
Drucknormalspannungen  entstehen  senkrecht  zur  Faser.  Daher  wurde  der  Wert  Gc‘  im  fol¬ 
genden  aus  einer  unsymmetrischen  Druckscherprobe  [33,  34,  35,  36,  37]  ermittelt. 


3,5  mm 


Gewebelagen 


Kerbe 


interlaminarer 
Riss 


Hingeleitete  Druckkraft  bewirkt  eine 
interlaminare  Schubspannung,  die  zu 
einer  Verlangerung  der  initiierten 
interlaminaren  Risse  an  den  Kerben 
fuhrt. 


Abb.  4-11:  Druckscherversuch 
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Beim  Druckscherversuch  [33]  wird  eine  3,5  mm  breite  C/C-Probe  verwendet,  die  unsymmet- 
risch  an  jeder  Seite  mit  einer  Kerbe  von  1,75  mm  Lange  versehen  ist.  Die  eingeleite  Druck- 
kraft  induziert  einen  interlaminaren  Riss,  der  jeweils  an  der  Spitze  der  Kerben  beginnt.  Der 
Riss  pflanzt  sich  zwischen  den  Schichten,  entweder  in  der  Kohlenstoffmatrix  oder  in  der  Fa- 
ser-/Matrix-Grenzflache  fort.  Um  die  Bruchenergie  des  Interfaces  zu  bestimmen,  wurde  fol- 
gender  Ablauf  durchgefuhrt,  der  in  Abb.  4-12  schematisch  dargestellt  ist. 

Anhand  eines  Dreipunktbiegeversuches  wurden  Initialrisse  (interlaminare  Risse)  in  der  a- 
symmetrisch  gekerbten  2D-C/C-Probe,  im  folgenden  als  Originalprobe  bezeichnet,  erzeugt. 
Da  sich  die  mikroskopische  Risslangenermittlung  schwierig  gestaltet,  war  eine  Kalibrierkur- 
ve  fair  die  Abhangigkeit  Compliance/Risslange  zu  erstellen.  Diese  Kalibrierkurve  wurde  er¬ 
zeugt,  indem  bestimmte  Risslangen  a  aus  einer  Vergleichsprobe  ausgesagt  wurden  und  die 
zugehorige  Compliance  dieser  Probe  ermittelt  wurde.  AnschlieBend  wurde  der  Druckscher- 
versuches  an  der  Originalprobe  durchgefuhrt  und  deren  Compliance  ermittelt.  Die  zugehori¬ 
ge  Risslange  wurde  anhand  der  Kalibrierkurve  der  Vergleichsprobe  ermittelt  (Abb.4-12). 


Ermitteln  einer  Kalibrierkurve:  Compliance  C  in 
Abhangigkeit  von  der  Risslange  a  in  2D-C/C 


Rissinitiierung  an  jeder  Kerbe  durch  Dreipunktbiegung  (2D- 
C/C-Probe),  Druckkraft  der  Kerbe  jeweils  gegeniiber 


Druckscherversuch  ergibt  Compliance  C  fur  den 
Rissfortschritt  des  interlaminaren  Risses 


Ermitteln  der  zugehorigen  Risslange  a  des  inter¬ 
laminaren  Risses  anhand  der  Kalibrierkurve 


Berechnen  der  Bruchenergie  des  Interfaces  Gc'  fur 

C/C  karbonisiert  und  C/C  graphitiert 


Abb.4-12:  Ablaufdiagramm  zur  Bestimmung  des  Interfacebruchenergie 
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4.4.1  Rechnerische  Bestimmung  der  Bruchenergie  aus  dem  Druckscherversuch 

Fur  die  Bestimmung  der  Bruchenergie  aus  dem  Druckscherversuch  ist  es  notwendig,  einen 
Zusammenhang  zwischen  der  Probennachgiebigkeit  (Compliance),  die  sich  aus  dem  Kraft- 
verlangerungsdiagramm  ergibt,  der  Interfacerisslange  a  und  der  Bruchenergie  Gcl  aufzustel- 
len.  Die  Energiefreisetzungsrate  ist  die  elastische  Energie,  die  frei  wird,  wenn  der  Riss  um 
das  infinitesimale  Flachenelement  d  A  verlangert  wird: 


(4-16) 


d  U  ist  hierbei  die  gespeicherte  elastische  Energie.  Bei  gerade  durch  die  Probe  verlaufende 
Risse  wird  A  wie  folgt  ermittelt: 

A  =  aBp  (4-17) 

wobei  Bp  die  Probendicke  ist  und  a  die  Rissverlangerung  darstellt.  Der  Zusammenhang  zwi¬ 
schen  der  elastischen  Energie  U  und  der  Compliance  C  ergibt  sich  fur  sprode  Materialien 
wie  folgt: 


1  P2C 
U  —  —  PAl  —  — — 

2  2 


mit  der  Compliance  C  (reziproke  Federkonstante) 


C 


M 


(4-18) 


(4-19) 


Die  Bruchenergie  Gc'  ergibt  sich  damit  zu 

2  Bp  da 


(4-20) 


Die  Probennachgiebigkeit  C  kann  aus  dem  Kraftverlangerungsdiagramm  anhand  eines  Stei- 
gungsdreieckes  ermittelt  werden.  Die  zu  dieser  Compliance  zugehorige  Rissverlangerung 
d  a  muss  wurde  wie  bereits  erwahnt  anhand  einer  Kalibrierkurve  bestimmt. 
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4.4.2  Zusammenfassung  des  Druckscherversuches 


Fur  das  karbonisierte  C/C  wurde  die  Bruchenergie  im  Interface  zu  16,3  J/m2  bestimmt  [33]. 
Fur  das  C/C-Material,  dass  bei  einer  Temperatur  von  ca.  1800°C  warmebehandelt  wurde  und 
als  graphitiertes  Material  beschrieben  wird,  wurde  eine  Interfacebruchenergie  von  10,8  J/m2 
bestimmt  [33]. 

Die  so  erhaltenen  Werte  stellen  nicht  die  exakten  Interfacebruchenergien  dar,  da  der  Riss 
sich  zwischen  den  Schichten,  entweder  in  der  Kohlenstoffmatrix  oder  in  der  Faser-/Matrix- 
Grenzflache  fortpflanzt.  Hierbei  ist  allerdings  zu  bemerken,  dass  es  bei  dieser  Berechnung 
nur  auf  das  Rissverhalten  ankommt,  d.h.  ob  der  Matrixriss  entweder  zum  Faserbruch  fuhrt 
oder  beim  Erreichen  der  Faser  abzweigt.  Um  diese  Aussage  zu  treffen,  ist  es  unerheblich,  ob 
sich  der  Riss  exakt  entlang  des  Interfaces  oder  teilweise  noch  in  der  Kohlenstoffmatrix  nach 
dem  Abzweigen  fortbewegt.  Die  ermittelten  Werte  konnen  somit  fur  die  Ermittlung  des  Riss- 
fortschrittsverhaltens  genutzt  werden. 

Die  Verhaltnisse  Gc7Gcf  ermitteln  sich  somit  zu  1,087  fur  den  karbonisierten  kohlenstofffa¬ 
serverstarkten  Kohlenstoff  und  0,72  fur  den  zusatzlich  warmebehandelten.  Daraus  ergeben 
sich  folgende  Bedingungen  fur  ein  Interfaceversagen  in  Abhangigkeit  von  der  Behandlungs- 


temperatur: 


G,  \6,3J/m 2 

>  —> - ! - r- 

G  \5J/m2 


1,087 


fur  karbonisierten  kohlenstofffaserverstarkten  Koh¬ 
lenstoff 


> 


Gj_  >  10,8 Jim  _ 
Gp  15  J/m2 


fiir  zusatzlich  warmebehandelten  (graphitierten)  koh¬ 
lenstofffaserverstarkten  Kohlenstoff 


Die  Verhaltnisse  der  Bruchenergien  haben  sich  infolge  der  Graphitierung  des  Materials  er- 
heblich  verringert.  Dieses  Verhalten  lasst  auf  die  Tendenz  zum  Rissabzweigen  im  zusatzlich 
warmebehandelten  Kohlenstofffaserverstarkten  Kohlenstoff  schlieBen. 


Um  die  in  der  Einleitung  beschriebene  Problemstellung  zu  losen,  werden  im  Anschluss  die 
Verhaltnisse  der  Energiefreisetzungsraten  Gd/Gp  fur  verschieden  herstellungsbedingte  Eigen- 
schaften  des  kohlenstofffaserverstarkten  Kohlenstoff  ermittelt.  Dies  wird  mittels  einer  Fini- 
ten  Element  Methode  erfolgen. 
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4.5  Berechnung  des  Energiefreisetzungsratenverhaltnisses  G,i/G„  fur  kohlenstofffa¬ 

serverstarkten  Kohlenstoff-Einzelfaserverbund 

4.5.1  Allgemeines 

In  Abschnitt  2  wurden  die  theoretischen  Ansatze  beschrieben,  die  mit  Hilfe  der  Energiefrei- 
setzungsrate  den  Rissfortschritt  in  faserverstarkten  Verbundmaterial  untersuchen.  Wie  be- 
reits  ausfiihrlich  diskutiert,  beschreiben  diese  Ansatze  jedoch  immer  nur  eine  oder  wenige 
Eigenschaften,  die  wahrend  eines  Herstellungsprozesses  auftreten.  Um  jedoch  eine  vollstan- 
dige  Analyse  einer  groBen  Anzahl  von  herstellungsbedingten  Einfliissen  auf  einen  Werkstoff 
zu  untersuchen,  bedarf  es  einer  anderen  Methode.  Die  bisherigen  analytischen  Ansatze  lie- 
fern  keine  Aussage  daruber,  wie  sich  herstellungsbedingte  Schrumpfspannungen  in  der  Mat¬ 
rix  auf  den  Rissverlauf  auswirken.  Weiterhin  wurden  bisher  Reibungseffekte  zwischen  Faser 
und  Matrix  vemachlassigt.  Die  teilweise  Graphitierung  der  Matrix  bei  einer  zusatzlichen 
Warmebehandlung  muss  ebenfalls  anhand  numerischen  Losung  angegangen  werden.  Nach- 
dem  nun  das  Finite  Element  Verfahren  anhand  bekannter  analytischen  Losungen  verifiziert 
wurde  kann  nun  die  eigentliche  Erweiterung  zum  kohlenstofffaserverstarkten  Kohlenstoff 
stattfinden.  Dies  Modellierung  wird  im  folgenden  beschrieben. 

4.5.2  Modellierung  des  Einzelfaserverbundes 

Die  bisherigen  Berechnungen  sollen  auf  den  kohlenstofffaserverstarkten  Kohlenstoff  (C/C) 
erweitert  werden.  Dazu  wird  das  bereits  vorgestellte  Modell  fur  einen  Einzelfaserverbund 
mittels  des  Pre-/Postprozessors  PATRAN  verandert.  Ein  Modell  fur  einen  Einzelfaserver¬ 
bund  mit  einem  Faserdurchmesser  von  6,95  pm  und  einen  Gesamtdurchmesser  (mit  Matrix- 
mantel)  von  8,62  pm  wird  im  folgenden  verwendet.  Dies  fuhrt  zu  einem  rechnerischen  Fa- 
servolumengehalt  von  65  %.  Zwischen  Faser  und  Matrix  wird  vorerst  eine  perfekte  Haftung 
angenommen.  Folgende  Materialdaten  [1]  liegen  der  Rechnung  zugrunde. 


K.ohlenstofffaser 

Tenax  HTA  nach  DIN  65382 

Glaskohlenstoffmatrix 

E|| 

225000  MPa 

35000  MPa 

Ei 

15000  MPa 

35000  MPa 

Hz|| 

15000  MPa 

14580  MPa 

V||l=Viz=V2|| 

0,2 

0,2 

«ll 

1,0*1 0'6  1/K 

3,5*  1C6  1/K 

U± 

8,0*1 0'6  1/K 

3,5*  10^  1/K 

1 

'ab.  4-8:  Materialkennwerte  Faser  und  Matrix 
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Es  wird  der  Verlauf  der  Energiefreisetzungsrate  bei  Penetration  des  Matrixrisses  durch  die 
Grenzflache  zwischen  Faser  und  Matrix  ermittelt.  Die  untere  Modellkante  wurde  dabei  wie- 
derum  als  Rissgrund  modelliert,  wobei  der  Riss  anfangs  durch  das  Festhalten  der  Rissufer- 
knoten  in  z-Richtung  (Faserrichtung)  geschlossen  gehalten  wurde.  Der  Riss  wird  durch  ein 
schrittweises  Loschen  der  Verschiebungseinschrankungen  in  Richtung  des  Interfaces  ver- 
groBert.  Im  folgenden  werden  nun  nacheinander  der  Einfluss  von  verschiedenen  herstel- 
lungsbedingten  Eigenschaftsanderungen  auf  das  Bruchverhalten  (Energiefreisetzungsrate) 
ausfuhrlich  untersucht  und  diskutiert.  Es  folgt  zuerst  die  Untersuchung  von  herstellungsbe- 
dingten  Matrixschrumpfspannungen  in  der  Matrix. 

4.6  Ermitteln  Ga/Gn  unter  Einfluss  von  herstellungsbedingten  Matrixschrumpf- 
spannungen  fiir  den  Einzelfaserverbund 

4.6.1  Modellierung  des  Einflusses  von  Matrixschrumpfspannungen 

Zuerst  soil  der  Einfluss  der  herstellungsbedingten  Matrixschrumpfvorgange  (wurde  bereits  in 
der  Einleitung  beschrieben)  auf  den  Rissfortschritt  untersucht  werden.  Messungen  [1]  besta- 
tigen,  dass  nach  der  Karbonisierung  Matrixschrumpfspannungen  mit  einer  unbekannten 
GroBe  vorliegen.  Bei  der  Karbonisierung  wird  die  Polymer-Matrix  in  eine  Glaskohlenstoff- 
matrix  umgewandelt.  Dies  geht  mit  einer  linearen  Schrumpfimg  der  Matrix  von  20%  (Volu- 
menschrumpfung  von  50%)  einher,  soweit  dies  nicht  durch  die  Fasem  verhindert  wird.  Ein 
Teil  dieser  Schrumpfung  findet  spannungsfrei,  ein  Teil  jedoch  durch  Aufbau  von  Schrumpf- 
spannungen  statt.  Das  Schrumpfen  unter  Aufbau  von  Spannungen  wurde  mit  ca.  1-2%  ange- 
nommen.  Im  folgenden  wird  der  Einfluss  von  Matrixschrumpfspannungen  auf  die  Risserwei- 
terung,  ausgelost  durch  eine  Matrixschrumpfung  von  1%  bzw.  2%,  untersucht.  Es  wird  ein 
unmittelbarer  Vergleich  zwischen  einem  C/C-Einzelfaserverbund  mit  0%-  und  1%  bzw.  2%- 
Matrixschrumpfen  notwendig,  um  den  Einfluss  auf  die  Energiefreisetzungsrate  und  damit 
dem  Rissverlauf  zu  bestimmen.  Da  das  Versagen  aufgrund  der  Erfahrungen  vorangegange- 
ner  Versuche  [1]  bei  entsprechcnden  unterschiedlichen  Spannungen  auftritt,  muss  die  Riss- 
erweiterung  in  beiden  Modellen  fur  1%-  bzw.  2%-Matrixschrumpfen  (karbonisiert)  bzw.  0% 
Matrixschrumpfen  (graphitiert)  durch  unterschiedliche  Zugkrafte  belastet  werden. 

Fur  das  unidirektionale  C/C-Biindel  mit  der  Kohlenstofffaser  Tenax  HTA  nach  DIN  65382 
ergibt  sich  bei  einer  Behandlungstemperatur  von  1000°C  nach  der  Karbonisierung  (ent- 
spricht  dem  Fall  mit  Matrixschrumpfspannungen)  eine  Faserbiindelfestigkeit  von  264  MPa 
bezogen  auf  den  Faserquerschmtt  von  0,455  mm  . 
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Da  das  Modell  in  der  Einheit  „Meter“  modelliert  wurde,  ist  es  folglich  in  der  Realitat  um  den 
Faktor  106  vergroBert  worden.  Daraus  ermittelt  sich  eine  Last  von  10,02  kN,  die  als  Last  fur 
den  Einzelfaserverbund  mit  der  Querschnittsflache  von  insgesamt  0,586  mm2  dient. 

Bei  einer  Behandlungstemperatur  von  1800°C  nach  der  Graphitierung  (entspricht  dem  Fall 
ohne  Matrixschrumpfspannungen)  ergibt  sich  eine  Festigkeit  des  C/C-Bundels  von  ca.  888 
MPa.  Dies  resultiert  in  einer  deutlich  hoheren  Last  von  33,6  kN  fur  den  Einzelfaserverbund. 
Auch  hier  muss  wieder  berucksichtigt  werden,  dass  das  Modell  in  der  Einheit  „Meter“  er- 
stellt  ist.  Es  werden  zwei  Berechnungen  durchgefuhrt: 

1)  Modellierung  des  kohlenstofffaserverstarkten  Kohlenstoffs,  der  bis  1800°C  behandelt 
wurde,  wodurch  die  Matrixschrumpfspannungen  relaxiert  sind.  Durch  das  Abkiihlen  auf 
20°C  werden  thermische  Spannungen  infolge  unterschiedlicher  Warmeausdehnungsko- 
effizienten  (Tab.  4-8)  berucksichtigt.  AnschlieBend  erfolgt  die  Belastung  durch  die 
Bruchspannung  von  264  MPa. 

2)  Modellierung  des  kohlenstofffaserverstarkten  Kohlenstoffs,  der  bis  1000°C  behandelt 
wurde.  Dieser  besitzt  Matrixschrumpfspannungen.  Der  Einfluss  dieser  Schrumpfspan- 
nungen  wird  durch  ein  Matrixschrumpfen  von  1%  modelliert.  AnschlieBend  erfolgt  ein 
Abkiihlen  auf  20°C,  um  thermische  Spannungen  aufgrund  unterschiedlicher  Warmeaus- 
dehnungskoeffizienten  (Tab.  4-8)  zu  beriicksichtigen,  und  ein  Zugversuch  mit  einer 
Bruchspannung  von  888  MPa. 

Die  Berechnung  1)  konnte  in  einem  Lastfall  mit  Hilfe  eines  linear  elastischen  Losungsan- 
satzes  (sog.  SOL  101)  des  Finiten  Element  Programms  NASTRAN  durchgefuhrt  werden. 

Um  jedoch  das  Matrixschrumpfen  gemaB  Berechnung  2)  zu  simulieren,  waren  zwei  Lastfal- 
le  notwendig.  Dazu  wurde  den  Faserelementen  und  Matrixelementen  eine  jeweils  unter- 
schiedliche  Anfangstemperatur  zugeordnet.  Um  in  der  Matrix  eine  Schrumpfung  von  1% 
„aufzupragen“  war  bei  einem  Warmeausdehnungskoeffizienten  von  3,5*10 6  1/K  eine  simu- 
lierte  Abkiihlung  von  2857  K  notwendig.  Folglich  erhielt  die  Matrix  die  Ausgangstemperatur 
mit  Hilfe  der  TEMP-Karten  von  3857°C.  Die  Faserelemente  wurden  mit  1000°C  beauf- 
schlagt. 

Dieser  Ausgangszustand  wurde  mit  unterschiedlichen  Temperaturen  fur  Faser  und  Matrix 
durch  TEMPERATURE  (INITIAL)  mit  NASTRAN  bestimmt.  AnschlieBend  erfolgte  im  ers- 
ten  Lastfall  die  simulierte  Abkiihlung  auf  1 000°C,  damit  konnte  die  Matrix  „schrumpfen“. 
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Im  zweiten  Lastfall  wurde  dann  auf  20°C  abgekuhlt.  Diese  Abkiihlung  wurde  dann 
gleichzeitig  in  Faser  und  Matrix  durchgefiihrt. 


4.6.2  Berechnung  der  Energiefreisetzungsrate  Gn  fur  Penetration 


Die  Berechnung  1  und  2  wurde  fur  einen  Rissfortschritt  im  Einzelfaserverbundrand  radial 
nach  innen  durchgefiihrt.  Im  folgenden  wird  nun  beispielhaft  erklart,  wie  die  Energiefreiset- 
zungsraten  aus  der  FEM-Berechnung  ermittelt  werden. 

Bei  jedem  infinitesimalen  Rissfortschritt  wurde  aus  dem  von  NASTRAN  erzeugten  .f06-File 
die  Gesamtenergie  des  Modells  W  gemafi  Abschnitt  2.12.3  ermittelt.  Der  Rissfortschritt  von 
a/aj=0,738  nach  0,816  ergab  ein  W  von  8616,876  J  bzw.  8626,255  J  aus  der  FEM-Rechnung 
(unter  der  Beriicksichtigung,  dass  das  Modell  die  Dimension  „Meter“  besitzt).  Die  Genauig- 
keit  der  Ausgangsgrofien  von  drei  Dezimalstellen  war  notwendig,  um  auch  geringe  Ande- 
rungen  der  Gesamtenergie  (anhand  der  Differenz)  zu  beriicksichtigen. 

Die  Differenz  betragt  9,379  J,  was  der  Energiefreisetzungsrate  fur  die  konzentrische  Risser- 
weiterung  AA=  1,484  m2  entspricht.  Die  Energiefreisetzungsrate  ergibt  sich  nach  Multiplica¬ 
tion  mit  dem  Faktor  2  gemaB  2.12.3  von  12,642  J/m2  fur  den  C/C-Einzelfaserverbund  mit 
1%  Matrixschrumpfen  nach  der  Karbonisierung.  Diese  Werte  sind  fur  beide  Berechnungen 
fur  das  0%-  und  1%  Matrixschrumpfen  in  Tab.  4-9  zusammengefasst. 


C/C  nach  zusatzlicher  Warmebehandlung 

(0%  Schrumpfen) 

C/C  nach  der  Karbonisierung 

(1%  Schrumpfen) 

a/a, [ 1 ] 

G(Matrix)  [J/m2] 

G(Faser)  [J/m2] 

G(Matrix)  [J/m2] 

G(Faser)  [J/m2] 

0,638 

0,949 

11,787 

0,738 

0,992 

12,309 

0,816 

1,019 

12,642 

0,877 

0,924 

1 1,489 

0,924 

0,868 

10,762 

0,960 

0,684 

8,501 

0,988 

0,484 

5,999 

1,013 

2,439 

15,180 

1,040 

4,229 

16,135 

1,069 

5,227 

16,173 

1,100 

6,351 

14,457 

1,132 

7,204 

13,052 

Tab.  4-9:  Energiefreisetzungsraten  fur  Rissfortschritt  in  die  Faser  in  Abhangigkeit 

des  Matrixschrumpfens 
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Die  Energiefreisetzungsraten  mit  der  Einheit  J/m2  werden  in  der  Abb.  4-13  fur  Penetration 
(Risseindringen  in  die  Faser)  fur  die  beiden  Berechnungen  dargestellt.  Anhand  der  Abb.4-13 
sieht  man  eine  deutliche  Erhohung  der  Energiefreisetzungsrate  fur  Penetration  Gp  fur  das  un- 
relaxierte  Material.  Die  zusatzliche  Warmebehandlung  oberhalb  von  1000°C  bewirkt  eine 
deutliche  Verringerung  der  Energiefreisetzungsrate  fur  Penetration.  Dies  bedeutet  eine  Rela¬ 
xation  und  damit  Erholung  des  Materials 


O 
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Abb.  4-13:  Energiefreisetzungsraten  eines  penetrierenden  Matrixrisses  in  C/C  fur  0%-  und 

1%  Schrumpfung 

Die  Energiefreisetzungsraten  nehmen  beim  Erreichen  des  Interfaces  ab.  Dies  ist  damit  zu  er- 
klaren,  dass  G=l/E  entspricht.  Mit  dem  Ubergang  in  ein  Medium  hoherer  Steifigkeit  (Koh- 
lenstofffaser)  erfolgt  somit  eine  rapide  Abnahme  der  Energiefreisetzungsrate. 

4.6.3  Verhaltnis  der  Energiefreisetzungsraten  fur  Penetration 

Um  weitere  Aussagen  uber  das  Rissfortschrittsverhalten  bei  aufgebrachter  Belastung  zu  er- 
halten,  ist  es  notwendig,  die  errechneten  Energiefreisetzungsraten  far  die  beiden  unterschied- 
lich  behandelten  Materialien  ins  Verhaltnis  zu  setzen. 
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Abb.  4-14  zeigt,  dass  das  Verhaltnis  der  Energiefreisetzungsraten  G(1%)/G(0%)  in  der  Mat¬ 
rix  bei  aufgebrachter  Last  fur  die  beiden  Falle  jeweils  mit  einem  Faktor  von  ca.  12  nahezu 
konstant  ist.  Erst  nach  Erreichen  des  Interfaces  fallt  diese  in  der  Faser  bei  weiterem  Rissfort- 
schritt  ab. 

Die  Energiefreisetzungsrate  ist  proportional  zum  Quadrat  der  auftretenden  Gesamtspannung. 
Diese  Gesamtspannung  in  der  Matrix  setzt  sich  zusammen  aus  einem  Schrumpfanteil  und  ei¬ 
nem  Lastanteil  zusammen: 

G  -O2  (4-21) 

®  ~  ® Schnimpfen  +  ®  Last  (4*22) 


Der  Betrag  von  aSchrumpfen  tritt  nur  in  der  Matrix  auf,  erst  beim  Erreichen  des  Interfaces  wird 
dieser  0.  Damit  ergeben  sich  im  Matrixbereich  hohe  Betrage  fur  die  Energiefreisetzungsrate 
Gp.  Diese  bewirken  beim  karbonisiertem  Material  Matrixversagen.  Im  graphitierten  Material 
sind  die  Schrumpfanteile  in  Gp  nicht  vorhanden,  wodurch  ein  Matrixversagen  unwahrschein- 
lich  wird. 

4.6.4  Vergleich  von  Penetration  und  Deflection 

Fur  die  Berechnung  mit  1%-  und  0%-Matrixschrumpfen  muss  das  Verhaltnis  Gd/Gp  ermittelt 
werden  und  anschlieBend  anhand  der  bereits  nach  Gupta  und  Martinez  [11,  22]  bzw.  He  und 
Hutchinson  [5]  berechneten  Werte  verglichen  werden.  Fur  den  Fall  der  Rissabzweigung  ins 
Interface  wird  der  Riss  nach  Beriihren  der  Faser  urn  90°  abgelenkt.  In  der  FEM-Berechnung 
erfolgt  dies  durch  das  aufeinanderfolgende  Losen  der  Knoten  im  Interface,  wobei  die  or¬ 
thotrope  Faser  von  der  isotropen  Matrix  gelost  wird.  Die  Energiefreisetzungsraten  ermitteln 
sich  nach  der  bereits  erwahnten  Methode.  Die  Ergebnisse  fur  die  Abzweigung  sind  in  Tab.  4- 
1 1  unter  Gd  zusammengefasst.  AnschlieBend  werden  noch  die  charakteristischen  Verhaltnis- 
se  Gd/Gp  fur  den  Fall  mit  1%-Matrixschrumpfen  und  0%-Matrixschrumpfen  ermittelt. 
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Gp[J/m5] 

G„[J/mJ] 

a/ai 

0%  Matrixschrumpfen 

0=880  MPa 

1%  Matrixschrumpfen 

0=264  MPa 

0%  Matrixschrumpfen 

0=880  MPa 

1%  Matrixschrumpfen 

0=264  MPa 

1,013 

2,439 

15,180 

1,040 

4,229 

16,135 

1,069 

5,227 

16,173 

1,100 

6,351 

14,457 

1,132 

7,204 

13,052 

1,015 

0,213 

1,360 

1,044 

0,303 

2,480 

1,075 

0,335 

2,491 

1,106 

0,383 

2,617 

Tab.4-1 1 :  Energiefreisetzungsrate  fur  Rissabzweigung  und  Risspenetration  am  Interface 


G  bzw.Gd  [J/m2] 


C/C  karbonisiert 


C/C  graphitiert 


a/a,[l] 


Einzelfaserverbund 

C/C  karbonisiert 

0=264  MPa 

C/C  graphitiert 

0=888  MPa 

Gd 

1,36  J/m2 

0,213  J/m2 

Gr 

15,18  J/m2 

2,439  J/m2 

Gd/Gp 

0,0896 

0,0875 

Abb.  4-15:  Energiefreisetzungsrate  G  fur  Rissabzweigung/Risspenetration 
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4.6.5  Ermitteln  von  Gd/Gn 

Um  Aussagen  iiber  den  Verlauf  eines  Matrixrisses  zu  treffen,  muss  das  folgende  Verhaltnis 


G^<G± 

G{  Gp 


(4-23) 


betrachtet  werden,  bei  dem  Rissabzweigen  auftritt.  Die  linke  Seite  Gc'/Gcf  gibt  das  Verhaltnis 
der  Interfacebruchenergie  Gc'  zu  der  Bruchenergie  der  Faser  Gcf  wieder.  Das  Verhaltnis  der 
Energiefreisetzungsraten  Gd/Gp  ist  bereits  fur  Faserverbunde  mit  isotropen  Werkstoffkompo- 
nenten  gcmaB  He  und  Hutchinson  [5]  bekannt  und  fur  orthotrope  Materialien  aus  Gupta  und 
Martinez  [1 1,  22]  bekannt.  Diese  wurden  fur  den  kohlenstofffaserverstarkten  Kohlenstoff  be¬ 
reits  aus  den  analytischen  Ansatzen  bestimmt.  Allerdings  vemachlassigen  die  bisherigen 
Theorien  den  Einfluss  der  unterschiedlichen  Warmeausdehnungskoeffizienten  und  den  Ein- 
fluss  von  Matrixspannungen  aufgrund  der  Volumenschrumpfung.  Aus  den  Tabellenwerten 
nach  Tab.  4-1 1  kann  nun  das  Verhaltnis  Gd/Gp  unter  Beriicksichtigung  des  Matrixschrump- 
fens  und  damit  der  unterschiedlichen  Warmebehandlung  ermittelt  werden: 


Gd/Gp 

fur  C/C  nach  zusatzlicher 

Warmebehandlung 

(0%-Matrixschmmpfen) 

Gd/Gp 

fur  C/C  nach  der 

Karbonisierung 

( 1  %-Matrixschrumpfen) 

0,0875 

0,0896 

Tab.4-12:  Verhaltnis  der  Energiefreisetzungsraten  Gd/Gp  fur  0%-  und  1%-Matrixschrumpfen 


Die  Werte  in  Tab.4-12  zeigen,  dass  die  zusatzliche  Warmebehandlung  nur  eine  geringe  Ver- 
anderung  des  Energiefreisetzungsratenverhaltnis  beim  Einzelfaserverbund  bewirkt. 
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4.6.6  Zusammenfassung  der  Berechnungen  der  Energiefreisetzungsraten  Gd  bzw.  Gn 
des  Einzelfaserverbundes 


Ergebnisse  fur  Gr 

Abb.  4-15  zeigt,  dass  die  Energiefreisetzungsraten  bei  Rissabzweigen  fur  das  unrelaxierte, 
schrumpfspannungsbehaftete  Material  nach  der  Karbonisierung  (1%  Matrixschrumpfen,  o  = 
264  MPa  )  trotz  kleinerer  Last  groBer  sind  als  fur  das  zusatzlich  warmebehandelte  Material 
(0%  Matrixschrumpfen  o  =  888  MPa). 

Die  Energiefreisetzungsrate  Gp  des  karbonisierten  Materials  (1%  Matrixschrumpfen)  in  der 
Faser  ist  ebenfalls  groBer  als  die  im  zusatzlich  warmebehandelten,  graphitierten  Material. 
Ein  Faserbruch  wird  somit  im  karbonisierten  Material  bei  gleicher  Bruchenergie  der  Faser 
begunstigt.  Die  Bruchenergie  der  Faser  Gcf  wurde  gemaB  Tab. 4-5  und  4-6  mit  15  J/m2  ermit- 
telt.  Der  Wert  der  Energiefreisetzungsrate  Gp  iiberschreitet  mit  15,18  J/m2  gemaB  Tab.  4-9 
fur  das  karbonisierte  Material  die  Bruchenergie  von  15/m2.  Wohingegen  die  ermittelte  Ener¬ 
giefreisetzungsrate  fur  den  zusatzlich  warmebehandelten  Werkstoff  mit  2,439  J/m2  weit  un- 
terhalb  der  Bruchenergie  liegt. 

GemaB  Ahlbom,  Chou  und  Okura  [41]  wird  die  kritische  Bruchenergie  fur  eine  Glaskohlen- 
stoffmatrix  nach  einer  Behandlung  bei  600°C  mit  16  J/m2  und  bei  2000°C  mit  37  J/irf  ange- 
geben.  Werte  gemaB  [42,  43]  liegen  etwas  hoher.  Bei  der  Berechung  mit  0%- 
Matrixschrumpfen  und  1%  Matrixschrumpfen  werden  diese  Bruchenergien  nicht  erreicht.  Es 
zeigt  sich  aber  deutlich,  dass  mit  zunehmender  Schrumpfspannung  in  der  Matrix  die  Ener¬ 
giefreisetzungsrate  stark  anwachst  und  bei  einem  bestimmten  Wert  die  Bruchenergie  der 
Matrix  iiberschreitet. 

Auffallend  ist  auch,  dass  die  Energiefreisetzungsrate  fur  1%-Matrixschrumpfen  nach  Passie- 
ren  des  Interfaces  und  Eindringen  in  die  Faser  nach  kurzem,  annahemd  konstanten  Verlauf 
abfallt.  Im  Gegensatz  dazu  steigt  G  bei  dem  zusatzlich  warmebehandelten  Material  (0%  Mat¬ 
rixschrumpfen)  an. 
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Bereits  in  Arbeiten  von  Eric  Martin  [17,  44]  zeigt  sich,  dass  bei  der  Beanspruchung  von  Ein- 
zelfaserverbunden  mit  geringen  Lasten  ein  groBer  Einfluss  auf  den  Verlauf  der  Energiefrei- 
setzungsrate  durch  die  thermischen  Spannungen  als  Folge  der  unterschiedlichen  Warmeaus- 
dehnungskoeffizienten  gegeben  wird.  Der  Einfluss  auf  die  Energiefreisetzungsrate  zeigt  sich 
bei  den  von  mir  durchgefuhrten  Berechnungen  im  Fall  ohne  Matrixrestspannungen  nicht,  da 
hier  eine  wesentlich  groBere  Last  aufgepragt  wurde  und  der  Einfluss  der  Warmeausdeh- 
nungsverhaltnisse  somit  in  den  Hintergrund  gerat. 

Weiterhin  ist  der  Verlauf  in  der  Nahe  des  Interface  von  groBer  Bedeutung.  Wie  bereits  von 
den  analytischen  Ansatzen  von  He  und  Hutchinson  [5]  bzw.  Gupta  und  Martinez  [1 1,  22]  aus 
Abschnitt  2  bekannt,  besteht  am  Interfacebereich  eine  Singularitat.  Die  Gleichung  2-14  fur 
0>X  <0,5  und  2-67  fur  0>y>-0,5  sagen  hier  fur  Aa->0  einen  Verlauf  fur  Gp  gegen  0  voraus. 
Der  Grenzwert  fur  Gp  wird  unabhangig  a,  wenn  X  und  y  den  exakten  Wert  von  'A  annehmen. 
Fur  diesen  Sonderfall  existiert  ein  exakter  Grenzwert  fur  Gp.  Aus  Gleichung  2-12  ergibt  sich 
fur  den  in  dieser  Arbeit  untersuchten  kohlenstofffaserverstarkten  Kohlenstoff,  dass  X  nur 
wenig  von  'A  abweicht.  Daher  ist  bei  einer  Annaherung  mit  Aa— >0  im  Interface  ebenfalls  ein 
Wert  Gp  zu  ermitteln. 

Weiterhin  sagen  beide  Theorien  aus,  dass  es  fur  das  Verhaltnis  Gd/Gp  einen  Grenzwert  gibt, 
bzw.  fur  den  Sonderfall,  wenn  der  Dundursparameter  a=0  (Fasermaterial  =  Matrixmaterial) 
ist.  Das  Verhaltnis  Gd/Gp  ist  namlich  gemaB  Gleichung  2-20  unabhangig  von  a. 

Ergebnisse  fur  Gd/Gp 

Die  Energiefreisetzungsratenverhaltnisse  Gd/Gp  betragen  0,0896  nach  der  Karbonisierung  bei 
1000°C  (1%  Matrixschrumpfen)  und  0,0875  nach  der  zusatzlichen  Warmebehandlung  (0% 
Matrixschrumpfen).  Aus  dem  Vergleich  Gd/Gp  fur  0%-  und  1%-Matrixschrumpfen  erkennt 
man,  dass  die  zusatzliche  Warmebehandlung  des  karbonisierten  Materials  und  die  damit  ein- 
hergehende  Relaxation  der  Matrixschrumpfspannungen  einen  (wenn  auch  geringen)  Einfluss 
zeigt.  Beide  Verhaltnisse  sind  im  Vergleich  mit  den  von  Gupta  mit  0,1 18  und  He  und  Hut¬ 
chinson  mit  0,248  bestimmten  Werte  sehr  niedrig. 

Offenbar  spielt  hierbei  der  Unterschied  zwischen  Einzelfaserverbund  und  Quasimehrfaser- 
verbund  eine  wesentliche  Rolle.  Um  die  Wirklichkeit  besser  nachzuahmen,  ist  daher  eine 
Erweiterung  zum  Mehrfaserverbund  notwendig.  Aus  diesem  Grande  wird  in  Abschnitt  4.7 
eine  Randbedingungen  zusatzlich  eingefuhrt,  die  eine  Randbedingung  im  Quasimehrfaser- 
verbund  darstellt. 
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Zusatzlich  sei  nochmals  auf  den  theoretischen  Ansatz  von  Petrovic  [9]  hingewiesen.  Im  Ein- 
zelfaserverbund  ist  gemafi  Petrovic  [9]  keine  Abhangigkeit  von  thermischen  Lasten  auf  das 
Energiefreisetzungsratenverhaltnis  Gd/Gp  zu  erwarten,  wenn  die  mechanische  Last  und  die 
thermische  Last  sich  in  Phase  befinden  und  eine  perfekte  Haftung  im  Interfacebereich  vor- 
liegt. 

4.6.7  Ergebnisse  fur  Spannungsverteilungen  am  Interface  im  Einzelfaserverbund 

Die  Ergebnisse  aus  den  Spannungsanalysen  werden  nun  erganzend  zu  der  Analyse  der  Ener- 
giefreisetzungsrate  dargestellt.  Hierbei  wird  der  Verlauf  von  Axialspannung  in  der  Faser, 
Radialspannung  und  Schubspannung  entlang  des  Interfaces  vom  Matrixriss  aus  bis  an  den 
Modellrand  betrachtet. 

Axialspannungen  in  der  Faser  im  Einzelfaserverbund 

Die  Darstellung  der  Axialspannungen  in  der  Faser  gemaB  Abb.  4-16  wurde  fur  0%- 
Matrixschrumpfen  (relaxiertes  Material),  1  %-Matrixschrumpfen  (  karbonisiertes  Material  ) 
und  2%-Matrix-schrumpfen  erstellt.  Die  Kurven  fur  2%-Matrixschrumpfen  dienen  als  Ver- 
gleich,  um  den  Einfluss  groBer  Matrixschrumpfspannungen  zu  verdeutlichen. 


Abb.4-16:  Axialspannungen  (Faser)  im  Einzelfaserverbund 
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Aufgrund  der  Annahme  der  perfekten  Haftung  tritt  in  Rissnahe  eine  Mode-I  Spannungssin- 
gularitat  fur  alle  drei  Falle  auf.  Der  Riss  ist  dort  folglich  einer  Zugbelastung  ausgesetzt.  Der 
Spannungsverlauf  andert  sich  im  Interface  iiber  den  betrachteten  Bereich.  Langs  des  Inter¬ 
faces  fallen  die  Axialspannungen  fur  1%-  und  2%-Matrixschrumpfen  in  den  Druckbereich. 
Nach  einem  Druckmaximum  zwischen  1,5  und  3  jim  Interfacelange  verringem  sich  diese 
Axialspannungen  leicht.  Sie  gehen  aber  nicht  in  den  Zugbereich  iiber.  Der  Riickgang  der 
1%-Kurve  auf  „Null“  ist  zufallig. 

Die  Axialspannungen  fur  das  relaxierte  Material  bleiben  aufgrund  der  hohen  mechanischen 
Beanspruchung  iiber  das  ganze  Interface  im  Zugbereich. 

Radialspannungen  im  Einzelfaserverbund 

Die  Radialspannungen  fur  0%-,  1%-  und  2%-Matrixschrumpfen  sind  in  Abb.  4-17  darge- 
stellt. 


Abb.4-17:  Radialspannungen  im  Einzelfaserverbund 

Charakteristisch  fur  den  Verlauf  fur  alle  drei  Kurven  ist  wiederum  die  Spannungssingularitat 
an  der  Rissspitze,  die  iiber  den  weiteren  Verlauf  des  Interfaces  abklingt.  Die  Singularity 
wird  auch  hier  durch  die  Annahme  der  perfekten  Haftung  zwischen  Faser  und  Matrix  im 
C/C-Einzelfaserverbund  verursacht. 
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Auffallend  ist  ebenfalls  der  charakteristische  Unterschied  zwischen  dem  relaxierten  Material 
und  dem  unrelaxierten  (karbonisierten)  Material,  indem  Matrixschrumpfung  nach  der  Kar- 
bonisierung  aufgetreten  ist.  Die  an  der  Rissspitze  als  Zugbeanspruchung  ausgepragte  Span- 
nungssingularitat  fallt  im  weiteren  Verlauf  des  Interfaces  in  den  Druckbereich  ab.  Das 
Druckmaximum  befmdet  sich  zwischen  0,75  und  2,25  |im  Interfacelange  von  der  Rissspitze 
entfemt. 

Im  Gegensatz  dazu  unterliegt  der  relaxierte  C/C-Einzelfaserverbund  iiber  den  gesamten  In- 
terfacebereich  einer  Zugbeanspruchung  zwischen  Faser  und  Matrix,  was  zu  einer  Begunsti- 
gung  von  Interfacerissen  fuhrt.  Die  Erklarung  hierfiir  liegt  in  den  unterschiedlichen  War- 
meausdehungskoeffizienten  in  Faser  und  Matrix  bei  identischer  Querkontraktionszahl.  Beim 
relaxierten  Material  kontrahiert  die  Faser  radial  starker  als  die  Matrix  beim  Abkuhlen.  Beim 
unrelaxierten  Material  kompensieren  die  Matrixschrumpfspannungen  den  Effekt  der  unter¬ 
schiedlichen  Warmeausdehnungskoeffizienten.  Das  Vorhandensein  von  Druckspannungen 
iiber  den  Interfacebereich  bei  dem  karbonisierten  C/C  reduziert  die  Anfalligkeit  zur  Bildung 
von  Interfacerissen. 

Schubspannungen  im  Einzelfaserverbund  in  der  Grenzflache 

Betrachtet  man  die  Schubspannungsverlaufe  in  Abb.  4-18,  so  erkennt  man  ebenfalls  einen 
bedeutenden  Einfluss  der  Matrixschrumpfung. 


Abb.4-18:  Schubspannungen  im  Einzelfaserverbund 
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Ein  Schrumpfen  der  Matrix  um  1%  und  2%  verursacht  groBe  Schubspannungen  im  Inter- 
facebereich.  Die  Kurve  mit  0%-Schrumpfen  fallt  deutlich  geringer  aus.  Die  Verlaufe  sind 
aber  in  ihrer  Charakteristik  gleich.  Alle  besitzen  eine  Singularity  an  der  Rissspitze.  Der 
Schubspannungsverlauf  steigt  bereits  kurz  nach  der  Rissspitze  im  Interface  auf  eine  Maxi¬ 
mum,  was  allerdings  bei  dem  relaxierten  Material  sehr  gering  ausfallt.  Nach  dem  Maximum 
fallen  die  Schubspannungen  zum  oberen  Modellrand  hin  und  damit  mit  zunehmender  Inter- 
facelange  auf  den  Wert  0. 
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4.7  Berechnung  des  Energiefreisetzungsratenverhaltnisses  Ga/Gn  fur  kohlenstofffa- 

serverstarkten  Kohlenstoff-Ouasimehrfaserverbund 

4.7.1  Allgemeines 

Die  bisherigen  Berechnungen  wurden  ausschliefllich  fur  einen  Einzelfaserverbund  ermittelt. 
Die  Versuche  wurden  allerdings  an  Mehrfaserverbunden  durchgefuhrt.  Es  ist  somit  eine  Er- 
weiterung  des  Modells  notwendig.  Im  folgenden  wurde  eine  zusatzliche  Randbedingung  ein- 
gefuhrt.  Diese  Randbedingung  stellt  eine  der  Randbedingungen  eines  Mehrfaserverbundes 
dar.  Versagt  ein  Einzelfaserverbund,  der  in  einem  Mehrfaserverbund  eingebettet  ist,  dann 
geschieht  dies  unter  anderen  Randbedingungen  als  beim  Versagen  eines  einzelnen  Einzelfa- 
serverbundes.  Beim  Versagen  des  Einzelfaserverbundes  in  einem  Mehrfaserverbund  kann 
davon  ausgegangen  werden,  dass  die  radiale  Verschiebung  auf  der  MatrixauBenflache  in 
Rissnahe  identisch  mit  der  Verschiebung  in  einer  Entfemung  vom  Riss  ist.  Diese  konstante 
radiale  Verformung  wird  nun  in  einer  Berechnung  aufgezwungen.  Diese  behinderte  Radial- 
verschiebung  wird  im  folgenden  als  Quasimehrfaserverbund  bezeichnet. 


Abb.4-19:  Quasimehrfaserverbund  mit  aufgepragter  konstanter  Radialverschiebung 
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4.7.2  Modellierung  eines  Quasimehrfaserverbundesf  Behinderung  der  radialen  Ver- 
schiebung) 

Die  radiale  Aufweitung,  die  im  Einzelfaserverbund  in  der  Nahe  des  Risses  groBer  als  im  un- 
gestorten  Bereich  ist,  muss  verhindert  werden.  Dazu  war  es  notwendig,  eine  Berechnung  oh- 
ne  aufgepragte  Radialverschiebungen  durchzufuhren  und  die  radiale  Verschiebung  im  unge- 
storten  Bereich  zu  ermitteln.  Die  Verschiebung  muss  gemaB  Abb.4-19  radial  fur  den  gesam- 
ten  Einzelfaserverbund  (mit  konstanter  GroBe)  aufgepragt  werden.  Dadurch  wird  eine  Rand- 
bedingung  ahnlich  wie  in  einem  Mehrfaserverbund  erzeugt  (im  folgenden  Quasimehrfa- 
serverbund  genannt).  Fur  die  u.a.  Falle  wurden  die  entsprechenden  Berechnungen  durchge- 
fuhrt  und  wiederum  die  Energiefreisetzungsraten  fur  Penetration  und  Rissabzweigung  be- 
rechnet  und  ins  Verhaltnis  gesetzt. 

a)  kohlenstofffaserverstarkter  Kohlenstoff  nach  der  Karbonisierung  (1%  Matrixschrump- 
fen) 

b)  kohlenstofffaserverstarkter  Kohlenstoff  nach  einer  zusatzlichen  Warmebehandlung  (0% 
Matrixschrumpfen) 

4.7.3  Energiefreisetzungsrate  fur  den  schrumpfspannungsbehafteten  kohlenstofffa- 
serverstarkten  Kohlenstoff  nach  der  Karbonisierung  bzw.  Graphitierung 

Die  Tab.  4-13  stellt  die  Energiefreisetzungsraten  Gd  und  Gp  bei  der  Berechnung  mit  Matrix- 
schrampfen  von  1%  dar. 


a/ai  [ 1  ] 

Gp  [J/m2]  (1%)  c=264MPa 

Gd  [J/m2]  (1%)  cr=264MPa 

1,013 

13,179 

1,015 

1,906 

1,044 

2,255 

1,075 

2,204 

1,106 

2,237 

Tab.  4-13  Energiefreisetzungsraten  des  Quasimehrfaserverbundes  nach  der  Karbonisierung 

(1%  Matrixschrumpfen) 


Die  Tab.  4-14  zeigt  die  Werte  fur  den  Fall  0%-Matrixschrumpfspannungen. 
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a/ai  [  1  ] 

Gp  [J/m2]  (0%)  0=888  MPa 

Gd  [J/m2]  (0%)  0=888  MPa 

1,013 

2,514 

1,015 

0,318 

1,044 

0,499 

1,075 

0,572 

1,106 

0,662 

Tab.4-14:  Energiefreisetzungsraten  des  Quasimehrfaserverbundes  nach  der  zusatzlichen 


Warmebehandlung  (Graphitierung,  0%  Matrixschrumpfen) 


Gcl 


Zusatzliche 

Randbedingung 

C/C  karbonisiert 
0=264  MPa 

C/C  graphitiert 
0=888  MPa 

Gd 

1,906  J/m2 

0,318  J/m2 

13,179  J/m2 

2,514  J/m2 

Gd/Gp 

0,145 

0,126 

Abb.4-20:  Energiefreisetzungsraten  des  Quasimehrfaserverbundes  fur  1  %-Matrixschrumpfen 
nach  der  Karbonisierung  und  fur  0%  Matrixschrumpfen  nach  der  Graphitierung 

Nach  der  Karbonisierung  ist  eine  Restspannung  in  der  Matrix  vorhanden  (Abb.  4-20).  Diese 
entsteht  aufgrund  von  1%  Matrixschrumpfen.  Diese  Schrumpfspannungen  konnen  erst  durch 
eine  anschlieBende  zusatzliche  Temperaturbehandlung  relaxiert  werden.  Aus  der  Berechnung 
unter  Berucksichtigung  von  1%  Matrixschrumpfen  erhalt  man  fur  die  Energiefreisetzungsra- 
te  Gp  fur  einen  moglichen  Faserbruch  einen  Wert  von  13,179  J/m2.  Fur  die  Energiefreiset- 
zungsrate  fur  Rissabzweigung  am  Interface  wird  ein  Wert  von  Gd=l,906  J/m2  ermittelt. 
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Fur  das  graphitierte  Material  sind  die  Energiefreisetzungsraten  nach  der  zusatzlichen  War- 
mebehandlung  ebenfalls  in  Abb.  4-20  dargestellt.  Hierbei  betragt  die  Matrixschrumpfspan- 
nung  0%.  Lediglich  thermische  Spannungen  aufgrund  unterschiedlicher  Warmeausdeh- 
nungskoeffizienten  spielen  hier  eine  Rolle.  Die  Finite  Element  Rechnung  ergibt  fur  die  Riss- 
abzweigung  eine  Energiefreisetzungsrate  Gd  von  0,318  J/m2.  Die  Energiefreisetzungsrate  Gp, 
die  einen  Faserbruch  erzeugt,  wird  mit  2,514  J/m2  ermittelt. 

4.7.4  Zusammenfassung  der  Berechnung  des  Quasimehrfaserverbundes 

Im  Gegensatz  zu  den  Rechnungen  des  Einzelfaserverbundes,  die  bereits  erstellt  worden  sind, 
wird  durch  die  zusatzlichen  Randbedingung  ein  Quasimehrfaserverbund  realisiert. 

Aus  Abb.  4-20  wird  der  Gd/Gp-Wert  von  0,318/2,514=0,1265  fur  das  relaxierte  Material  und 
ein  Gd/Gp-Wert  von  1,906/13,179=0,1446  far  das  karbonisierte  Material  ermittelt.  Die  Aus- 
wertung  zeigt,  dass  sich  die  Verhaltniswerte  (Deflection/Penetration)  um  ca.  14,31%  von 
0,1265  auf  0,1446  vergrdBert,  wenn  Matrixspannungen  vorhanden  sind.  Diese  Erhohung  ist 
deutlicher  als  bei  den  bisherigen  Modellen. 

Nach  der  Karbonisierung  betragt  die  Energiefreisetzungsrate  fur  Gp  13,179  J/m2.  Die  zusatz- 
liche  Warmebehandlung  lasst  die  Energiefreisetzungsrate  deutlich  sinken  (Gp  =  2,514  J/m2). 
Der  Wert  von  Gp  =  2,514  J/m2  ist  zu  gering  um  eine  Schadigung  der  Faser  nach  der  Graphi- 
tierung  zu  erzeugen.  Die  deutlich  hohere  Energiefreisetzungsrate  von  Gp=13,179  J/m2  er- 
scheint  fur  die  Kohlenstofffaser  nach  der  Karbonisierung  deutlich  kritischer  zu  sein.  Die  E- 
nergiefreisetzungsraten  in  der  Matrix  sind  in  dem  schrumpfspannungsbehafteten  Material 
ebenfalls  sehr  groB,  da  G^o2  ist.  Daher  ist  ein  Versagen  der  Faser  nach  der  Karbonisierung 
infolge  von  Matrixrissen  sehr  wahrscheinlich.  Eine  Anderung  im  Bruchmodus  in  Folge  der 
Graphitierung  ist  zu  erwarten,  da  weniger  Matrixrisse  auftreten  und  die  Energiefreisetzungs¬ 
rate  fur  Penetration  relativ  gering  ist.  Hierbei  sei  jedoch  kritisch  zu  bemerken,  dass  im  Ein- 
zelfaserverbund  eine  Risspenetration  zum  Versagen  des  Einzelfaserverbundes  fuhrt.  Dage- 
gen  muss  der  Mehrfaserverbund  noch  nicht  versagen,  wenn  ein  Einzelfaserverbund  als  Teil 
des  Gesamtverbundes  versagt. 

Im  folgenden  Abschnitt  wird  nun  nicht  mehr  eine  perfekte  Haftung  zwischen  Faser  und  Mat¬ 
rix  angenommen.  Die  Kraftiibertragung  erfolgt  durch  Reibung  im  Interface. 
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4.8  Untersuchung  von  Reibungseinfliissen  im  Interface 

4.8.1  AHgemeines 

In  den  bisherigen  Untersuchungen  wurde  zwischen  Faser  und  Matrix  eine  perfekte  Haftung 
vorausgesetzt.  Dadurch  ist  eine  Spannungsubertragung  in  der  Faser-Matrix-Grenzflache  ge- 
wahrleistet.  Krafte  zwischen  den  beiden  Medien  korinen  jedoch  auch  durch  Reibung  iibertra- 
gen  werden.  Die  auf  die  Faser  einwirkende  radiale  Druckspannung  aus  dem  Abkiihlprozess 
ermoglicht  iiber  das  Coulombsche  Gesetz  Reibungskrafte  zu  iibertragen,  die  das  Verhalten 
des  Einzelfaserverbundes  erheblich  beeinflussen.  Aus  diesem  Grunde  wurde  ein  Modell 
entwickelt,  mit  dem  Reibung  im  Interface  modelliert  werden  konnte.  NASTRAN  bietet  hier 
die  Moglichkeit,  dies  in  Form  von  GAP-Elementen  [45]  zu  modellieren.  Mit  Hilfe  dieser 
GAP-Elemente  ist  es  moglich,  zwischen  Faser  und  Matrix  die  Kraftubertragung  durch  Rei¬ 
bung  zuzulassen.  Hierbei  soil  der  Unterschied  zur  perfekten  Haftung  zwischen  Faser  und 
Matrix  untersucht  werden.  Als  erstes  werden  die  Bereiche  ermittelt,  in  denen  Kraftubertra¬ 
gung  zwischen  Faser  und  Matrix  aufgrund  von  Haftreibung  bzw.  Gleitreibung  erfolgt.  Dabei 
wird  zwischen  dem  Einzelfaserverbund  und  dem  Quasimehrfaserverbund  unterschieden, 
wobei  die  GroBe  der  Matrixschrumpfspannungen  den  wesentlichen  Einflussfaktor  darstellt. 

4.8.2  Modellierung  mit  Hilfe  des  FEM-Programms  NASTRAN 

Fur  die  folgenden  Berechnungen  war  eine  Verdoppelung  der  Modellange  notwendig,  da 
Storspannungen  am  Modellende  ausklingen  sollen.  Bei  der  Modellierung  der  Reibung  (ge- 
ringere  Spannungsubertragung  in  der  Faser/Matrix-Grenzflache)  erfordert  dies  eine  groBere 
Modelllange.  Die  Seitenlange  (Hohe)  im  Interfacebereich  wurde  gleich  gewahlt,  um  eine 
gleichmaBige  Vemetzung  zu  erhalten.  Dadurch  ist  es  einfacher  die  Gleitlange  im  Interface¬ 
bereich  zu  bestimmen.  Je  nach  GroBe  der  auftretenden  Schubspannungen  im  Verhaltnis  zur 
Coulombschen  Reibung  kommt  es  in  der  Grenzflache  zwischen  Faser  und  Matrix  entweder 
zur  Ausbildung  von  einer  Gleitreibungszone  oder  einer  Ilaftreibungszone. 

Fur  die  Gleitreibungszone  (Gleitlange)  gilt: 

o rz  ^  ^r,  (4’24) 

Fiir  die  Haftreibungszone  gilt: 

°n  <  V-Orr  (4‘25) 
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Es  sollen  im  folgenden  zwei  grundsatzliche  Probleme  in  Bezug  auf  Reibung  zwischen  Faser 
und  Matrix  [46,  47, 48,  49]  behandelt  werden: 

•  Einfluss  des  Reibungskoeffizienten  zwischen  Faser  und  Matrix  auf  die  Gleitlange 

•  Einfluss  des  Reibungskoeffizienten  zwischen  Faser  und  Matrix  auf  Gp 

4.8.3  Einfluss  der  Reibung  auf  die  Gleitlange 

Flierzu  wurde  der  Reibungskoeffizient  p  in  einem  Bereich  von  0,2  und  0,8  variiert.  Anhand 
des  ,f06-Files,  das  von  NASTRAN  erzeugt  wird,  kann  die  Gleitlange  ermittelt  werden.  Fur 
jeden  Reibungskoeffizienten  wurde  eine  separate  Rechnung  durchgefiihrt. 

4.8.4  Einfluss  der  Reibung  auf  die  Gleitlange  bei  1%  und  2%-  Matrixschrumpfen 
beim  Einzelfaserverbund 


Fur  den  Lastfall  mit  0=264  MPa  und  1%  Matrixschrumpfspannungen  wurde  die  Gleitlange 
in  Abhangigkeit  vom  Reibungskoeffizienten  p  fur  den  Einzelfaserverbund  und  den  Quasi- 
mehrfaserverbund  [mit  behinderten  radialen  Verschiebungen]  ermittelt.  Die  Gleitlangen 
werden  graphisch  und  tabellarisch  dargestellt.  Es  ist  zu  erwarten,  dass  herstellungsbedingte 
Einflusse,  die  eine  groflere  Matrixschrumpfspannung  hervorrufen,  eine  deutliche  Anderung 
in  den  Gleitlangen  verursachen.  Dieser  Einfluss  wird  mit  einer  Vergleichsrechnung  ermittelt. 
Dazu  wurden  die  Matrixschrumpfspannungen  verdoppelt  (2%  Matrixschrumpfen)  und  wie- 
derum  Berechnungen  fur  den  Einzelfaserverbund  und  den  Quasimehrfaserverbund  durchge- 
ftihrt.  Tab.  4-15  gibt  die  Lange  der  Gleitreibungszone,  im  folgenden  als  Gleitlange  bezeich- 
net,  die  in  einem  kohlenstofffaserverstarkten  Kohlenstoff  in  Abhangigkeit  vom  Reibungsko¬ 
effizienten  p  zwischen  Faser  und  Matrix  auftreten,  wieder. 


Gleitlangen  fur  den  Einzelfaserverbund  [pm] 

Reibungskoeffizient  p 

1%  Matrixschrumpfen 

2%  Matrixschrumpfen 

0,2 

40,00 

40,00 

0,4 

40,00 

22,52 

0,5 

31,83 

17,00 

0,6 

26,35 

14,70 

0,8 

19,48 

10,87 

Tab.  4-15:  Gleitlange  fur  Einzelfaserverbund  mit  1%  und  2  %Matrixschrumpfen 
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Abb.4-21:  Gleitlange  fur  Einzelfaserverbund  in  Abhangigkeit  vom  Reibungskoeffizienten 

fur  1%  -  und  2%  -  Matrixschrumpfen 


Es  zeigt  sich  eine  deutliche  Abhangigkeit  der  Gleitlange  vom  Reibungskoeffizienten.  Die 
gestrichelte  Linie  zwischen  p  =  0,2  und  0,4  deutet  an,  dass  Gleiten  liber  die  gesamte  Modell- 
hohe  stattfindet.  Dieser  Bereich  unmittelbar  am  oberen  Modellrand  muss  fur  eine  weitere 
Betrachtung  ausgeklammert  werden. 


4.8.5  Einfluss  der  Reibung  auf  die  Gleitlange  bei  1%  und  2%-  Matrixschrumpfen 
beim  Quasimehrfaserverbund 


Fiir  den  Quasimehrfaserverbund  wird  eine  Verringerung  der  Gleitlangen  erwartet,  da  die  ra- 
diale  Aufweitung  durch  die  umliegende  Verbundstruktur  verhindert  wird.  Die  Berechnung 
fur  den  Quasimehrfaserverbund  ist  im  folgenden  dargestellt. 


Gleitlangen  fiir  den  Quasimehrfaserverbund  [|im] 

Reibungskoeffizient  |i 

1%  Matrixschrumpfen 

2%  Matrixschrumpfen 

0,2 

40,00 

38,97 

0,4 

34,43 

20,96 

0,5 

27,39 

16,61 

0,6 

22,87 

13,74 

0,8 

16,70 

09,83 

Tab.4-16:  Gleitlangen  fur  Quasimehrfaserverbund  mit  1%-  und  2%-Matrixschrumpfen 
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Tab.  4-16  gibt  die  Gleitlangen  zwischen  Faser  und  Matrix  in  Abhangigkeit  vom  Reibungs- 
koeffizienten  (I  an.  Hierbei  wurde  wiederum  der  Prozentsatz  des  Matrixschrumpfens  von  1% 
und  2%  angenommen,  um  den  Einfluss  herstellungsbedingter  Schrumpfspannungen  nach  der 
Karbonisierung  bei  ca.  1000°C  darzustellen.  Bei  Betrachtung  der  Gleitlangen  gemaB  Abb.  4- 
22  des  Quasimehrfaserverbundes  sieht  man  deutlich  eine  Reduzierung  der  Gleitlangen  infol- 
ge  der  verhinderten  radialen  Aufweitung  und  damit  der  groBeren  radialen  Druckspannung. 


Abb  .4-22:  Gleitlange  fiir  Quasimehrfaserverbund  in  Abhangigkeit  vom  Reibungskoeffizien- 

ten  fur  1%-  und  2%-  Matrix schrumpfen 

Im  folgenden  sol!  nun  der  Einfluss  der  Matrixschrumpfspannungen  auf  die  Energiefreiset- 
zungsraten  bei  Penetration  (Faserbruch)  ermittelt  werden. 

4.8.6  Einfluss  der  Reibung  auf  die  Energiefreisetzungsrate  fiir  Penetration  G„ 

Nach  dem  Auftreten  von  Rissabzweigung  an  der  Faseroberflache  ist  es  durch  die  Span- 
nungsiibertragung  als  Folge  von  Reibung  trotzdem  moglich,  dass  Faserversagen  (Penetrati¬ 
on)  auftritt.  Dieses  mogliche  Verhalten  wird  nun  untersucht.  Fiir  die  Ermittlung  der  Energie¬ 
freisetzungsrate  wurde  die  Methode  nach  Kaninen  [30]  verwendet.  Hierbei  wurde  jeweils  die 
Last  AF  ermittelt,  die  bendtigt  wird,  um  den  Riss  (mit  der  Rissoffnung  Az)  wieder  zu  schlie- 
Ben. 
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Der  Rissfortschritt  Aa  wurde  hierbei  jeweils  vom  Interface  aus  in  die  Faser  in  der  gleichen 
Ebene  des  Anfangsrisses  angenommen.  Die  Gegeniiberstellung  von  Penetration  und  Ab- 
zweigung  ist  nicht  sinnvoll,  da  Abzweigung  (durch  die  fehlende  Haftung)  hierbei  nicht  rele¬ 
vant  ist.  Meine  weiteren  Betrachtungen  beschranken  sich  daher  auf  den  Rissfortschritt  als 
Pentrationsriss.  Die  Energiefreisetzungsrate  G  wird  nach  folgender  Formel  fur  das  FEM- 
Modell  berechnet: 


G  =  —  FAu  (4-25) 

2A A 

wobei  AA  die  Flache  der  konzentrischen  Risserweiterung  und  F  die  Knotenkraft  darstellen, 
die  notwendig  ist,  um  die  Rissoffnung  Au  auf  „NULL“  zu  kompensieren.  Die  Einzelwerte 
wurden  fur  die  jeweiligen  Reibungskoeffizienten  aus  jeweils  zwei  Rechnungen  ermittelt. 
Nach  der  ersten  Rechnung  (Ergebnis:  F)  wurde  der  Riss  um  eine  gewisse  Lange  Aa  erweitert, 
was  einer  konzentrischen  Erweiterung  AA  bei  rotationssymmetrischer  Anordnung  entspricht. 
Damit  wurde  die  Rissoffnung  Au  berechnet. 


4.8.7  Ermittlung  von  Gn  in  Abhangigkeit  vom  Reibungskoeffizienten  fur  den  Einzel- 
faserverbund 


Bei  der  folgenden  Rechnung  wird  die  Abhangigkeit  der  Energiefreisetzungsrate  Gp  vom 
Reibungskoeffizienten  und  zusatzlich  der  Effekt  der  Matrixschrumpfspannungen  beriicksich- 
tigt.  Matrixschrumpfen  von  1%  bzw.  2%  wurde  beriicksichtigt.  Die  Ergebnisse  sind  fiber 
dem  Reibungskoeffizienten  aufgetragen. 


Energiefreisetzungsrate  Gp  [J/nf  ] 

Reibungskoeffizient  (1  [  1  ] 

1%  Matrixschrumpfen 

2%  Matrixschrumpfen 

0,2 

0,078 

0,162 

0,4 

0,099 

0,322 

0,5 

0,112 

0,385 

0,6 

0,125 

0,469 

0,8 

0,148 

0,641 

Tab. 4-1 7:  Energiefreisetzungsraten  fur  Einzelfaserverbund  mit  1%  und  2%- 

Matrixschrump  fen 
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Abb.4-23:  Energiefreisetzungsraten  in  Abhangigkeit  von  Li  fur  den  Einzelfaserverbund  bei 

1%-  und  2%-Matrixschrumpfen 

Abb.4-23  zeigt  eine  deutliche  Zunahme  der  Energiefreisetzungsrate  bei  VergroBerung  des 
Reibungsfaktors  zwischen  Faser  und  Matrix.  Wenn  man  mehr  Reibung  zulasst,  wird  folglich 
die  Wahrscheinlichkeit  eines  Faserbruches  erhoht.  Die  Energiefreisetzungsrate  steigt  bei 
groBeren  Matrixschrumpfspannungen  ebenfalls  deutlich  an.  Die  Tendenz  zum  Faserbruch 
bei  hohen  Matrixschrumpfspannungen  ist  somit  eher  gegeben.  Vergleicht  man  die  Ergebnis- 
se  mit  den  Ergebnissen  aus  den  Berechnungen  mit  perfekter  Haftung,  so  ist  ersichtlich,  dass 
die  Energiefreisetzungsraten  deutlich  geringer  sind.  Dies  ist  damit  zu  erklaren,  dass  nur  noch 
ein  geringer  Anteil  der  Spannungen  im  Interface  infolge  der  Reibung  iibertragen  wird  und 
damit  nur  noch  eine  geringe  Energiefreisetzungsrate  fur  Risspenetration  resultiert. 

Um  Aussagen  fur  den  Quasimehrfaserverbund  zu  erhalten,  wurden  die  Rechnungen  auch  fur 
diesen  durchgefiihrt.  Die  Ergebnisse  sind  im  folgenden  dargestellt. 

4.8.8  Ermittlung  der  Abhangigkeit  der  Energiefreisetzungsrate  von  den 
Reibungskoeffizienten  fur  den  Quasimehrfaserverbund 

Im  folgenden  wird  nun  der  Einfluss  von  Matrixschrumpfung  beim  Quasimehrfaserverbund 
diskutiert. 
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Es  wurde  wie  bereits  beschrieben  die  Energiefreisetzungsrate  fur  den  Spannungszustand  im 
Mehrfaserverbund  ermittelt,  wo  die  Rissflache  im  Vergleich  zum  Verbundquerschnitt  ver- 
nachlassigbar  klein  ist.  Dies  bedeutet,  dass  die  Rissebene  die  gleiche  radiale  Verformung  er- 
fahrt  wie  eine  Ebene  weit  von  dem  Riss  entfemt.  Fur  die  Simulation  des  Quasimehrfaserver- 
bundes  mussten  daher  wieder  pro  Energiefreisetzungsrate  eine  Vorrechnung  erfolgen,  mit 
der  die  Verschiebung  am  oberen  Modellrand,  d.h.  weit  weg  vom  Riss  ermittelt  wurde.  Diese 
radiale  Verschiebung  wurde  dann  auf  die  gesamte  AuBenflache  aufgepragt. 


Energiefreisetzungsrate  Gp  [J/m2] 

Reibungskoeffizient  p  [1] 

1%  Matrixschrumpfen 

2%  Matrixschrumpfen 

0,2 

0,076 

0,605 

0,4 

0,109 

1,273 

0,5 

0,145 

1,723 

0,6 

0,169 

2,218 

0,8 

0,206 

3,563 

Tab.4-18:  Energiefreisetzungsraten  fur  Quasimehrfaserverbund  mit  1%  und  2%  Matrix- 

schrumpfen 


Abb.4-24:  Energiefreisetzungsraten  in  Abhangigkeit  von  |i  fur  den  Quasimehrfaserverbund 
bei  1%  Matrixschrumpfen  und  2%  Matrixschrumpfen 
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Es  zeigt  sich,  dass  die  Energiefreisetzungsraten  fiber  denen  des  Einzelfaserverbundes  liegen. 
Diese  Tendenz  wurde  schon  bei  dem  Modell  mit  perfekter  Haftung  ermittelt.  Die  fiir  den 
Einzelfaserverbund  ermittelte  Tendenz,  namlich  den  Anstieg  der  Energiefreisetzungsrate  mit 
der  Zunahme  des  Matrixschrumpfens,  wird  auch  hier  wieder  bestatigt.  Die  Steigung  der 
Kurven  mit  2%  Matrixschrumpfung  liegt  deutlich  fiber  der  1%-Kurve,  so  dass  auch  noch  die 
Abhangigkeit  vom  Reibungskoeffizienten  vergroGert  wird. 

Die  geringen  Werte  fur  die  Energiefreisetzungsrate  bei  1%  sind  darauf  zuriickzufuhren,  dass 
die  Radialkomponente  der  Druckspannungen,  die  auf  das  Interface  wirkt,  geringer  ist  als  bei 
2%-Matrixschrumpfung.  Damit  kann  nur  ein  geringe  Kraftiibertragung  zwischen  Faser  und 
Matrix  erfolgen,  was  sich  in  einer  niedrigen  Energiefreisetzungsrate  auBert. 

4.8.9  Ergebnisse  der  Untersuchung  des  Einflusses  der  Reibung  zwischen  Faser  und 
Matrix  auf  die  Eigenschaften  von  kohlenstofffaserverstarktem  Kohlenstoff 

Ergebnisse  der  Energiefreisetzungsraten  bei  1%  und  2%  Matrixschrumpfen 

Die  Darstellung  der  Energiefreisetzungsraten  bei  Penetration  Gp  gemafi  Abb.4-23  zeigt  eine 
lineare  Abhangigkeit  zwischen  Gp  und  dem  Reibungskoeffizienten  |_l.  Mit  zunehmendem 
Reibungskoeffizienten  nimmt  die  Energiefreisetzungsrate  ebenfalls  zu.  Charakteristisch  ist 
auch,  dass  die  Energiefreisetzungsrate  fur  2%  Matrixschrumpfen  deutlich  hoher  ist  als  fur 
1%  Matrixschrumpfen.  Bei  hohen  Schrumpfspannungen  wird  somit  ein  Faserbruch  eher  er- 
moglicht  als  bei  niedrigen.  Ein  groBer  Reibungskoeffizient  zwischen  Faser  und  Matrix  be- 
wirkt  ebenfalls  eine  Tendenz  zum  Faserbruch.  Fiir  den  Quasimehrfaserverbund  zeigt  sich 
ebenfalls  eine  lineare  Abhangigkeit  in  dem  Bereich  zwischen  ji=0,4  und  |i=0,8.  Die  Charak- 
teristik  der  Verlaufe  entsprechen  denen  des  Einzelfaserverbundes.  Lediglich  die  Hohe  der 
Energiefreisetzungsraten  fibertrifft  die  des  Einzelfaserverbundes. 

Betrachtung  der  Spannungen  in  Abhangigkeit  vom  Matrixschrumpfen  und  Reibung 
im  Interface 


Es  werden  nun  die  Auswirkungen  der  Anwesenheit  von  Schrumpfspannungen  in  der  Matrix 
bzw.  Veranderung  des  Reibungskoeffizienten  auf  die  Radial  spannungen  und  Axialspannun- 
gen  im  Interfacebereich  dargestellt. 
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Radialspannungen 

Zunachst  soli  der  Einzelfaserverbund  betrachtet  werden.  In  Abb.  4-24  ist  die  Radialspannung 
dargestellt.  Die  Kurven  geben  die  Abhangigkeit  der  Radialspannung  entlang  des  Interfaces 
in  Abhangigkeit  vom  Reibungskoeffizienten  an.  Charakteristisch  ist  zunachst,  dass  die  Radi¬ 
alspannung  in  einer  diskreten  Entfemung  vom  Riss  im  Druckbereich  ist. 


Abb.4-25:  Radialspannungen  im  Interface 

Weiterhin  besitzen  die  Spannungsverlaufe  fur  p=0,6  und  0,8  einen  charakteristischen  Ver- 
lauf.  Die  Gleitlange  wird  namlich  durch  ein  lokales  Spannungsminiinum  am  Ubergang 
von  der  Gleitreibungszone  zur  Haftreibungszone  gekennzeichnet. 

Fur  p=0,6  betragt  die  Gleitlange  gemaB  Tab.  4-15  fur  1%  Matrixschrumpfen  26,35pm.  Fiir 
2%  Matrixschrumpfen  wird  sie  auf  14,7pm  reduziert.  Fiir  p=0,8  folgt  eine  Gleitlange  von 
19,48pm  bzw.  10,87pm.  Weiter  erkennt  man,  dass  eine  Verdoppelung  des  Matrixschrump- 
fens  eine  deutliche  Vergrofierung  der  Radialspannungen  in  den  Druckbereich  hervorruft.  Es 
zeigt  sich  auch,  dass  im  Gleitbereich  (innerhalb  der  Gleitlange)  eine  Abhangigkeit  der  Radi¬ 
alspannung  vom  Reibungskoeffizienten  vorhanden  ist.  Kleine  Reibungskoeffizienten  fiihren 
zu  einer  Verschiebung  der  Spannung  in  Richtung  des  Zugbereiches  (also  zu  kleineren 
Druckspannungen  hin).  Daraus  erkennt  man,  dass  ein  Interface  mit  geringer  Reibung  zwi- 
schen  Faser  und  Matrix  dem  Matrixriss  ein  Fortschreiten  im  Interfacebereich  ermoglicht  und 
die  Wahrscheinlichkeit  eines  Faserbruches  verringert. 
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Axial  spannun  gen.  Faser 

Die  Betrachtung  der  Abb.4-26  gibt  einen  Eindruck  von  den  Axialspannungen  in  der  Faser 
entlang  des  Interfaces.  Auch  hieraus  lassen  sich  wieder  anhand  der  lokalen  Minima  die 
Gleitlangen  ablesen. 


Der  Vergleich  zwischen  den  Axialspannungsverlaufen  im  unmittelbaren  Bereich  des  Riss- 
grundes  ist  hier  von  groBer  Bedeutung.  Die  Berechnung  mit  einer  Matrixschrumpfspannung 
von  2%  und  einem  Reibungskoeffizienten  von  0,8  ergibt  hierbei  die  groBte  Zugspannung  an 
der  Rissspitze  und  damit  der  groBten  Spannungssingularitat.  Rechnungen  mit  geringer  Mat- 
rixschrumpfung  und  kleinen  Reibungskoeffizienten  fuhren  hier  zu  kleinen  Werten. 
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4.8.10  Zusammenfassung  der  Rechnungen  mit  Reibung  zwischen  Faser  und  Matrix 


Die  Anbindung  im  Faser-/Matrixbereich  als  perfekte  Haftung  bewirkt  eine  deutliche  Ten- 
denz  zum  Faserversagen.  Die  Moglichkeit,  im  Interfacebereich  eine  Kraftubertragung  durch 
Reibung  zu  ermoglichen,  verringert  die  Wahrscheinlichkeit  des  Faserbruches.  Hierbei  wird 
eher  die  Tendenz  zum  Debonding  und  damit  Abzweigen  des  Matrixrisses  und  Ubergang  in 
einen  Riss  entlang  des  Interfaces  gegeben. 

Im  folgenden  wird  nun  eine  weitere  herstellungsbedingte  Eigenschaft  des  kohlenstofffaser- 
verstarkten  Kohlenstoffs  untersucht.  Neben  den  hier  betrachteten  Matrixschrumpfspannun- 
gen  nach  der  Karbonisierung,  wird  zusatzlich  eine  Graphitierung  der  Matrix  fur  moglich 
gehalten.  Die  Graphitierung  der  Matrix  soli  daher  in  den  folgenden  Seiten  naher  betrachtet 
werden.  Da  die  Matrixorthotropie  nur  schwer  quantifizierbar  ist,  wird  ein  bestimmter  Be- 
reich  untersucht.  Auch  hierbei  spielt  die  Auswirkung  auf  den  Rissfortschritt  wieder  eine  gro- 
Be  Rolle. 
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4.9  Auswirkungen  von  Strukturanderungen  in  der  Matrix  (Graphitierung)  auf 
Gh/Gp  nach  einer  zusatzlichen  Warmebehandlung 

Bisher  wurden  nur  Einflusse,  wie  z.B.  die  Kraftubertragung  durch  Reibung  im  Interface  bzw. 
die  herstellungsbedingten  Matrixspannungen,  die  aufgrund  des  Matrixschrumpfens  entste- 
hen,  beriicksichtigt.  Eine  weitere  Ursache  fur  eine  Anderung  des  Rissfortschrittsverhalten 
kann  in  einer  Strukturanderung  der  Matrix  liegen.  Bisher  wurde  davon  ausgegangen,  dass  die 
Matrix  nach  der  zusatzlichen  Warmebehandlung  isotrop  bleibt.  Dies  ist  aber  in  der  Realitat 
nicht  der  Fall.  Es  ist  davon  auszugehen,  dass  die  Matrix  nach  der  zusatzlichen  Warmebe¬ 
handlung  auch  eine  gewisse  Orientierung  besitzt.  Eine  Graphitierung  der  Matrix  tritt  nach 
bisherigen  Auswertungen  von  Aufnahmen  mit  polarisiertem  Licht  nur  lokal  auf.  Um  eine 
Aussage  zu  erhalten,  inwieweit  die  Orthotropie  der  Matrix  einen  Einfluss  hat,  wurde  die 
Matrix  in  den  folgenden  Berechnungen  komplett  orthotrop  angenommen. 

4.9.1  Einfluss  einer  moglichen  Eigenschaftsanderung  in  der  Matrix 

Fur  diese  Berechnung  wurde  eine  Orthotropie  der  Matrix  beriicksichtigt.  Als  Grundvoraus- 
setzung  wurde  wiederum  perfekte  Haftung  angenommen.  Der  Gleitmodul  wurde  variiert.  Fur 
p.zr  wurden  Werte  zwischen  10000  MPa  und  18000  MPa  in  2000  MPa-Schritten  gewahlt.  Das 
Orthotropieverhaltnis  [Steifigkeit  in  Faserrichtung/Steifigkeit  senkrecht  zur  Faserrichtung] 
wurde  mit  1/2,  2/3,  3/2  und  2/1  betrachtet.  Dabei  wurde  einmal  die  Steifigkeit  in  Faserrich¬ 
tung  (also  in  z-Richtung)  und  einmal  die  Steifigkeit  in  radialer  Richtung  (also  senkrecht  zur 
Faser)  konstant  gehalten.  Als  Belastung  wurde  die  Bruchspannung  von  888  MPa  eingeleitet. 
Thermische  Spannungen  wurden  durch  das  Abkiihlen  von  1 800°C  auf  20°C  beriicksichtigt. 
Als  Ergebnis  wurde  das  Verhaltnis  Gd/Gp  errechnet.  Weiterhin  wurde  der  Gp-Wert,  also  die 
Energiefreisetzungsrate,  getrennt  betrachtet.  Hierbei  wurde  eine  Normierung  auf  die  gemafi 
Honjo  und  Shindo  [8]  bzw.  Ochiai  und  Osamura  [9]  ermittelte  Bruchenergie  der  Faser  Gcf 
vorgenommen. 

4.9.2  Vorgehensweise 

Um  einen  Gd/Gp-Wert  zu  erhalten,  miissen  jeweils  drei  Rechnungen  erstellt  werden.  In  der 
ersten  Ausgangsrechnung  wird  Risserweiterung  in  der  Matrix  errechnet.  Die  Bifurkation 
lasst  wieder  zwei  mogliche  Arten  des  Rissfortschrittes  zu.  Eine  zweite  Rechnung  simuliert 
folglich  den  Rissfortschritt  in  die  Faser  [Penetration].  Eine  dritte  Rechnung  dient  zur  Be- 
stimmung  der  Energiefreisetzungsrate  fur  Rissabzweigung  am  Interface  [Deflection]. 


119 


Rissausbreitung  in  kohlenstofffaserverstarktem  Kohlenstoff  (C/C) 
Thorsten  Weber  _ 


DLR  -  Institut  fur  Werkstoff-Forschung 


Im  folgenden  wurden  jeweils  vier  verschiedene  Werte  fur  das  Orthotropieverhaltnis  angege- 
ben.  Der  Gleitmodul  variiert  monoton  wachsend  mit  5  Werten.  Es  werden  jeweils  zwei  ge- 
trennte  Parameteruntersuchungen  fur  die  Matrixorthotropie  durchgefuhrt. 

Eine  Parameteruntersuchung  halt  den  E-Modul  Ezm  der  Matrix  konstant  und  variiert  den  E- 
Modul  senkrecht  dazu,  also  Em.  Die  zweite  Parameteruntersuchung  variiert  den  E-Modul  der 
Matrix  in  axialer  Richtung  Ezm  und  halt  Em,  konstant.  Fur  die  u.a.  Diagramme  waren  somit 
2x5x4x3=120  Berechnungen  notwendig. 

4.9.3  F.„„=const 

Fur  die  Berechnung  mit  konstantem  Matrixelastizitatsmodul  in  axialer  Richtung  Ezm  wird 
zum  einen  der  Einfluss  auf  das  Verhaltnis  der  Energiefreisetzungsraten  Gd/Gp  untersucht  und 
zum  anderen  die  Auswirkung  auf  die  Energiefreisetzungsrate  bei  Penetration  Gp  in  Bezug 
auf  die  Faserbruchenergie  betrachtet. 

4.9.3.1  Untersuchung  der  Verhaltnisse  der  Energiefreisetzungsraten 

Abb.  4-27  zeigt  die  dreidimensionale  Darstellung  der  Abhangigkeit  des  Energiefreisetzungs- 
ratenverhaltnisses  vom  Schubmodul  der  Matrix  und  des  Orthotropieverhaltnisses  (Verhaltnis 
E-Modul  in  axialer  Richtung  zum  E-Modul  in  radialer  Richtung).  Eine  Anderung  der  Mat- 
rixeigenschaft  hinsichtlich  einer  VergroGerung  der  Schubsteifigkeit  und  radialer  Steifigkeit 
bewirkt  eine  erhebliche  Verringerung  des  Energiefreisetzungsratenverhaltnisses  und  damit 
die  Tendenz  zum  Faserbruch. 

Ein  kleiner  Schubmodul  und  eine  Matrix,  die  in  axialer  Richtung  steifer  ist,  lassen  das  Ver¬ 
haltnis  Gd/Gp  stark  anwachsen.  Die  Tendenz  zu  einer  Rissabzweigung  ware  somit  bei  An- 
nahme  einer  solchen  Anderung  der  Matrixeigenschaft  groBer. 

Abb.  4-28  zeigt,  dass  mit  zunehmender  Schubsteifigkeit  das  Verhaltnis  Gd/Gp  abnimmt.  Eine 
VergroGerung  der  Schubsteifigkeit  (VergroGerung  des  Schubmoduls  der  Matrix)  um  80% 
bewirkt  einer  Verringerung  des  Verhaltnisses  der  Energiefreisetzungsraten  Gd/Gp  von  10,9  % 
bei  einem  Orthotropieverhaltnis  von  Ezm/E1TO=0,5.  Eine  Abnahme  des  Verhaltnisses  der  E- 
nergiefreisetzungsraten  hat  eine  Tendenz  zum  Faserbruch  zur  Folge.  Das  o.a.  Diagramm 
zeigt  jedoch,  dass  sehr  groBe  Anderungen  in  der  Schubsteifigkeit  produziert  werden  miissen, 
um  eine  markante  Anderung  des  Energiefreisetzungsratenverhaltnisses  durchzufuhren. 
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Orthotropieverhaltnis  E?m/En71  [1] 


Abb.  4-29:  Verhaltnis  Gd/Gp  in  Abhangigkeit  von  der  Orthotropie  in  der  Matrix 

Eine  abnehmende  Steifigkeit  in  radialer  Richtung  in  der  Matrix  (Erm  wird  kleiner)  vergroBert 
die  Wahrscheinlichkeit,  dass  der  Riss  ins  Interface  abzweigt,  da  Gd/Gp  groBer  wird.  Eine 
Vervierfachung  des  Orthotropieverhaltnisse  (Ezm/Em  mit  Ezm=const)  fuhrt  zu  einem  Anstieg 
von  23,1%  im  Verhaltnis  der  Energiefreisetzungsraten  gemaB  Abb.4-29. 

4.93.2  Auswirkung  auf  die  Energiefreisetzungsraten  bei  Penetration 

Betrachtet  man  die  Energiefreisetzungsrate  fur  Penetration  in  Relation  auf  die  Faserbruch- 
energie  absolut,  so  zeigt  sich,  dass  mit  zunehmender  Schubsteifigkeit  Penetration  und  damit 
Faserbruch  erzeugt  wird. 

Aus  Abb.  4-30  erkennt  man,  dass  eine  Erhohung  der  Schubsteifigkeit  der  Matrix  schnell  zu 
einer  Steigerung  der  Energiefreisetzungsrate  Gp  ffihrt. 

Absolut  gesehen  wirkt  hier  ein  groBer  Orthotropiefaktor  (also  kleines  EJ  noch  zusatzlich 
fur  eine  groBere  Neigung  zur  Penetration. 
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O  1.02 


Schubmodul  [MPa] 

Abb.  4-30:  Verhaltnis  Gp/Gcf  in  Abhangigkeit  vom  Schubmodul  der  Matrix 

Erm=const 

Die  Berechnung  unter  der  Annahme,  dass  der  Elastizitatsmodul  der  Matrix  in  radialer  Rich 
tung  konstant  bleibt  dient  als  Referenzrechnung. 


SillsStl 


Schubmodul  [MPa] 


Abb.  4-31:  Verhaltnis  Gd/Gp  in  Abhangigkeit  von  den  Eigenschaften  der  Matrix 
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Abb.  4-32:  Verhaltnis  Gp/Gcf  in  Abhangigkeit  von  den  Eigenschaften  der  Matrix 

4.9.5  Ergebnisse  aus  der  Matrixorthotropie 

Die  Auswertung  der  Abb.  4-27  bis  4-32  gibt  einen  nachhaltigen  Eindruck,  wie  sich  eine 
Graphitierung  der  Matrix  auf  das  Bruchverhalten  des  kohlenstofffaserverstarkten  Kohlen¬ 
stoff  auswirkt.  Die  zusatzliche  Warmebehandlung  nach  der  Karbonisierung  beeinflusst  daher 
das  Werkstoffverhalten  nicht  nur  durch  die  Erholung  der  herstellungsbedingten  Schrumpf- 
spannungen,  sondern  auch  durch  eine  Eigenschaftsanderung  der  Matrix.  Die  Ergebnisse 
werden  im  folgenden  dargestellt. 

E7m=const 

Anhand  der  Abb.  4-27  bis  4-30  zeigt  sich  folgende  Grundtendenz  ab: 

•  Wenn  die  Matrixeigenschaften  hin  zu  einer  groBeren  Schubsteifigkeit  [groBer  Gleit- 
modul]  tendieren,  so  muss  bei  kleinen  Orthotropieverhaltnissen  [Steifigkeit  in  axialer 
Richtung  kleiner  als  in  radialer  Richtung]  mit  einer  deutlichen  Tendenz  zum  Faser- 
bruch  gerechnet  werden. 

•  Dahingegen  wird  ein  Rissabzweigen  des  Matrixrisses  in  den  Interfacebereich  eher 
begiinstigt,  wenn  die  Schubsteifigkeit  der  Matrix  sehr  gering  ist  und  zusatzlich  das  Or- 
thotropieverhaltnis  (Steifigkeit  in  radialer  Richtung  sehr  viel  kleiner  als  in  axialer 
Richtung)  sehr  groB  ist. 
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Im  folgenden  sind  nun  noch  einmal  die  Einfliisse  der  Schubsteifigkeit  und  der  Orthotropie 
einzeln  aufgefuhrt. 


Orthotropiefaktor 

Schubmodul 

Wirkung 

kleiner 

groB 

Faserbruch 

groBer 

klein 

Interfacebruch 

Tab.4-19  Ergebnisse  aus  der  Orthotropiebetrachtung 


Neben  den  Ergebnissen  des  Energiefreisetzungsratenverhaltnisses  Gd/Gp  bzw.  des  Verhalt- 
nisses  aus  Energiefreisetzungsrate  bei  Penetration  zur  der  Bruchenergie  der  Faser  sind  auch 
die  auftretenden  Radialspannungen  von  groBem  Interesse.  Diese  werden  im  Anschluss  dar- 
gestellt. 

Auswirkung  auf  die  Spannungen  im  Interfacebereich  bzw.  an  der  Rissspitze  fur 
Ezm^const 

Die  Anderung  von  in  Abhangigkeit  von  der  Matrixorthotropie  ist  in  Abb.  4-33  dargestellt. 
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Abb.  4-33:  c„  in  Abhangigkeit  von  der  Matrixorthotropie 
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Offenbar  bewirkt  eine  VergroBerung  des  Orthotropiefaktors  infolge  eines  kleinen  Matrix-E- 
Modul  in  radialer  Richtung  eine  Erhohung  der  Radialspannungen  langs  des  Interfaces.  Die 
Druckbeanspruchung  im  Interface  wird  dadurch  reduziert.  Es  besteht  fur  Matrixrisse  nun 
die  Chance  sich  im  Interface  auszubreiten.  Dagegen  fuhrt  eine  Verringerung  des  Orthotro¬ 
piefaktors  zu  groBen  Druckspannungen  im  Interfacebereich,  was  zu  einer  Behinderung  von 
Interfacerissen  fuhrt  und  Penetrationsrisse  (Sprodbruch)  begiinstigt. 

Erm=const 


Die  Ergebnisse  fur  E^const  geben  einen  Uberblick  iiber  das  Verhalten  des  kohlenstofffa- 
serverstarkten  Kohlenstoffs  unter  der  Annahme,  dass  die  Matrixeigenschaften  nach  der  zu¬ 
satzlichen  Warmebehandlung  nur  in  axialer  Richtung  sich  geandert  haben  und  in  radialer 
Richtung  annahemd  konstant  geblieben  sind. 

Abb.  4-31  zeigt,  dass  groBe  Schubsteifigkeit  der  Matrix  zu  einem  kleinen  Gd/Gp  fuhrt  und 
damit  zu  sprodem  Verhalten,  da  der  Riss  mit  hoher  Wahrscheinlichkeit  zum  Faserbruch 
fuhrt.  Faserbruch  und  damit  Schwachung  des  kohlenstofffaserverstarkten  Kohlenstoffs  zeigt 
sich  deutlich  an  der  o.a.  Abb.4-32. 

Man  erkennt  anhand  der  Normierung  auf  die  Bruchenergie  der  Faser  deutlich,  dass  jede  Ver- 
steifung  der  Matrix  in  axialer  Richtung  zu  einer  Uberschreitung  der  Bruchenergie  Gcf  fuhrt 
und  somit  einen  Sprodbruch  einleitet. 

4.9.6  Zusammenfassung  der  Ergebnisse  aus  der  Betrachtung  einer  Graphitierung  der 
Matrix  nach  der  zusatzlichen  Warmebehandlung 

Die  Ergebnisse  zeigen  deutlich,  dass  die  beobachtete  Graphitierung  der  Matrix  nach  der  zu¬ 
satzlichen  Warmebehandlung  einen  deutlichen  Effekt  auf  das  Werkstoffverhalten  des  koh¬ 
lenstofffaserverstarkten  Kohlenstoffs  hat. 

Neben  den  Matrixschrumpfspannungcn,  die  nach  der  Karbonisierung  auftreten  und  im  Laufe 
einer  zusatzlichen  Warmebehandlung  abklingen  und  der  Beschaffenheit  des  Interfaces 
(Kraftiibertragung  Reibung  zwischen  Faser  und  Matrix  oder  perfekte  Haftung)  spielt  eine 
Graphitierung  der  Matrix  bei  hohen  Temperaturen  eine  nicht  zu  vemachlassigende  Rolle  bei 
der  Untersuchung  des  Bruchverhaltens. 
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5.  Diskussion 


Das  am  Anfang  der  Arbeit  geschilderte  Verhalten  von  kohlenstofffaserverstarktem  Kohlen¬ 
stoff  wahrend  des  Herstellungsprozesses  wurde  anhand  von  analytischen  Ansatzen  betrachtet 
und  anschlieBend  mit  Hilfe  der  Methode  der  Finiten  Elemente  einer  genauen  Analyse 
unterzogen. 

Um  das  Verhalten  zu  analysieren,  wurden  umfangreiche  Berechnungen  der  Spannungen  im 
Interfacebereich  durchgefuhrt  und  die  Energiefreisetzungsrate  bei  Penetration  und  Rissab- 
zweigen  am  Interface  durchgefuhrt.  Anhand  von  Druckscherversuchen  [33,  34]  war  es  mog- 
lich  die  Bruchenergie  des  Interfaces  im  kohlenstofffaserverstarktem  Kohlenstoff  zu  ermit- 
teln.  Ausgehend  von  den  Untersuchungen  von  Ochiai,  Osamura,  Honjo  und  Shindo  [7,8]  an 
TiC-beschichteten  Kohlenstofffasem  konnte  mit  Hilfe  eines  iterativen  Verfahrens  die  Bruch¬ 
energie  der  Kohlenstofffaser  ermittelt  werden.  Im  folgenden  wird  nun  auf  Grundlage  der  in 
dieser  Arbeit  erbrachten  Ergebnisse  das  Verhalten  wahrend  des  Herstellungsprozesses  disku- 
tiert. 


5.1  Verhalten  des  kohlenstofffaserverstarkten  Kohlenstoffs  aufgrund  der 
herstellungsbedingten  Matrixschrumpfspannungen 


Charakteristisch  fur  die  Festigkeitsentwicklung  nach  den  verschiedenen  Schritten  bei  der 
Herstellung  von  kohlenstofffaserverstarktem  Kohlenstoff  ist  das  Auftreten  eines  Festig- 
keitsminimums  unmittelbar  nach  der  Karbonisierung  bei  ca.  1 000°C.  Die  auftretende  Bruch- 
flache  ist  dabei  glatt,  was  auf  ein  Eindringen  des  Matrixrisses  in  die  Faser  hindeutet,  die  dar- 
aufhin  sprode  versagt.  Die  zusatzliche  Warmebehandlung  oberhalb  1000°C  fuhrt  zu  einem 
anderen  Bruchbild.  Hierbei  treten  deutlich  Rissabzweigungsmerkmale  auf.  Die  Fasem  versa- 
gen  erst  nach  Auftreten  eines  Debondings  zwischen  Faser  und  Matrix.  Anhand  des  Energie- 
kriteriums  kann  nun  ermittelt  werden,  ob  das  Auftreten  von  Schrumpfspannungen  unmittel¬ 
bar  wahrend  der  Karbonisierung  hierfur  verantwortlich  gemacht  werden  kann.  Faserbruch 
tritt  gemaB  He  und  Hutchinson  dann  auf,  wenn  folgende  Ungleichung  erfullt  ist: 


G{ 


(5-1) 
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Die  bisher  in  der  Literatur  ermittelten  Werte  fur  das  Energiefreisetzungsratenverhaltnis 
Gd/Gp  fur  Werkstoffkombinationen  aus  isotropen  Komponenten  nach  He  und  Hutchinson  [5] 
wurden  unter  der  Annahme  eines  mittleren  Schubmoduls  fur  die  orthotrope  Faser  ein  Wert 
von  0,248  ermittelt.  Dieser  Wert  liefert  allerdings  nur  einen  Anhaltswert.  Die  Erweiterung 
von  Gupta  und  Martinez  [11,  22]  fur  orthotrope  Werkstoffkombinationen  lieferte  mit  0,118 
einen  wesentlich  besseren  Wert.  Nachteil  der  analytischen  Ansatze  ist  jedoch,  dass  keinerlei 
Moglichkeit  einer  Matrixschrumpfspannung  beriicksichtigt  wird.  Die  Finite  Element  Rech- 
nung  lieferte  daher  wesentlich  genauere  Werte. 

Fur  den  Quasimehrfaserverbund  lautet  das  Bruchkriterium  unmittelbar  nach  der  Karboni- 
sierung  bei  1000°C  unter  Einbeziehung  der  Energiefreisetzungsraten  Gd  mit  1,906  J/m2  und 
Gp  13,179  J/m2  aus  Tab.  4-13  und  der  Bruchenergien  fur  das  Interface  zu  16,3  J/m2  gemaB 
Abschnitt  4.4.2  und  fur  die  Faser  mit  15  J/m2  gemaB  4.3.6.  wie  folgt: 

1,906  4  16,3-4 

m  <  m 

13,79-4  15^t 

m  m 

0,1446  <  1,087 

Das  Bruchkriterium  zeigt  deutlich  die  Tendenz  zum  Faserbruch  unmittelbar  nach  der  Karbo- 
nisierung  an. 

Nach  der  zusatzlichen  Warmebehandlung  oberhalb  von  1000°C  folgt  das  Energiekriterium 
unter  Einbeziehung  der  entsprechenden  Energiefreisetzungsraten  (Gd/Gp  gemaB  Tab  4-14) 
und  Bruchenergien  (Gc'  gemaB  Abschnitt  4.4.2  und  Gcf  aus  Abschnitt  4.3.6): 

0,318^4  10,8^4 

_ m_  <  m 

2,514-4  15-4 

m  m 

0,1265  <0,72 

Ausgehend  von  diesem  Energiekriterium  wird  fur  den  zusatzlich  warmebehandelten  kohlen- 
stofffaserverstarkten  Kohlenstoff  ebenfalls  ein  Faserbruch  vorausgesagt.  Allerdings  zeigt 
sich  jedoch  deutlich  die  Tendenz  in  Richtung  Debonding  zwischen  Faser  und  Matrix  ab,  da 
das  Bruchenergieverhaltnis  sich  erheblich  verringert  hat. 
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Gd/Gp  [1] 

GcVGc'M 

He  und  Hutchinson  [5] 

0,248 

Gupta  und  Martinez  [1 1,  22] 

0,118 

C/C  (karbonisiert) 

0,1446 

1,087 

C/C  (zusatzlich  warmebehandelt) 

0,1265 

0,72 

Tab.5-1  Energiefreisetzungsraten  /  Bruchenergien 

Vergleicht  man  die  absolute  Bruchenergie  der  Faser  von  15  J/m2  mit  der  Energiefreiset- 
zungsrate  der  Faser  unmittelbar  nach  der  Karbonisierung  mit  13,79  J/m2  bzw.  nach  der  zu- 
satzlichen  Warmebehandlung  von  2,514  J/m2,  so  zeigt  sich  deutlich  der  Einfluss  der  Relaxa¬ 
tion  des  Materials.  Unter  der  Voraussetzung  von  1%  Matrixschrumpfen  erreicht  die  Energie- 
freisetzungsrate  mit  92%  annahemd  den  Wert  den  Bruchenergie  der  Faser,  was  einen  Faser- 
bruch  ermoglicht.  Die  zusatzliche  Warmebehandlung  oberhalb  von  1000°C  lasst  die  Energie- 
freisetzungsrate  auf  einen  Wert  von  16,8%  der  Faserbruchenergie  fallen,  was  einen  Faser- 
bruch  ausgelost  durch  Matrixrisse  ausschlieBt. 

Die  Betrachtung  der  Energiefreisetzungsraten  fur  Risseindringen  in  die  beim  Einzelfaser- 
verbund  unterstreicht  diese  Aussage  zusatzlich.  Fur  das  restspannungsbehaftete  Material 
(unmittelbar  nach  der  Karbonisierung  bei  1000°C)  wird  gemaB  Tab.  4-9  eine  Energiefreiset- 
zungsrate  fur  Penetration  von  15,18  J/m2  ermittelt.  Diese  liegt  eindeutig  iiber  der  Bruchener¬ 
gie  der  Faser.  Die  anschlieBende  Warmebehandlung  des  Einzelfaserverbundes  lasst  die  E- 
nergiefreisetzungsrate  auf  2,439  J/m2  fallen,  was  nur  noch  16,3  %  der  Faserbruchenergie  ent- 
spricht. 

Die  Auswertung  der  Spannungsverlaufe  entlang  des  Interfaces  fur  den  Einzelfaserverbund 
zeigt  ebenfalls  deutlich  den  Einfluss  des  Matrixschrumpfens  bzw.  der  Matrixrestspannungen 
auf  das  Bruch verhalten.  So  treten  gemaB  Abb.4-16  die  maximalen  Spannungen  (Zugbereich) 
in  axialer  Richtung  und  damit  in  Faserrichtung  fur  grofie  Matrixschrumpfspannungen  in  der 
Nahe  der  Matrixrissspitze  auf.  In  diesem  Bereich  klingen  die  Axialspannungen  mit  zuneh- 
mender  Warmebehandlung  und  damit  Relaxation  (0%  Matrixschrumpfen)  auf  50%  des  Ma- 
ximalwertes  (2%  Matrixschrumpfen)  ab.  Erst  im  weiteren  Verlauf  des  Interfaces  fallen  sie  in 
den  Druckbereich  ab.  Das  zusatzlich  warmcbehandelte  Material  liefert  die  geringsten  Axial¬ 
spannungen  und  damit  die  geringste  Belastung  der  Faser  in  der  Nahe  des  Matrixrissspitze. 
Im  weiteren  Verlauf  des  Interfaces  bleiben  die  Axialspannungen  im  Zugbereich. 
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Die  Betrachtung  von  Abb.  4-17  gibt  die  Debondingtendenz  des  Matrixrisses  bei  dem  rela- 
xierten,  zusatzlich  warmebehandelten  kohlenstofffaserverstarkten  Kohlenstoff  wieder.  Die 
Radialspannungen  im  Interface  befinden  sich  fur  das  relaxierte  Material  (0%  Matrix- 
schrumpfen)  im  Zugbereich.  Damit  wird  ein  Debonding  zwischen  Faser  und  Matrix  eingelei- 
tet.  Im  Gegensatz  dazu  wird  ein  moglicher  Rissfortschritt  zwischen  Faser  und  Matrix  fur  das 
karbonisierte  Material  durch  die  Anwesenheit  von  starken  Druckspannungen  verhindert. 

Fazit:  Die  durchgefuhrten  FEM-Rechnungen  konnten  somit  das  Vorhandensein  von  Matrix¬ 
schrumpfspannungen  als  eine  Erklarung  fur  das  ermittelte  Versagensverhalten  von  kohlen¬ 
stofffaserverstarktem  Kohlenstoff  gelten  lassen. 

5.2  Verhalten  des  kohlenstofffaserverstarkten  Kohlenstoffs  unter  Beriicksichtigung 
einer  Kraftiibertragung  im  Interface  durch  Reibung 

Unter  Beriicksichtigung  von  Reibung  zwischen  Faser  und  Matrix  [46,  47,  48,  49]  ergibt  sich 
fur  das  Verhalten  von  kohlenstofffaserverstarktem  Kohlenstoff  eine  Anderung  der  Bruchcha- 
rakteristik.  Ein  Faserbruch  wird  aufgrund  der  geringen  Energiefreisetzungsraten  Gp,  die  nur 
noch  ca.  23,8%  der  Bruchenergie  der  Faser  erreicht,  nahezu  unwahrscheinlich. 


Quasimehrfaserverbund 

Faserbruchenergie  Gc* 

1 5  J/m2 

15  J/m2 

Gp-Bereich  in  Abhangigkeit  von  |i  fur  1%  Schrumpfen 

0,078-0,148  J/m2 

0,076-0,2065  J/m2 

Gp-Bereich  in  Abhangigkeit  von  g  fur  2%  Schrumpfen 

0,162-0,641  J/m" 

Tab. 5-2:  Energiefreisetzungsraten  bei  reibungsbehaftetem  Interface 

Eine  Eigenschaftsanderung  des  Interfaces  von  der  perfekten  Haftung  zu  einer  Kraftiibertra- 
gung  durch  Reibung  im  Interface  fuhrt  somit  ebenfalls  zu  einem  Anderung  des  Bruchverhal- 
tens.  Selbst  hohe  Matrixschrumpfspannungen  aufgrund  eines  2%-Matrixschrumpfens  erzeu- 
gen  keine  kritischen  Energiefreisetzungsraten,  die  die  Faserbruchenergie  iiberschreiten  konn¬ 
ten.  Das  Zulassen  von  geringer  Reibung  (|i=0,8)  fuhrt  zu  einer  Anderung  des  Bruchmodus. 
Somit  kann  ein  quasiduktiles  Verhalten  des  Werkstoffes  dadurch  begiinstigt  werden,  dass 
man  Reibung  zwischen  Faser  und  Matrix  zulasst,  wodurch  der  Einfluss  der  herstellungsbe- 
dingten  Matrixschrumpfspannungen  nahezu  unerheblich  wird. 
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5.3  Verhalten  des  kohlenstofffaserverstarkten  Kohlenstoffs  unter  Beriicksichtigung 
einer  Graphitierung  der  Matrix  nach  der  zusatzlichen  Warmebehandlung  o- 
berhalb  von  1000°C 


Neben  den  Restspannungen,  die  aufgrund  des  herstellungsbedingten  Matrixschrumpfens  ein- 
treten,  und  der  Kraftiibertragung  im  Interface  (perfekte  Haftung  bzw.  Gleitreibungs-  und 
Haftreibungszone)  spielt  die  Eigenschaftsanderung  der  Matrix  eine  wesentliche  Rolle.  Diese 
Eigenschaftsanderung  der  Matrix,  die  sich  durch  die  zusatzliche  Warmebehandlung  oberhalb 
von  1000°C  in  einer  teilweisen  Graphitierung  der  ursprunglichen  Glaskohlenstoffmatrix  au- 
Bert,  beeinflusst  ebenfalls  das  Bruchverhalten  des  kohlenstofffaserverstarkten  Kohlenstoffs. 

Anhand  der  Ergebnisse  aus  4.9.5  und  den  zugehorigen  Abb.  4-27  und  4-30  zeigt  sich,  dass 
die  Rissabzweigung  am  Interface  begiinstigt  wird,  wenn  die  Steifigkeit  der  Matrix  in  radialer 
Richtung  stark  abnimmt. 

Eine  Graphitierung  der  Matrix,  die  sich  in  einer  Abnahme  des  Elastizitatsmoduls  in  radialer 
Richtung  auswirkt,  wird  somit  neben  den  bereits  diskutierten  Effekten  eine  Anderung  des 
Bruchmodus  in  Richtung  Rissabzweigung  begiinstigen. 

Diese  Tendenz  wird  zusatzlich  durch  eine  mogliche  Verringerung  der  Schubsteifigkeit  des 
Matrixwerkstoffes  verstarkt. 

Die  im  Interface  auftretenden  Radialspannungen  gemaB  Abb.  4-33  unterstreichen  diesen  Ef- 
fekt  zusatzlich.  Die  Abnahme  des  Elastizitatsmoduls  der  Matrix  in  radialer  Richtung  verrin- 
gert  die  Druckbeanspruchung  des  Interfaces  und  gibt  somit  einer  Rissausbreitung  im  Inter¬ 
face  gegeniiber  einem  Faserbruch  den  Vorzug. 
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6.  Zusammenfassung 

Kohlenstofffaserverstarkter  Kohlenstoff  wird  in  der  Luft-  und  Raumfahrttechnik  fiir  ther- 
misch  hochbeanspruchte  Bauteile  verwendet.  Die  bisher  dort  fiir  niedrige  Einsatztemperatu- 
ren  vorgesehene  Werkstoffe  der  Titan-  und  Aluminiumlegierungen  besitzen  ein  duktiles 
Verhalten,  d.h.  sie  sind  schadenstolerant  und  Versagen  nicht  „katastrophal“.  Man  versucht 
diese  Eigenschafit  auch  fur  Faserverbundwerkstoffen  mit  keramischer  und  sproder  Matrix  zu 
realisieren.  Eine  Moglichkeit  in  sproden  Faserverbundwerkstoffen  eine  Quasiduktilitat  zu 
verwirklichen  ist  das  Zulassen  von  mehrfachen  Matrixrissen,  die  zu  einem  allmahlichem 
Verlust  des  Beitrages  der  Matrix  an  der  Last  fiihren.  Voraussetzung  fur  einen  Erhalt  der  Fes- 
tigkeit  und  Risszahigkeit  ist  jedoch,  dass  diese  Matrixrisse  an  der  Faseroberflache  abzwei- 
gen. 

Die  vorliegende  Arbeit  ist  eine  Parameters tudie,  die  das  Ziel  verfolgt,  das  Phanomen  der 
Rissausbreitung  durch  Penetration  in  die  Fasem  und  durch  Rissabzweigung  an  der  Fa- 
ser/Matrix  Grenzflache  fiir  C/C  zu  beschreiben.  Dies  wird  fur  C/C  durchgefuhrt,  der  mit  Hil- 
fe  von  Polymer-Impragnation  und  Pyrolyse  hergestellt  wird.  Dieser  Werkstoff  ist  auBerst 
komplex,  da  Fasem,  Matrix  und  Faser/Matrix-Grenzflache  bei  der  Pyrolyse  und  nachfolgen- 
der  Graphitierung  starke  Verandemngen  erfahren. 

Wesentliche  Verandemngen,  deren  Einfluss  auf  das  Risserweiterungsverhalten  beschrieben 
wird,  sind  als  Folge  der  Karbonisierung  (bis  1000°C) 

die  Umwandlung  der  Polymermatrix  in  eine  Glaskohlenstoffmatrix 
einhergehend  mit  einer  starken  Schrumpfung  der  Matrix 
eine  eventuelle  Faserschadigung 
und  als  Folge  der  Graphitierung  oberhalb  1000°C: 

eine  teilweise  Graphitierung  der  Glaskohlenstoffmatrix 
eine  verstarkte  Graphitierung  der  C-Faser 
eine  mogliche  Erholung  der  Faserfestigkeit 

Die  Parameterstudie  versucht  einen  Zusammenhang  zwischen  diesen  Parametem  und  expe- 
rimentellen  Ergebnissen,  erzielt  am  C/C  nach  der  Karbonisierung  bis  1000°C  (Bruchart: 
sprode)  (Bild  2-4)  und  erzielt  nach  einer  weiteren  Behandlung  bis  1800°C  (Bruchart:  ge- 
mischt)  (Bild  2-4),  zu  erstellen. 

Das  Rissausbreitungsverhalten  dieser  zwei  unterschiedlich  hergestellten  kohlenstofffaserver- 
starkten  Kohlenstoffe  infolge  der  o.a.  Faktoren  wird  mittels  des  Versagenskriteriums  nach 
He  und  Hutchinson  [5]  beschrieben. 
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Fiir  dieses  Kriterium  ist  die  Kenntnis  der  WerkstoffkenngroBen  Gc'  (die  Bruchenergie  des  In¬ 
terfaces  )  und  Gcf  (die  Mode  I  Bruchenergie  der  Faser)  notwendig.  Die  Bruchenergie  des  In¬ 
terfaces  wurde  fur  das  karbonisierte  Material  und  das  zusatzlich  warmebehandelte  Material 
anhand  von  Druckscherversuchen  bestimmt.  Daraus  ergab  sich  Gc‘  mit  16,3  J/m2  nach  der 
Karbonisierung  und  10,8  J/m2  nach  der  Graphitierung.  Die  Bruchenergie  der  C-Faser  wurde 
auf  der  Grundlage  von  experimentellen  Ergebnissen  der  Festigkeit  von  TiC-beschichteten 
Fasem  [8]  mit  Hilfe  von  Einzelfaserverbundmodellen  unter  Beriicksichtigung  von  thermi- 
schen  Spannungen  ermittelt.  Gcf  wurde  daraus  zu  15  J/m2  bestimmt.  Alle  Berechnungen  wur- 
den  an  Einzelfaserverbunden  durchgefuhrt. 

Ein  Riss  in  der  Matrix  wird  nach  dem  Energiekriterium  nach  He  und  Hutchinson  [5]  beim 
Auftreffen  auf  das  Interface  abzweigen,  wenn 


G'c  GP 


(6-1) 


wobei  Gd  die  Energiefreisetzungsrate  fur  einen  abzweigenden  Riss  und  Gp  die  Energiefrei- 
setzungsrate  fur  einen  Matrixriss,  der  in  die  Faser  eindringt,  darstellen. 

Die  wichtigsten  Ergebnisse  sind: 

-  Das  Kriterium  sagt  fur  beide  Werkstoffe  Penetration  von  Matrixrissen  in  die  Faser  vor- 
aus  (Sprodbruch).  Am  ausgepragtesten  ist  dies  fur  das  karbonisierte  C/C. 

Die  Penetration  von  Matrixrissen  in  die  C-Fasem  des  karbonisierten  C/C  wird  im  we- 
sentlichen  durch  die  hohe  Schrumpfspannung  in  der  Matrix  verursacht. 

Der  geringe  Wert  fur  Gp  fur  das  graphitierte  C/C  (bei  der  Bruchspannung  von  888  MPa) 
deutet  darauf  hin,  dass  Versagen  dieses  Werkstoffes  nicht  durch  Matrixrissausbreitung, 
sondem  eher  durch  Faserbruch  verursacht  wird. 

Der  Vergleich  der  beiden  moglichen  Kraftubertragungsarten  im  Interface,  namlich  per- 
fekte  Haftung  und  reibungsbehaftetes  Interface,  zeigen  erheblichen  Einfluss  auf  das 
Rissfortschrittsverhalten.  Durch  die  Kraftubertragung  im  Interface  durch  Reibung  wird 
nur  noch  eine  geringe  Energiefreisetzungsrate  fur  Penetration  ermittelt. 

Die  Umwandlung  der  Matrix  von  der  urspriinglichen  isotropen  Glaskohlenstoffmatrix 
(nach  der  Karbonisierung)  zu  einer  teilweise  orthotropen,  graphitierten  Kohlenstoffmat- 
rix  (nach  der  zusatzlichen  Warmebehandlung)  beeinflusst  ebenfalls  die  Art  der  Rissaus¬ 
breitung.  Eine  Rissabzweigung  im  Interface  wird  begiinstigt,  wenn  die  Steifigkeit  der 
Matrix  in  radialer  Richtung  (senkrecht  zur  Faserachse)  abnimmt.  Eine  Verringerung  der 
Schubsteifigkeit  des  Matrixwerkstoffes  wirkt  hierbei  ebenfalls  gunstig  auf  die  Rissab¬ 
zweigung. 
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